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I INTRODUCTION 


The potential role of an occultation satellite in the National Weather 
Satellite system has been, investigated under NASA Contracts NAS 1-99^2 
and NAS 1-99^3= Detailed reports of the findings are contained in the 
two attached volumes*. "Orbit Determination for a Weather Occultation 
Satellite," which describes the outcome of computer simulations of orbit 
dynamics and data processing procedures that would be used with an occul- 
tation system, and "An Occultation Satellite System for Determining 
Pressure Levels in the Atmosphere," which describes the radio and atmos- 
pheric physics of the occultation system, reports on the outcome of a 
simulation using occultation system data together with satellite infrared 
sounding data to determine the atmospheric pressure reference, and describes 
a ground based test of the system in Hawaii and the implications for a 
satellite occultation system. Following is a brief discussion of these 
two documents and their .implications . The intent is not to summarize 
the reports, the reports already contain adequate summaries, but rather 
to relate the findings to each other and to the present needs of the 
U.S. weather program. 

The physical basis of a radio occultation satellite system is that 
radio waves travel more slowly in the atmosphere than they do in a vacuum. 
Additionally, if the atmospheric density varies across the path of the 
radio waves, the radio path will be bent. These two effects combine to 
make the effective radio path between two satellites longer if the path 
goes through the atmosphere. The amount of this path increase is a function 
of the atmospheric density and hence provides the opportunity for measure- 
ment of atmospheric variations responsible for the earth* s weather. 

The use of radio occultation to measure atmospheres is not new. The 
most precise measurements of the atmospheres of Mars, Venus, and now 
Jupiter have been obtained by radio occultations of Mariner and Pioneer 
spacecraft. The equipment used in these efforts has precision far 



greater than that required for useful information on the Earth's atmrj ir- 
phere, The physics of radio occultation and the equipment performance 
have been well established,, What remains to be answered is whether the 
information that can be provided on a continuous operational basis by an 
earth radio-occultation system is useful to the world's meteorologists 
and whether the occultation system is the least costly way of obtaining 
that information o 

The work described in the attached reports concludes that the infor- 
mation supplied by an occultation system would indeed be of great use 
to the meteorologists „ However, whether the occultation system is the 
least costly way of obtaining that information is still to be resolved. 
The authors feel that is is, But it is clear from reaction to radio 
occultation proposals, that the cost and performance of other systems 
designed to obtain equivalent information must be more clearly defined. 

It is hoped that the upcoming tests of global observation systems and 
weather prediction technique will help to better define the performance 
and cost of these alternati ves „ 
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There are many ways the radio occultation technique could be 
employed to obtain useful weather information. The particular system 
evaluated in the following reports reflects several years of inter- 
action between the authors and members of the meteorological community. 

The proposed system would provide only information that is very hard to 
obtain by other measurement techniques „ To do this, the system makes 
use of data readily available from already proven satellite sensors, 
from existing ground based observations, and from historical records. 

In other words, the proposed system is intended to be a good match to 
other components of the weather program. 

The occultation system consists of two satellites in the same orbit. 
One satellite, the ,s master fr satellite, would carry the main Instrumentation 
and d^ta recording equipment j in an operational system this satellite would 
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also carry other sensors involved in the overall weather observation system. 
A append, much smaller "slave" satellite, would be located in the same orbit 
but behind the "master" satellite sufficiently far that radio signals be- 
tween the two would intersect the atmosphere in the vicinity of the 5 OO 
millibar altitude. A transmitter in the "master" satellite would send a 
radio signal to the receiver in the "slave" satellite. The "slave" 
satellite would shift the signals in frequency to avoid interference and 
transmit them back to a raido receiver in the "master" satellite. At the 
"master", the radio path length would be recorded and relayed to the earth 
by the normal satellite telemetry system,, 

The system would use the orbit best suited to the other meteorological 
sensors, most likely a sun-synchronous polar orbit. The occultation data 
would provide continuous information in the subsatellite orbit track. 

Used in conjunction with other information, the occultation system 
would provide weather prediction programs with continuous measurements of 
the world wide 5 OO millibar pressure reference level with better than 
24 meter accuracy. To achieve this, information required from other sources 
is continuous atmospheric temperature profiles obtained from satellite 
infrared sounding instruments, historic seasonal averages of lower atmos- 
pheric water vapor content, and a few daily surface or upper air pressure 
level measurements wherever it is most convenient to obtain them. 

The major concerns that have been raised about the proposed occulta- 
tion system are: 

•When the inaccuracies of infrared soundings are combined with the 
expected occultation inaccuracies and the complexities of inversion 
procedures, will the results have sufficient accuracy for meteorologi- 
cal applications? 

•Will the effects of water vapor in the lower atmosphere cause serious 
errors in the interpretation of the occultation measurements? 

•Will complex radio ray paths from inhomogenieties in the lower 
atmosphere cause loss of data or distortion of results? 

•Can the satellites motions in orbit due to inhomogenieties in the 
earth's gravity field, air drag, and satellite thrusters, be 
separated from changes in apparent distance caused by atmospheric 
pressure variations? 
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Analyses of these questions are presented in great detail in the 

attached reports. The conclusions are summarized briefly below; 

Inversion Accuracy 

•To evaluate the accuracy of the occultation system, computer inver- 
sion techniques were used to derive pressure references using simu- 
lated infrared and occultation data. Actual atmospheric temperature 
profiles were used to derive infrared measurements and from these 
"derived ,f temperature profiles were obtained by inversion techniques. 
(The derived profiles have no absolute altitude reference.) The 
actual and "derived" temperature profiles were supplied by Dr. 
William Smith of the National Oceanic and Atmospheric Administra- 
tion (NOAA). Next the actual temperature profiles were used to 
calculate atmospheric density profiles, and raytracing programs 
calculated "actual" occultation raypath measurements. To these 
measurements were added 0„5 M rms random errors to simulate 
occultation measurement errors. The noisy occultation data and 
"derived" temperature profiles were used in a computer inversion 
program to obtain estimates of the 5 OO mb levels. This was done 
for a range of longitudes and weather conditions. The estimated 
and actual levels were compared yielding less than 2^ M rms errors, 
which meets the requirements of the weather program. Checks showed 
that the majority of the error comes from the temperature inver- 
sion inaccuracy rather than the occultation system inaccuracy 
and thus the occultation pressure reference results will improve 
as the infrared sounding techniques improve. 


Water Vapor Effects 

•Radio waves are effected by water vapor in the atmosphere as well 
as by the air density itself. At the higher altitudes water vapor 
density is not enough to cause a significant effect; but at lower 
altitudes, especially in the temperate regions in the summer, it 
must be considered. In the above runs the effects of water vapor 
were included and different methods of correction were evaluated. 

To include the effects, actual water vapor profiles for each test 
day were included in the occultation raytracing procedure. The 
occultation results and derived temperature profiles were obtained 
as before. This information was then used together with different 
estimates of water vapor in the inversion program to obtain estimates 
of the pressure level. Ignoring water vapor all together, caused 
systematic pressure level errors in the summer temperate regions 
that exceeded the desired precision by three to four times. Using 
the historical average seasonal water vapor profiles for these 
regions reduced the errors to the 2b M rms requirements. And, as 
might be expected, using the actual dally water vapor profiles 
in the inversion, reduced the errors to levels obtained with dry 
atmospheres . 
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■♦ir'jb u h oven crude climactic data is sufficient to correct for the 

effect of water vapor on the occultation pressure-reference results. 
In an operational weather system, better than climactic water vapor 
information would be available from other components of the observa- 
tion or forecasting system and would be used in the occultation data 
inversion. 

Multipath 

•Multiple ray paths due to water vapor layers in the lower atmos- 
phere is a phenomenon observed during system demonstrations in 
Hawaii. In these demonstrations equipment similar to that propo- 
sed for the occultation measurements was operated between two 
islands. Amplitude and phase path length measurements were made. 

The phase measurements confirmed equipment precision and propaga- 
tion physics. However, during a number of days deep long-term 
amplitude fading was observed. Computer raytracing was carried 
out using atmospheric refractivity profiles obtained by airplane 
flights during the demonstrations. These computations confirmed 
that on clays of deep fades conditions required for persistant 
deep fading were satisfied. These conditions are that there be 
two and only two ray paths between 1 transmitter and receiver, that 

. the two signals be very nearly equal in amplitude, and that the 
paths differ in length by a small number of wavelengths. When 
these conditions are satisfied, two signals can precisely cancel 
each other periodically, causing deep fades. In the Hawaiian 
geometry these conditions can occur fairly often. More than two 
rays, unequal amplitudes, or a rapid variation of path length 
difference maKe* the probability of precise cancellation very, 
very small. 

The ray tracing programs were used with the Hawaiian atmospheric 
refractivity profiles to determine the probability of multipath 
for occultation satellite geometries. The probability of multi- 
path occurrence will be less than that observed in Hawaii for 
two reasons: first the primary occultation raypaths will be much 

higher in the atmosphere, where there is less water vapor to 
distort them; and secondly, the humidity and layering conditions 
in Hawaii are much more severe than is typical of other geographic 
areas and seasons. 

The probability of more than one ray path reaching the occultation 
satellite is thus reduced in comparison with the Hawaiian obser- 
vations. Furthermore, when multiple paths do occur, the probability 
that they will cause deep or prolonged fading is very small; none 
of the three conditions for such fading is satisfied in the satellite 
geometry. For the satellites three or more rays are involved in any 
multipath occurrence not just two; this results from the transmitter 
and receiver being located above the refracting layers, rather than 
in them. Secondly the multiple rays that reach the satellite pass 
long distances through aignif leant ly different parts of the atraos- 
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phere, undergoing different focusing and absorption, which changes 
their relative amplitudes <> Thirdly, the phase relationships 
between the multiple rays will change rapidly as the satellites 
(and ray paths between them) move at orbital speeds over the 
weather patterns* En the Hawaiian experiments, the weather moved 
through the ray paths at a much lower velocity. The probability 
of three or more rays arriving with the right amplitudes and phases 
to closely cancel each other for any length of time will thus be 
very, very small* 

The ray tracing programs also showed that when more than one ray 
existed between transmitter and receiver, the difference in path 
lengths between the rays was typically less than one-half meter. 
This would indicate that if the equipment was designed to measure 
path length in the presence of multiple rays, it would yield 
usefully accurate results even when more than one ray is received. 
It was not clear from the ray tracing why this occurred, and more 
simulations should be done before a system is designed to take 
advantage of the phenomenon 0 

The analysis concluded that, while multipath could be expected 
to occur, it would not occur with high probability, that when it 
did occur, the probability of deep fading would be remote, and 
that results obtained during multiple rays would still meet the 
desired measurement accuracy. 

Orbit Analysis 


•In the occultation technique, the basic measurement is the effective 
length of the radio path between two orbiting satellites. When this 
path passes through the atmosphere, its effective length is increased 
by bending and retardation. To determine atmospheric pressure to 
meteorologically useful accuracy, changes in this increase have to 
be measured to an accuracy of 0«5 M. To accomplish this, of course, 
requires knowledge of the relative movement of the satellites to the 
same accuracy* 

To test the predictability of satellite spacing, first computer models 
were used to simulate the motions of the two satellites. All pertur- 
bations of consequence were included; where actual disturbances 
were unknown, e*g*, the higher order components of the gravity field 
or air drag, random numbers were used to simulate the effect with 
the expected magnitude* 

Next this "actual 11 motion was used to derive the information that 
would be available to the proposed occultation system. In the 
proposal system the satellites would be launched together and then 
caused to drift apart over a period of two weeks until they reach 
the occultation separation. During this drift period, continuous 
path length measurement between the satellites would be available 
undistorted by atmospheric bending. After this time, essentially 
continuous data would still be available but now actual motions 
of the satellites would be mixed with path lengthening caused by 
the atmospheric bending, the data for the meteorologists. 
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Orbit modeling programs then used the drift period data to find 
coefficients for the orbit elements of both satellites, for all 
significant elements of the geopotential field and for other fac- 
tors (such as air drag) effecting the satellite motions. These 
coefficients were used in the prediction of satellite motion 
during occultation. The predicted and "actual" were compared to 
find the kind of errors that would result with this technique. 

The resulting errors were all acceptable for the meteorological 
application, but for different reasons. 

Perturbations due to the high order components of the earth* s 
gravitational field were predictable well within the 0.5 M 
accuracy. The drift-period data (actually accurate to 0.03 M) 
was sufficient to determine the geopotential coefficients to much 
better accuracy than required. Once these fixed coefficients are 
determined, relatively crude knowledge of the satellites* positions 
in the geopotential field during occultation is sufficient to predict 
the geopotential perturbations of both satellites. 

On the other hand, relative drifts of the two satellites due to 
different drag coefficients could easily accumulate to a steadily 
increasing separation error exceeding several meters in days. 

However, there is no way to confuse such a steady increase in path 
length with an actual meteorological phenomenon. (is the atmosphere 
all over the globe expanding at a steady rate?) Data obtained during 
occultation could be processed so as to attribute any such steady 
path length growth to drift in the orbit and the data corrected 
accordingly. 

In a similar way, worst-case unknown air drag variation can cause 
an error in prediction of variation of the orbits* eccentricities, 
which can amount to a meter or so once-per-orbi t oscillation error 
in the satellite separation prediction. Again, there is not an 
atmospheric phenomena that would cause a worldwide smooth atmos- 
pheric bulge of the form indicated by the eccentricity variation. 

The data would be processed to filter out the eccentricity harmonics 
from the path length measurement and to correct eccentricity 
predictions „ 

It is assumed that the satellites would use attitude control systems 
that would fire gas jets only rarely (unloading the momentum in a 
momentum wheel for instance). Furthermore, control of relative 
positioning of the satellites would be accomplished by gas thrusters 
similar to those used for synchronous communications satellites. In 
this event, changes in relative satellite drift rates from these 
causes would occur at known times and could also be easily filtered 
from the data. 

For any of these systematic errors, in actual operation it may be 
advantageous to determine the corrections by comparing occultation 
measurements against known atmospheric conditions rather than by 
statistical filtering. Such calibration data is readily available 
from rawinsonde measurements over the populated areas of the northern 
hemisphere. 
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These simulations thus showed that an occultation system can detect 
the variations caused by the weather patterns with the required 0.5 M 
accuracy,, All errors larger than this are systematic and of a form 
that can readily be separated from actual meteorological information. 
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The work in the following reports thus indicates that an occultation 
system would be able to obtain worldwide pressure reference information 
to the accuracy required by meteorological forecasters. 

The original intent of the study was to proceed from this analysis 
to an interface analysis for including the occultation system in existing 
satellite projects. The conclusions, however , are that the configuration 
proposed originally in 1968 (Report no. RL 5-68) were an adequate solution. 
It would serve no purpose to conduct further engineering design tradeoffs 
unless a specific mission was identified. While the occultation equipment 
is simple, the details of including the equipment on a mission is very much 
dependent on the design of the spacecraft and launch vehicle. Antenna 
sizes and transmitter powers on both master and slave satellites depend on 
the spacecraft attitude control technique and precision, orbit adjustment 
needs, and on space, power, and weight capabilities. 

Judging from past design exercises and from knowledge of present 
transponder costs, launch vehicle cost, and cost of minisatellites, the 
cost to add occultation to an operational weather satellite may be expected 
to be in the $4 to ^6 million dollar range. But a study of an actual 
implementation is needled for an accurate estimate. 

For the transponder design two guidelines have resulted from the 
research. The originally proposed equipment had two techniques for 
measuring the path length, a coherent phase path counter accurate to 
0.03 M change in path length, and a modulation technique determining 
absolute group path length to 0,5 M accuracy, (Both techniques are used 
in current spacecraft tracking equipment.) It is recommended that ambigui- 
ties of 100 meters or more be allowed in the group path measurement if 
this simplifies equipment design, Ground information will be able to 
resolve such ambiguities c 
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In addition, the group path measurement should account for the 
expected characteristics of multipath when it does occur; namely, that the 
lengths of the different paths will differ by several wavelengths but by less 
than the desired 0„5 M accuracy. Thus a modulation technique that under 
these conditions will measure the group path length of the strongest 
signal or of any average of the signals will give useful results even 
during multipath. Amplitude modulation techniques will yield this result, 
frequency or phase modulation will not. 

o o o o o 

The above comments summarize the status of information on the occupa- 
tion techniques. The attached reports in extensive detail show that it 
is feasible for these techniques to obtain meteorlogical ly useful infor- 
mation. This technique should be compared in cost and accuracy with other 
means available to obtain similar data. 
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ABSTRACT 


The microwave occultation technique has proven useful in sensing 
the neutral atmosphere of terrestrial planets such as Mars and Venus. 

An adaptation of this technique can provide meteorological information 
for the Earth. 

A two-satellite microwave occultation system is described that 
will fix as an absolute function of altitude the pressure-temperature 
profile generated by a passive infrared sounder. This method would be 
successful in determining the altitude of the 300 mb pressure level to 
within 24 m rms, assuming the temperature errors prod viced by the IR 
sensor are not greater than 2 °K rms . Error caused by water vapor in 
the radio path can be corrected by climatological adjustments (by use 
of mean water-vapor profiles) if accurate water-vapor sensors are not 
available . 

A ground test of the proposed system is described. A microwave 
signal propagating between. two mountain tops was found to be subject 
to periods of intense fading. Computer analysis of the raypath between 
the transmitting and receiving stations indicates that multipath and 
defocusing were responsible for this fading. Because these phenomena 
are associated with lower altitudes than the closest approach altitude 
of an occultat ion-system raypath, it is unlikely that an operational 
pressure-reference-level system will be subject to the deep fades ob- 
served in the ground test . 
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Chapter I 


INTRODUCTION 

"Every science and every inquiry, and similarly every activity 
and pursuit j is thought to aim at some good." , 

— Aristotle 

To achieve 10-day weather forecasts with 80 percent accuracy (the 
same accuracy now available for overnight forecasts), new techniques 
must be developed. The partial differential equations that describe 
the behavior of the atmosphere must be written in a form that permits 
their solution by digital computers; faster computers must be developed 
to solve these equations in an acceptably short time. Data-acquisit ion 
systems must be designed and implemented, which will sense important 
atmospheric variables on a global basis and report their findings in 
essentially real time. • - ' 

In 1966, Stanford University published a report [1] that described 
a possible meteorological data-acquisit ion system. _ One of the proposed 
techniques described was the application of microwave occultation to 
the Earth’s atmosphere. Occultation (which is based on the fact that 
the atmosphere will bend and retard a radio wave passing through it) 
had been demonstrated with Mariner IV, the space probe that flew past 
Mars; an occultation system had received a microwave signal from Earth 
as the spacecraft moved behind the planet. When” the line of sight be- 
tween Mariner IV and Earth passed through the Martian atmosphere, the 
radio signal was bent and retarded (Fig. 1).. : By analyzing this bending 
and retardation, researchers at Stanford and at the Jet Propulsion Lab- 
oratory were able to deduce the density of the Martian atmosphere (21; 
model matching enabled them to estimate the atmospheric composition and 
the pressure and temperature profiles. 

SPINMAP [1] proposed that a similar system be utilized to sense 
the density of the Earth’s atmosphere, and thus supply meteorologists 
with temperature and pressure profiles of the atmosphere on a global 
basis. The SPINMAP configuration (Fig, 2) consisted of a mother satel- 
lite and several daughter satellites, all in the same near-polar orbit. 
The satellites would be spaced puch that the radio path between the 
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mariner 



TO EARTH 

Fig . 1 , TKS MARINER IV OCCULTATION EXPERI- 
MENT. 

mother and each of the daughters passed through the Earth's atmosphere. 
The fact that there was more than one daughter allowed the atmosphere 
to be "sampled" by rays that had different minimum altitudes; thus, a 
continuous "occultation" could be achieved. As the satellites proceeded 
in their orbit, they would continuously sense the atmosphere. The re- 
sults could be telemetered back to Earth, the data inverted, and pres- 
sure and temperature profiles generated. 

Because the Earth's atmosphere is more complex than that of Mars 
(or even that of Venus whose atmosphere was sensed by an occultation 
experiment on board Mariner V in 1988), objections were raised as to 
the desirability of using occultation satellites in preference to other 
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simpler techniques, such as passive infrared sounders. Pomalaza [3] 
discussed many of the objections and showed how empirical orthogonal 
functions can be used to invert occultation data to yield pressure 
profiles accurate to within 0.3 mb at the surface. 

In 1969, an infrared sensor was flown successfully aboard a NIMBUS 
satellite. The SIRS instrument (satellite infrared sounder) provided, 
as output, temperature as a function of pressure for values between 0.1 
and 1000 mb (the sensor actually provides only nine parameters, but em- 
pirical orthogonal functions are used to develop a least-squares temper- 
ature-pressure profile). The SIRS instrument, however, did not provide 
an altitude reference that would fix the temperature-pressure profile 
as a function of height. Although the hypsometric formula could be 
used to determine the relative altitude between any two pressure levels, 
a reference value was required that would permit these relative altitudes 
to be expressed as absolute altitudes. 

At a meeting of the Radio Occultation Working Group in Washington 
in September 1969, it was suggested that occultation satellites could 
be used to provide this pressure-reference altitude. This report will 
demonstrate the feasibility of a two-satellite occultation technique to 
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supplement the infrared sensor system by providing an accurate altitude 
reference that will serve to fix, as a function of height, the tempera- 
ture profile of the SIRS instrument. 

The results of ground tests made in Hawaii during June 1970 also 
will be described. These tests were conducted to estimate the likely 
effects of scintillation and fading on an occultation system. It was 
found that the microwave signal suffered periods of intense fading; 
however, it will be shown that the probable cause of the fading was 
multipath, which is a low-altitude phenomenon and, therefore, is not 
apt to affect the proposed pres sure -reference system which operates 
at a relatively high altitude («8km). 
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Chapter I I 


BASIC PRINCIPLES 

"Philosophy is written in this grand book — I mean the universe — 
which stands continually open to our gaze, but it cannot 
be understood unless one first learns to comprehend the language 
and interpret the characters in which it is written . 

It is written in the language of mathematics..." 

— Galileo Galilei 


A , Refraction 

The occultation technique is based on the fact that when a radio 
wave passes through the atmosphere it is bent and retarded . This hap- 
pens because the atmosphere has an index of refraction n that is 
greater than X, and this parameter will vary from place to place. 

The index of refraction, however, is not much greater than 1; a 
typical value at the surface of the Earth is 1.000314. Because the 
significant variations in the index of refraction are very small, it 
has been found convenient to define a quantity known as refractivity 
N: 

N = (n - 1) X 10 6 (2.1) 

g 

In other words, N is the fractional part of n, multiplied by 10 . 

Refractivity is a function of atmospheric density and, therefore, 
a function of atmospheric pressure and temperature (temperature, pres- 
sure, and density are related by the equation of state or ideal gas law, 
P = pRT, where P = pressure, p = density, T = absolute temperature, 
and R = the universal gas constant). In terms of temperature T, at- 
mospheric pressure P, and the partial pressure of water vapor e, re- 
fractivity can be written as [4] 

N = 77.6 | - 5.6 | + 3.75 X 10 5 (2.2) 

T 

where T is in degrees Kelvin, P and e are in millibars, and N is 
dimensionless. It can be seen from this equation that N is directly 
proportional to pressure and inversely proportional to some power of 
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temperature; in add it ion , it is a very strong function of water-vapor 
pressure because of the polar nature of the water molecule. The first 
term on the right-hand side of Eq . (2.2) is often called the "dry term," 
and the remainder are called the "wet terms." At low altitudes, where 
water-vapor pressures can be expected to be high, the wet terms make a 
significant contribution to the value of ]■$. At higher altitudes (above 
7 km), this contribution is very small and the dry term is principally 
responsible for the value of 

A radio wave in the atmosphere is retarded if N is greater than 
zero because the index of refraction n is the ratio between the speed 
of light in a vacuum c and the speed of light in the medium; therefore, 
N represents the fractional difference between the speed of the radio 
wave in the atmosphere and its speed in a vacuum. The higher the value 
of N, the greater the retardation. (It should be noted that a monochro- 
matic wave is being discussed and, therefore, the distinction between 
group velocity and phase velocity is not of concern. In any case, the 
atmosphere for the frequency range of interest (X-band) is essentially 
nondispersive . In this report, "velocity" of a wave refers to "phase 
velocity . ") 

In addition to retarding the wave, the atmosphere also causes it to 
bend because N (being a function of pressure, temperature, and water- 
vapor pressure) varies in magnitude, both horizontally and vertically. 
From Snell's Law, it is known that an electromagnetic "ray" tends to 
bend in the direction of increasing index of refraction. Mathematically, 
this can be written as [5] 


dX __ dp 

X P 


where X is the wavelength and p is the radius of curvature of the 
raypath (henceforth, p will refer to radius of curvature rather than 
density). The derivatives are taken in the direction of travel of the 
ray . 

The wavelength X is related to frequency v by X = v/v, where 
v is the velocity of the wave. Because v = c/n, 
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and ’ : • 

dA = - dn 
n v 

Substituting these expressions into Eq . (2.3) and rearranging results in 


” = ‘ 5 VTp <2 ' 4) 

which relates the radius of curvature of a ray to the index of refrac- 
tion n and the derivative dn/dp. This derivative can be expressed 
as the sum of two terms; one is proportional to the horizontal gradient 
of n, and the other is proportional to the vertical gradient: 


dn . 'q 

- — = Vn * "i — r 

dp I P 1 


idn . 1 dn . 

I — ! cos 0 + - ^ Bin 0 


(2.5) 


where r and 6 are polar coordinates with origin at the center of 
the Earth, and 0 is the angle between ~p (the radius of curvature 
vector) and r* (the radius vector). Note that, for r* pointed out-? 
ward, dn/dr is negative for a normal atmosphere. 

The radius of curvature p can now be written as 


p = 


|dn/dr | cos 0 - (l/r)(dn/d0) sin 0 


( 2 . 6 ) 


Therefore, the radius of curvature of a ray is found to be proportional 
to the index of refraction of the medium and inversely proportional to 
the gradient of the index of refraction. If the gradients are large, 
p will be small and the ray will curve sharply; if they are small and 
if n is large, the ray will follow a "straight" trajectory. 

Both the bending of the ray and its retardation tend to increase 
the apparent path length of the ray. If the atmosphere is not present, 
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the radio signal travels in a straight line at the speed of light in a 
vacuum. Bending and retardation caused by the atmosphere increase the 
time it takes for a ray to travel between the same endpoints and, thus, 
effectively increase the length of the raypath. 

This increase in path length is called the "phase defect. 

B. The Occultation System 

The proposed two-satellite occultation system (Fig. 3) would mea- 
sure the phase defect and use it to infer the altitude of a pressure 
level. 

The infrared sensor determines the pressure-temperature profile of 
the atmosphere over a particular portion of Earth; simultaneously, the 
occultation system propagates a radio wave through this same segment of 
atmosphere and measures the phase defect. This result is telemetered 
back to the ground 

(It should be noted that, in an operational system, the IR sensor 
can be mounted on one of the occultation satellites instead of on a 
third satellite, as illustrated in Fig. 3. This sensor has a resolu- 
tion cell similar in size to that of the occultation system, and the 
results can be correlated on the ground.) 

The pressure-temperature profile is used as input for a computer 
program that raytraces through the index -of-ref ract ion profile generated 
from the pressure-temperature curve. The results of this raytracing are 
compared to the phase defect measured by the satellites. The index— of- 
ref ract ion profile is then adjusted the proper vertical distance until 
the raytracing result obtained from the computer program matches the 
phase defect obtained from the satellites. 

Several of the parameters of the proposed system are significant 
and should be given special attention. Because an attempt is being made 


The origin of the term "phase defect" is obscure. It may be attributed 
to the fact that P as the path length and wavelength increase, the num- 
ber of zero crossings registered by the receiver decreases, resulting 
in a phase "defect' 1 at the receiver. 
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to measure an excess path length on the order of hundreds of meters over 
a path on the order of 8000 km, it is obvious that the system must be of 
high precision * In fact, only small variations (on the order of meters) 
in the phase defect are of interest. It is important, therefore, to be 
able to measure the phase defect as accurately as possible, hopefully to 
within 1 m or less, which means that the relative positions of the satel- 
lites must be known to within this distance. 

The factors affecting satellite position are atmospheric drag, solar- 
radiation pressure, and perturbations in the orbit attributed to anomalies 
in the Earth’s gravitational field. Drag and solar-radiation pressure 
will have periodic effects that can be removed from the data by digital 
filtering. The drag effect is a long-term phenomenon; the solar-radia- 
tion effect will have a period of half a revolution and should be remov- 
able without undue difficulty. 

Orbital perturbations caused by anomalies in the Earth’s gravita- 
tional field are a more difficult problem. Theoretical analysis [6] has 
shown that, by launching the mother and daughter satellites together and 
then letting the daughter drift away, the gravitational field can be 
mapped with sufficient accuracy to remove the effects of the anomalies. 

Another significant parameter is the nominal separation distance 
between the satellites. By increasing this distance, the raypath be- 
tween the satellites can be made to approach closer to the Earth’s sur- 
face at its lowest altitude. In addition, the Earth is somewhat oblate; 
the difference in radius measured at the equator and at the pole is ^20 
km. As will be shown in Chapter III, satellite separation and the ef- 
fects of the Earth's oblateness are important factors in system design. 

Figure 4 is a block diagram of the system. The mother satellite 
transmits a CW signal through the Earth’s atmosphere. The daughter sat- 
ellite receives the signal, multiplies its frequency by some rational 
fraction (such as 11/10), and retransmits it phase coherently to the 
mother satellite where it is multiplied by the inverse of the original 
multiplying factor (in this case, by 10/11). The signal now beats with 
the original frequency, and the beats are counted and recorded. If, in 
passing through the atmosphere-, the frequency of the signal was not af- 
fected, the result of beating will be zero; however, if the atmosphere 
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Fig. 4. BLOCK DIAGRAM OF THE OCCULTAT ION-SYSTEM HARDWARE. 

causes the apparent distance between the satellites to change, the 
frequency of the signal will have shifted slightly as a. result of the 
Doppler effect, and this shift will appear as a beat frequency. Half, 
the number of beats/sec produces the number of wavelengths of excess 
path length caused by the atmosphere (or, as noted above, caused by real 
path-length changes resulting from variations in satellite position). 

In addition to the system described above, which counts every wave- 
length change in distance between the satellites, the signal of the mas- 
ter satellite is phase modulated at «5MHz. This modulation is carried 
through the repeater satellite and received back at the master. Phase 
comparison of the received and transmitted signals yields an absolute- 
distance measurement of approximately 1/2 m accuracy with a basic ambi- 
guity of «60m; this ambiguity can be increased to approximately 3 km 
by adding a second modulation frequency if system considerations so 
require. The absolute-distance system provides the original absolute- 
distance reference and reestablishes the absolute reference whenever 
the more precise system loses count because of interference. The pos- 
sibility of signal loss from multipath is discussed extensively in later 
chapters . 

It should be noted that hardware similar to that described above 
has been flight tested and proven on numerous planetary and interplane- 
tary missions, including Mariner, Pioneer, Lunar Orbiter, and Apollo. 
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PRECEDING PAOT BLANK NOT FILMED 

Chapter III 

ERROR ANALYSIS OF THE PRESSURE-REFERENCE-LEVEL SYSTEM 

"Errare humamim est." 

— Ancient Saying 

A computer simulation of the occultation system was necessary to 
determine the feasibility of using occultation satellites to measure 
the altitude of a pressure level. This chapter describes the program 
used and the results obtained . 

A . The Pressure-Reference-Level Program 

The computer program was divided into two parts. In the first, 
a model atmosphere, consisting of the best temperature, pressure, and 
water-vapor information available, was input to the computer to con- 
struct an index -of -ref raction profile. A raytracing routine traced 
through this profile, simulating the result that would be obtained if 
a satellite system were to transmit through this model atmosphere. The 
result was a phase defect (the excess path length) and was called SR, 
meaning "satellite reading." 

In the second, the information used to construct the index-of- 
ref raction profile was degraded so as to simulate the errors inherent 
in the satellite infrared sensor (SIRS) which, in an operational system, 
would be the source for this data. (The particular species of degrada- 
tion employed varied from one run to another and will be described later 
in this chapter.) All information relating pressure or temperature to 
altitude was discarded, and the bottom-most temperature-pressure pair 
was arbitrarily fixed at some starting height. The hypsometric formula 
was used to assign corresponding altitudes to the remaining pressure- 
temperature pairs, an index-of -ref raction profile was constructed, and 
the ray tracing routine was employed to generate a value of the phase 
defect. (This value was denoted as PD1 .) The profile was then shifted 
vertically 1 km, the ray tracing routine was again applied, and a new 
value (PD2) was obtained. These two values were used to calculate the 
derivative of the phase defect with respect to altitude d0/dz, where 
d0/dz - PD2 -PDX. 
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Now, the SR of the first part of the experiment (the phase defect 
‘'measured" by the imaginary satellite system) was employed to construct 
an error term, 


A0 = SR - PD2 (3.1) 

This value A0 was used in conjunction with the derivative d0/dz to 
determine the shift Az necessary to bring the error term A0 to zero: 


Az 


A0 

d0/dz 


(3.2) 


In practice, the system is somewhat nonlinear, and more than one itera- 
tion is generally necessary before A0 is brought arbitrarily close to 
zero . 

When A0 had been made acceptably small through repeated raytrae- 
ing, application of Eq. (3,2), and the appropriate vertical shift of the 
index-of-ref raction profile, the altitudes of various pressure levels 
were compared to the altitudes of these same levels in the "measured" 
atmosphere profile used in the first part of the program. The differ- 
ences in altitude constituted the error of the system; therefore, if the 
500 mb level were at 5.005 km in the reference atmosphere and at 5.015 
km in the atmosphere generated by the shifting procedure, the error at 
500 mb would be 5.015 - 5.005 km = 0.010 km, or 10 m. Figure 5 is a 
flow chart of this program, and its listing is found in Appendix A. 

In the remainder of this chapter, the results of the final error 
analysis are discussed, in which simulated data from the satellite in- 
frared sounder (SIRS) were used to generate the reference atmosphere. 


B. Error Analysis, Using Simulated SIRS Data 

The simulated SIRS profiles [7] consisted of 206 profiles of tem- 
perature vs pressure, with temperatures at 100 pressure levels between 
1000,0 and 0.1 mb. A typical data set is shown in Fig. 6. Each level 
presents a simulated temperature "sensed" by the SIRS, a temperature 
"error" (how this sensed temperature differs from the true temperature), 
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Fig. 5. FLOW CHART OF THE PRESSURE -REFERENCE -LEVEL PROGRAM. 
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and the errosr in height resulting from the accumulated error in tempera- 
ture . In addition, each set tabulates rms errors for various ranges of 
pressure (such as TOO to 400 mb, 400 to 100 mb). 

Half of the data sets consist of simulated profiles for a 7-channel 
sensor; the other half consist of profiles for a 9-channel sensor, which 
are the only ones used in this error analysis. The data sets were divided 
into three groups (set 2, set 3, and set 4), according to the latitude of 
the stations, and these latitude ranges are indicated in Table 1. 

The pressure-reference "level program was run by combining the "sensed" 
temperature profile and the corresponding temperature errors so as to pro- 
duce the "true" temperature profile which became the "reference atmosphere." 
The original sensed -temperature profile was used as "measured data." An 
atmosphere was constructed from this profile and shifted vertically until 
the results of raytracing through the vertically shifted measured profile 
matched the results of ray tracing through the reference atmosphere. In 
both portions of the program, a random error of absolute magnitude (0.5m) 
was added to the phase defect to simulate the uncertainty in satellite 
spac ing. 

The resulting altitude error at each of nine pressure levels was re- 
corded. Figure 7 is a typical output sample. Statistical summaries for 
each of the three latitude ranges are tabulated in Tables 2 through 7. 

Tables 2, 3, and 4 list the results for a satellite separation of 7863 km f 
and Tables 5, 6, and 7 list the results for 7888 km. (The greater the sat- 
ellite separation, the lower is the closest approach altitude of the ray.) 
Table lb lists the closest approach altitudes for each of the cases shown 
in Tables 2 to 7; the results are also displayed in Figs. 8 and 9. 

A smaller number of profiles were run with water vapor introduced 
into the reference-atmosphere portion of the program. This water-vapor 
data consisted of the associated rawinsonde relative humidity profiles. 

A simplistic water-vapor correction was used in the second half (measured- 
data portion) of the program* This correction consisted of one of two 
humidity profiles, one representing a "high” value of water vapor and 
the other a "low" value (see Fig. 10), chosen according to the surface 
value of the rawinsonde water-vapor profile* All of the data sets from 
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Table 1 


SIGNIFICANT DATA FOE PRESSURE -REFERENCE -SYSTEM ERROR ANALYSIS, 
EMPLOYING SIMULATED SIRS TEMPERATURE PROFILES 


Latitude 

Set No, 

2 

3 

.. 

4 

Min 

Max 

25 N 
37 N 

1 

51 N 
76 N 


a. Latitude ranges 


Set 


Satellite 

Separation 



7863 km 

7888 km 

No. Prof 

Min Alt 

Max Alt 

No. Prof 

Min Alt 

Max Alt 

2 

14 

5 .78 

7.86 

14 

3.59 

5 .44 

3 

20 

7,55 

9.75 

19 

4.85 

6.81 

4 

18 

9.59 

12.56 

18 

7.29 

9.91 


b„ Dry Atmosphere 


Set 


Satellite 

Separat ion 



7863 km 

7888 km 

No. Prof 

Min Alt 

Max Alt 

No. Prof 

Min Alt 

Max Alt 

2 

6 

6 .25 

7.59 

6 

4.11 

4.95 

3 

8 

7.62 

9.96 

8 

4.99 

7.10 

4 

7 

9.70 

11.47 

7 

7.40 

9.28 


c. Wet atmosphere with simplistic climatic correction 



Satellite Separation 

Set 

7863 


No. Prof 

Min Alt 

Max Alt 

2 

6 

6.11 

7.65 

3 

8 

7.69 

10.00 

4 

7 

9 0 98 

11,85 


do Wet atmosphere with idea.l correction 
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Table 2 


SATELLITE SEPARATION: 7863.0000 KM (14 Stations); SET 2 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max 

500 1.000 

7.70 

14.05 

16.02 

-16 . 89 

28.33 

400.000 

8.18 

16.42 

18.35 

-21.21 

30.48 

300.000 

17.64 

28.40 

33.43 

-30.50 

53.93 

200.000 

24.65 

35.57 

43.28 

-25.81 

92.03 

140.000 

2.64 

33.55 

33.65 

-55 .26 

61.61 

100.000 

0.57 

27.10 

27.10 

-45 .36 

58.87 

70.000 

4.56 

28.17 

28.54 

-42 .50 

56.64 

50.000 

11.02 

38.20 

39.76 

-58.88 

91.92 

30.000 

10.55 

37.38 

38.84 

-49.27 

58.26 


Table 3 


SATELLITE SEPARATION: 7863.0000 KM (20 Stations); SET 3 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max • 

500 ^000 

1.21 

17.67 

17.71 

-20.61 

36.88 

400.000 

-3.15 

12.97 

13.35 

-24.05 

17.89 

300.000 

-0.07 

16.47 

16.47 

-43 . 97 

21.99 

200.000 

6.28 

25.77 

26.52 

-51,67 

41.15 

140.000 

4.90 

42.22 

42.50 

-83.59 

80.43 

100.000 

-2.32 

39.53 

39.59 

-73.74 

95.77 

70.000 

1.13 

28.35 

28.38 

-33.48 

87.17 

50.000 

9.72 

24.72 

26.57 

-28.99 

77.82 

30.000 

24.42 

40.21 

47.05 

-47.99 

87.49 


Table 4 


SATELLITE SEPARATION: 7863.0000 KM (18 Stations); SET 4 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max 

500.000 

12.38 

22.83 

25.97 

-27.09 

48.48 

400.000 

9.56 

25.29 

27.04 

-31.75 

42.22 

300.000 

-0.82 

25 .32 

25.33 

-49.44 

38.81 

200.000 

12.11 

22.26 

25,34 

-12.92 

72.49 

140.000 

14.80 

29.62 

33.11 

-21.11 

95.06 

100.000 

15.39 

33.36 

36.74 

-38.29 

105.33 

70.000 

11c 87 

29.70 

31.99 

-41.17 

76.05 

50.000 

1 o35 

19.38 

19.42 

-31.85 

30.24 

30 o 000 

-21,75 

36.20 

42.23 

-78,40 

43.18 
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Fig. 8. STANDARD DEVIATION AND RMS ERROR VS PRESSURE. Dry atmosphere, satellite separation: 7863 .Okm, 
Closest approach altitude of ray indicated by vertical line. 


Table 5 


SATELLITE SEPARATION: 7888,0000 KM (14 Stations); SET 2 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max 

500.000 

2.72 

10,82 

17.04 

-32.38 

30.76 

400.000 

3.20 

17.42 

17.71 

-28.07 

28.50 

300.000 

12.66 

15,84 

20.27 

-14.57 

34.83 

200.000 

19.67 

22,48 

29.87 

-13,67 

59.76 

140.000 

-2.34 

32.35 

32,44 

-47 , 95 

48.25 

100.000 

-4.41 

28,62 

28,96 

-62.93 

42.37 

70.000 

-0.42 

23.34 

23.34 

-47.52 

32.24 

50,000 

6.04 

29.24 

29.86 

-32.20 

68.65 

30.000 

5.57 

25.40 

26.00 

-44.91 

44.65 


Table 6 


SATELLITE SEPARATION: 7888.0000 KM (19 Stations); SET 3 


Level (mb) 

Mean 

St . Dev . 

rms 

Min 

Max 

500,000 

-3.31 

19,62 

19.90 

-43.18 

27.33 

400.000 

-7.95 

28.79 

29.87 

-68.50 

41.94 

300,000 

-5.63 

45 .71 

46.06 

-96.07 

69.94 

200,000 

-0.87 

43.21 

43.22 

-80.68 

71.70 

140.000 

-0,31 

38,60 

38,60 

-68.58 

68.34 

100,000 

-6.49 

37.56 

38.12 

-60.89 

73.51 

70.000 

-2.94 

30.74 

30.88 

-51.80 

72.04 

50.000 

4.29 

25.10 

25 .47 

-33.07 

51.54 

30.000 

15.71 

37.80 

40.93 

-47 .12 

89.89 


Table 7 


SATELLITE SEPARATION: 7888.0000 KM (18 Stations); SET 4 


Level (mb) 

Mean 

St , Dev . 

rms 

Min 

Max 

500.000 

14.07 

35.69 

38.36 

-29.69 

90.63 

400 .000 

11*25 

28,34 

30.49 

-17.20 

76.93 

300.000 

0.87 

27,49 

27.51 

-40.40 

56.26 

200.000 

13.81 

49.56 

51.45 

-56.19 

103.93 

140.000 

16.49 

58,22 

60.51 

-60.69 

133.87 

100.000 

17.08 

54.47 

57.08 

-40.43 

140.75 

70.000 

13.56 

50.62 

52.40 

-47,56 

127.72 

50.000 

3.04 

49.95 

50.04 

-79.01 

109.18 

30.000 

-20.06 

63.61 

66.70 

-134.16 

83.86 
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Fig. 9. STANDARD DEVIATION AND RMS ERROR VS PRESSURE. Dry atmosphere, satellite separation; 7888.0km, 
Closest approach altitude of ray indicated by vertical line. ' 





Fig.. 10. WATER VAPOR PROFILES, The above profiles were 
extrapolated to 25 km by lise of the equation e(z) = 
e(7) ■* exp [0 c54?S # (z-7)], where e is the water- 
vapor pressure, in millibars, and z is the altitude, 
in kilometers. 
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sets 3 and 4, and one from set 2 used the lower corrective profile; the 
remaining data sets in set 2 used the higher corrective profile. Tables 
8, 9, and 10 summarize the results for a satellite separation of 7863 km 
and Tables 11, 12, and 13 summarize the results for 7888 km. The corre- 
sponding closest approach altitudes are listed in Table lc, and the re- 
sults are illustrated in Figs. 11 and 12. 

To determine the effect of an improvement in the water-vapor correc 
tion technique (a more sophisticated approach than the simplistic method 
described above), the same data were run with the exact rawinsonde water 
vapor profile used in both the reference and measured atmospheres, which 
is equivalent to an ideal water-vapor correction. Results of this run 
(for the 7863 km satellite separation only) are tabulated in Tables 14 
through 16. The corresponding closest approach altitudes are listed in 
Table Id, and the results are displayed in Fig. 13. 

Finally, to determine the effects of the temperature errors on the 
system, the analysis of six low-latitude profiles (without water) was 
repeated,, this time using the reference temperature values (the values 
without error) in both the reference and measured atmospheres, which is 
equivalent to having an ideal -temperature sensor. The result of this 
run is illustrated in Fig. 14 for the 7863 km satellite separation only. 


C . Analysis of Results, Using Simulated SIRS Data 

In Figs. 8, 9, 11, and 12, the solid line represents the standard 

deviation of the error (in meters), and the dashed line represents the 

root mean square (rms) error; both are shown as a function of pressure. 

The difference between the standard deviation and the rms values at each 

pressure level is an indication of the size of the mean error at this 

level. The rms value is equal to the square root of the sum of the 

2 

squares of the standard deviation and the mean (rms = standard devia- 

2 2 
tion + mean ). 

A cursory examination of these figures reveals that the curve of 
altitude error vs pressure follows a typical pattern, regardless of 
whether the atmosphere is wet or dry or, if wet, how the water-vapor 
error is corrected. Typically, there are two minima and two maxima. 
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Table 8 


SATELLITE SEPARATION: 7863.0000 KM (6 Stations); SET 2 




Level (mb) 

Mean 

St. Dev. 

rms 


Max 

500.000 

13.12 

16.80 

21,31 

-7.95 

40.43 

400.000 

8.31 

20.21 

21.85 

-23 .56 

29.32 

300.000 

15.22 

26.48 

30.54 

-27.64 

40.54 

200.000 

34,51 

34.32 

48.67 

-2 . 97 

93.47 

140.000 

32.63 

47 . 46 , 

57.63 

-18.32 

80.53 

100 .000 

19.16 

31.17 

36.59 . 

-26.20 

59.18 

70.000 

6.13 

43.39 

43.82 

-56.23 

46.92 

50.000 

11.06 

59.10 

60.13 

-52.70 

83.33 

30.000 

15 .58 

26.30 

30.57 

-19.93 

46.52 


Table 9 


SATELLITE SEPARATION: 7863.0000 KM (8 Stations); SET 3 


Level (mb) 

Mean 

St o Dev . 

rms 

Min 

Max 

500.000 

39,18 

33.84 

51.77 

-5.46 

99.36 

400.000 

30.87 

30.54 

43.42 

-1 .78 

88,95 

300.000 

21.53 

29.34 

36.39 

-16.31 

68.32 

200,000 

38.97 

38,62 

54.86 

-3.64 

94.15 

140.000 

50.94 

48.82 

70,56 

-22.24 

121.11 

100.000 

44,02 

46.73 

64.20 

-35.13 

105 . 45 

70.000 

38.30 

34.71 

51 .69 

-13.50 

88.64 

50.000 

32.68 

31.24 

45 .21 

-21,36 

70.95 

30.000 

43.89 

40.11 

59.45 

-33.54 

85.94 


Table 10 


SATELLITE SEPARATION: 7863.0000 KM (7 Stations); SET 4 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max 

500.000 

25 .19 

27.94 


-20.77 

66.63 

400.000 

23.62 

25,98 


-26 . 80 

52.06 

300.000 

9.41 

8.46 


-3.93 

18.65 

200.000 

13.59 

14.00 

19.51 

-7.75 

35.83 

140.000 

17.46 

21.65 

27.81 

-13.57 

41.90 

100.000 

25.93 

25.36 

36.27 

-23.61 

56.55 

70.000 

29 c 48 

26.43 

39.59 

-20,16 

69.11 

50.000 

14,41 

24,78 

28.66 

-19.81 

43.91 

30.000 

-30.30 

47,90 

56.68 

-91,29 

39.23 
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Fig. 11. STANDARD DEVIATION AND RMS ERROR VS PRESSURE. Wet atmosphere with simplistic climatic correc- 
tion, satellite separation: 7863.0 km. Closest approach altitude of ray indicated by vertical line. 




Table 11 


SATELLITE SEPARATION: 7888.0000 KM (6 Stations); SET 2 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max 

500.000 

-53.11 

37.21 

64.85 

-97.58 

8.24 

400.000 

-57.93 

41.29 

71.14 

-122.42 

-2.86 

300.000 

-51.02 

45 .27 

68.21 

-126.50 

0.03 

200.000 

-31.72 

47.98 

57.52 

-68.77 

61.28 

140.000 

-33.55 

49.06 

59.44 

-89.27 

48.35 

100.000 

-47.07 

37.68 

60.29 

-77.58 

27.00 

70.000 

-SO. 11 

60 o 83 

85.52 

-142.15 

9.67 

50.000 

-55.17 

73.51 

91.91 

-151.57 

12.85 

30.000 

-50.65 

45.45 

68.73 

-118.79 

5 .25 


Table 12 


SATELLITE SEPARATION: 7888.0000 KM (8 Stations); SET 3 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max 

500.000 

51.29 

57.44 

77.00 

-12.32 

146.25 

400.000 

42.97 

56.47 

70.96 

-31.52 

135.84 

300.000 

33.63 

58.84 

67.77 

-51.88 

115.21 

200.000 

51.08 

63.97 

81.86 

-32 . 43 

135 . 99 

140.000 

63.04 

69.96 

94.18 

-30.39 

168.00 

100.000 

56.13 

72.17 

91.42 

-43.29 

150.06 

70 . 000 

50.40 

63.12 


-21.67 

126.30 

50.000 

50.12 

50.68 


-2.69 

117.84 

30.000 

55.99 

50.55 


-1.54 

125 .75 


Table 13 


SATELLITE SEPARATION: 7888.0000 KM (7 Stations); SET 4 


Level (mb) 

Mean 

St . Dev . 

rms 

Min 

Max 

500.000 

15,78 

48.91 

51.40 

-34.20 

99.85 

400.000 

14.21 

42.20 

44.53 

-37.26 

85.28 

300.000 

0.01 

46.03 

46.03 

-67.46 

67.63 

200.000 

4.18 

64.08 

64.22 

-89.11 

92.39 

140.000 

8.06 

60.09 

60.63 

-94.93 

75.12 

100.000 

16.53 

49.55 

52.23 

-68.07 

67.30 

70.000 

20.07 

44.89 

49.18 

-41.15 

72.56 

50.000 

5.00 

59.20 

59.41 

-65.16 

75.32 

30.000 

-39.71 

93.34 

101 0 44 

-145 . 74 

94.39 
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Table 14 

SATELLITE SEPARATION: 7863.0000 KM (6 Stations); SET 2 


■ 

Level (mb) 

Mean 

St. Dev. 


rms 

Min 

Max 

500.000 

10.37 

8.66 

13.51 

-0.26 

21.08 

400.000 

5 .56 

15 .49 

16.45 

-19.55 

28.87 

300.000 

12 .47 

22.52 

25.74 

“23.63 

45.27 

200.000 

31.76 

27.86 

42.25 

1.24 

71.97 

140.000 

29.93 

45 .5 9 

54.54 

-23.13 

80.99 

100.000 

16.42 

27.90 

32.38 

-18.49 

46.27 

70.000 

3.38 

38.68 

38.83 

-54.34 

56.15 

50.000 

8.31 

55.82 

56.43 

“49.15 

92.56 

30.000 

12.83 

24.88 

28.00 

-15.91 

55.76 


Table 15 


SATELLITE SEPARATION: 7863.0000 KM (8 Stations); SET 3 


Level (mb) 

Mean 

St. Dev. 

rms 

Min 

Max 

500.000 

13.91 

20.29 

24.60 

-18.24 

38.74 

400.000 

5 .59 

17.51 

18.38 

-19.58 

27.90 

300.000 

-3 .74 

18.36 

18.74 

-39.93 

24,18 

200.000 

13.70 

28.19 

31.34 

-20.49 

66.59 

140.000 

25.67 

39.40 

47.02 

-40.89 

92.33 

100.000 

18.75 

37.49 

41 . 92 

-53.79 

77.89 

70.000 

13.02 

26.43 

29.47 

-32.17 

61.07 

50.000 

12.74 

12.80 

18.06 

2.70 

42.34 

30.000 

18.61 

32.47 

37.43 

-46.31 

67.28 


Table 16 


SATELLITE SEPARATION: 7863.0000 KM (7 Stations); SET 4 


Level (mb) 

Mean 

— 
St „ Dev . 

rms 

Min 

Max 

500.000 

20.38 

26.00 

33.03 

-25 .66 

48.63 

400,000. 

18.81 

25.17 

31.42 

-31.70 

45.03 

300.000 

4.60 

11,40 

12.29 

-9.59 

19.66 

200.000 

8.78 

11 o 17 

14.20 

-8.24 

16.14 

140.000 

12.65 

20.65 

24,22 

-13.66 

34.70 

100.000 

21.12 

27.21 

34.45 

-28.50 

58.59 

70.000 

24.67 

29,83 

38.71 

-25.05 

71.14 

50.000 

9.60 

25 .67 

27.41 

-20.39 

45.94 

30.000 

-35.12 

44,69 

56.83 

-100.97 

18.59 
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Fig. 13. STANDARD DEVIATION AND RMS ERROR VS PRESSURE. Wet atmosphere with ideal water-vapor correc- 
tiOHf satellite separation: 7863.0 km. Closest approach altitude of ray indicated by vertical line 





One minima generally occurs in the neighborhood- Of -'300-' to 400* mb", 5 and 
the other in the neighborhood of 50 to 70 mb. -The- maxima are-hear^the 
30 mb level and between 100 and 200 mb (which,- for most prof i:iesV^*is : 
the pressure level of the tropopause). . : -cz iC. jrr< 

The pattern can be explained by considering that t he - minimum* at' 
the lower altitude, near the 300 mb level, is associated- with- - the: clos- 
est approach altitude of the raypath (the vertical line;’ in 'thett-igures ) . 
Because the system is most sensitive to atmospheric effects ^-'approxi- 
mately the lower 3 km of the raypath, the absolute value of error in 
this neighborhood tends to be minimized. '%"•••-■ course- 

The overall ' pattern is the result of -the errors -in^ -thei temperature 
prof lie, ' ; which generate errors- in alt itude when the tbmbera-tu^e-a'l-li%ude 
profile is^ Con's true tedlusing -the hypsometric- formul ■Thl s^ef ife'e t^ i's^- - 

illustrated.: by- the absence of any ext refta^dn;the=hc stilt stJgrf ’thesiid'e^Lor, 
temperature-profile case (Fig . '14);- where ‘ each- • -1 evel-h'as^theJ same e-rrorg ° 
as that- of the 300 to 400 mb level.. «- .• as ■cc.'cs.Viy eliminated by prop- 
erly T-he'^ magnitude of r-t he- error-iisr of primary^ interest Ion ij t&chri^'be "s'ebhl- 
in Figs. 8 and 9 that, generally, the rms and standard^ deviate. oh bellow 2 - 
300 mb are approximately 20 to 30 m. Figure 15 shows the 300 s mfe'^rrbr- 
for each of the profiles used in the generation of Fig-.-6> plotted 0 as® a 
function of rms temperature error in the 100 to 400 mb pressure range. 
Although there is no clear regression line, it should be noted- that- the 
upper left-hand portion of the graph is free of data points; which- indi- 
cates that a low rms value of temperature error in- the 100- to - 400" mb- 2 
range invariably results in a small altitude error at’ 'trie- 360'' r mb--rev%flV- 
It is reasonable, therefore, to expect that improvement^' in" the? perf dr-/ 
mance of the SIRS will lead to concommitant improvements -in -the^' perfor-' 
mance of the occultation system. It should be noted also thatP^a high" 
rms value of temperature error does not lead necessarily" to' Targe er- 
ror in the altitude of the 300 mb level, as witnessed by the large num- 
ber of points in the lower right-hand portion of the figure. 

Comparison of Figs. 8 and 9 reveals that it is possible to choose 
a satellite separation for each of the latitude ranges such that the 
rms and standard deviation at the 300 mb level are less than 26 m. 
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TEMPERATURE RMS ERROR (°K) 


Fig. 15. ALTITUDE ERROR AT 300 mb VS TEMPERATURE RMS ERROR 
Pressure range: 100 to 400 mb. 
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This does not suggest necessarily that there must be more than one 
satellite pair in the system, but it does illustrate that minimum ray 
altitude is a significant parameter that should be selected in such a 
way as to optimize system performance. In fact, the 7863.0 km separa- 
tion obtains acceptable results at all three latitude ranges. Should 
further analysis indicate that variable satellite separation is desir- 
able, a small oblateness can be introduced intentionally into the sat- 
ellite orbit. Note that, in each figure, the lower minimum occurs in 
the range of the closest approach altitude. 

It can be seen in Fig. 11 that. the simplistic water-vapor correc- 
tion produces remarkably good results. The standard deviation and rms 
at the 300 mb level and below are once again in the neighborhood of 20 
to 30 m. Figure 13 illustrates that the errors are caused by the fail- 
ure of the corrective profile to account adequately for the water vapor, 
and these results should be compared to the corresponding curves in Fig. 
8. Errors introduced by water vapor can be totally eliminated by prop- 
erly correcting for the water vapor. If exact corrections are not avail- 
able, even the most simplistic corrections (Fig. 11), in which a standard 
water-vapor profile corrects for water-vapor effects, can have dramatic 
results. Not correcting for water vapor will yield errors on the order 
of hundreds of meters . 

The effect of varying the minimum ray height in the presence of 
water vapor is illustrated in Fig. 12. The system operating at the 
lower ray height (greater satellite separation) yields generally poorer 
results than the system operating at the higher ray height (smaller sat- 
ellite separation). Because water-vapor concentration falls off rapidly 
with altitude (scale height on the order of 2 km), this is not a partic- 
ularly surprising result, as proper water-vapor correction is more im- 
portant as the ray passes into regions of high water-vapor content. 
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Chapter IV 


THE HAWAII EXPERIMENT 

"One would grow poor staying in one place always." 

— Poem of the Cid 


A major concern in the design of an occultation satellite system 
is the possible effect of atmospheric turbulence on system performance, 
consisting chiefly of scintillation of the microwave signal and the 
possible consequent loss of phase lock. In addition, multipath due to 
atmospheric water-vapor inhomogeneities could cause deep fading for 
several seconds or longer, again resulting in loss of phase lock. Al- 
though ambiguity resolution built into the system hardware would make 
it easy to reestablish the phase-path measurement after lock is regained, 
there will be a total loss of data while the signal is in deep fade. 

To assess the problems associated with scintillation and deep fade, 
a ground test was performed between two mountains in Hawaii during June 
1970. The choice of site was determined largely by such logistical fac- 
tors as the availability of power and shelter and ease of access. The 
Hawaii location, also provided a 150 km transmission path between stations 
at a relatively high altitude for a ground test (3.05 and 3.35 km above 
sea level). 

By proper choice of antenna beamwidths, the illuminated area between 
the transmitter and receiver was limited to altitudes between 2 and 4 km 
above sea level, thus avoiding interference from surface reflection. An 
operational pressure-reference system will have a closest approach ray 
altitude in the neighborhood of 8 km, which is much higher (and dryer) 
than the 2 to 4 km altitude sampled in this experiment. It was neces- 
sary, therefore, to extrapolate the results obtained in Hawaii to the 
higher altitude so as to assess the impact of the observed phenomenon 
on the proposed system. Ideally, 8 km mountains were preferred. None 
were available. 

The data acquisition and scintillation analysis were performed by 
a group from the U.S. Department of Commerce Office of Telecommunications 
[ 8 ]. 
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A . Description of Measurements 


Radio -propagation measurements were made in 1970 over a continuous 
two-week period from 0800 June 15 through 0500 June 29. The transmitter 
was located at the University of Hawaii Mees Solar Observatory, at the 
summit of Haleakala on the island of Maui. The receiver was located 150 
km to the southeast, at the Mauna Loa Observatory (at the 11,000 ft level) 
on the island of Hawaii. 

At the transmitter end of the path, three 1.25 m parabolic dishes ' 
(beamwidth 1.8°) transmitted signals at wavelengths very close to 3.4 cm 
(a slight frequency offset at each of the antennas enabled the receiver 
to distinguish between signals transmitted from each antenna). The an- 
tennas were spaced horizontally along a line perpendicular to the nominal 
transmission path, with 10 m between the first and second antennas and 
100 m between the first and the third. The receiver consisted of a sin- 
gle 2 m parabolic dish (1.2° beamwidth) and associated electronics, and 
the received signals were retransmitted to the Haleakala station via a 
frequency -modulated telemetry link. The beamwidths of the receiving and 
transmitting antennas were narrow enough to ensure that the Earth's sur- 
face would not be illuminated within the half -power points of the two 
antennas . 

Signal amplitudes from the three antennas, phase variability of the 
three transmission paths, and phase difference between the first and sec- 
ond and the first and third antennas were recorded continuously. Only 
the signal-amplitude measurements are discussed in this study. 

In addition to the radio measurements, an instrumented aircraft 
equipped with a microwave ref rac tome tez 1 was flown along the raypath to 
measure the refractivity of the atmosphere. Flights were made approxi- 
mately twice a day, on a schedule that ensured sampling from each period 
of the diurnal cycle. 

Figures 16 and 17 are horizontal and vertical views of the transmis- 
t 

sion path. 


t 


All figures in this chapter 


were taken from Ref. 8, 
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Fig, 16. HAWAII PROPAGATION PATH. The transmitter was located on 
Haleakala. Prevailing winds were from the northeast (upper right- 
hand corner). 


B. Significant Observations 

Perhaps the most interesting feature of the observed radio data was 
the presence of periods of deep fading, some as deep as -40 dB. Although 
these very deep fades would typically have a duration of only a few sec- 
onds, on occasion, the signal would drop below -30 dB for as long as 20 
sec. Bays with deep fades were also characterized by many incidents of 
less intense fading. Signal enhancements of up to +5 dB were occasion- 
ally observed. Typical amplitude records from a "quiet" day (one in 
which there was very little intense fading) and from a ’’noisy” day are 
shown in Fig. 18. 
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Fig. 17. VERTICAL SECTION OF THE HAWAII PROPAGATION PATH. The transmitter is on the left. 
The dashed lines indicate the beamwidths of the antennas. Note that only the region from 
2 to 4 km above sea level is illuminated by both the transmitting and receiving antennas. 






The amplitude record was divided into 30 min periods, each classified 
into one of four "fading types" designated A, B, C, and D, corresponding 
to 


Type 

A: 

<5dB peak-to-peak 

Type 

B: 

>5dB peak-to-peak 

Type 

C: 

> 5 dB fading range 

Type 

D: 

>5 dB fading range 


fading range 
fading range 

with 1 to 5 V-shaped fades >15dB 
with over 5 V-shaped fades > 15 dB 


The record of microwave fading over the two-week period of the experiment 
is summarized in Fig. 19 0 

It was necessary to modify the flight pattern of the instrumented 
aircraft because of the unexpected incidence of deep fading. The origi- 
nal plan ("Flight Pattern l”) called for the plane to fly an H-shaped 
pattern, taking two vertical profiles to study ray propagation, and a 
horizontal profile in between to be used as data for the scintillation 
analysis • 

Because the presence of deep fading made extra vertical profiles 
desirable, a second plan ("Flight Pattern II” ) was implemented, in which 
the aircraft flew a sawtooth pattern between the transmitter and the re- 
ceiver . 

The flight patterns are also summarized in Fig. 19. Figures 20 and 

l 

21 are typical N-profile records; flight 4 (Fig. 20) was made during a 
type B fading period and flight 12 (Fig. 21) during a type D fading pe- 
riod 0 


C • Analysis of the Data 

The primary objective in analyzing the data was to identify the 
mechanism responsible for the deep fades and, if possible, to determine 
whether this mechanism will be active at the 8 km altitude at which the 
pres sure -ref erence-levei system will operate. The method chosen was to 
raytrace through the ref ractivity profiles (N-profiles) recorded by the 
aircraft and to establish correlations between raytracing results and 
signal properties. 
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FADING ?¥F(E FLIGHT PATTERN 


SUMMARY OF FADING DATA 

HAWAII, JUNE 1970 9.6 GHz 
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C:>5dB with ito'5 "v" shaped fades >15<1B 
D: >5dB with over 5 "v" shaped fades >15dB 
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Fig. 19. RECORD OF MICROWAVE FADING AND INSTRUMENTED AIRCRAFT FLIGHTS 
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Fig. 20. TYPICAL N-PROFILES TAKEN BY THE AIRCRAFT ON A QUIET DAY. The nearly vertical lines are 

high-frequency residuals that result from applying a low-frequency filter to the adjacent N-profiles. 
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To raytrace through the Hawaii data, a raytracing program was re- 
quired that could accommodate a nonsymmetric atmosphere with data points 
at varied spacing (variable resolution data)- A similar program had been 
developed several years ago, called RAYTRACE, based on Snell's Law (see 
Chapter II). This early program was capable of raytracing through a non- 
symmetric atmosphere (one in which horizontal gradients were nonzero), 
with relatively low-resolution data. 

To adapt RAYTRACE to the Hawaii data, a new mode of ray tracing was 
required. RAYTRACE began tracing at an arbitrary "center" of a raypath, 
where the ray was assumed horizontal, and then traced its way out, first 
to the left and then to the right, stopping when it reached "satellite 
altitude/' For the Hawaii analysis, it was necessary for the raytracing 
program to begin at a "transmitter" and then follow a raypath through to 
the neighborhood of the "receiver, " each located on a "mountain" at a 
given range and altitude. 

It was desirable, also, for the program to raytrace from satellite- 
to-satellite through an atmosphere whose lower portion consisted of the 
Hawaii data. The variable-resolution feature of the program was needed 
for this phase of the analysis because it was necessary to supplement 
the Hawaii data at altitudes above 4 km with relatively low-resolution 
Standard Atmosphere data. 

Details of the resulting programs, HAWAII (s tat ion-to-stat ion ) and 
THRUWAY (satellite-to-satellite), and the original RAYTRACE are described 
in Appendix C, with a listing of THRUWAY in Appendix B. 

The N-profile data were analyzed in four stages. 

(1) Each available vertical profile was expanded spherically 
symmetrically, and the HAWAII program (stat ion-to-stat ion, 
raytrace between two mountains) simulated the performance 
of the propagation experiment. 

(2) A small but representative sample of these profiles was 
also expanded spherically symmetrically and used in the 
THRUWAY program (satellite-to-satellite) to simulate the 
performance of a satellite system transmitting through 
spherically symmetric atmospheres whose vertical profiles 
were the same as those obtained in the Hawaii experiment. 

(3) The vertical profiles from five of the airplane flights 
were combined to form nonsymmetric atmospheres (atmo- 
spheres with nonzero horizontal gradients). 
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(4) These nonsymmetric atmospheres were used with both the 
HAWAII and THRUWAY programs. 

Details of these analyses are presented in the following chapters. 
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Chapter V 


RAYTRACING STATION-TO-STATION THROUGH SYMMETRIC HAWA I I PROFILES 

"it a man ’will begin with certainties, he shall end in doubts; 
but if he will be content to begin with doubts, he shall end 
in certainties." 

— Francis Bacon 

The first stage in the raytracing analysis of the Hawaii data was 
to expand each of the refractivity profiles generated by the airplane 
into a spherically symmetric "atmosphere, 11 and ray trace through this 
stratified profile using the mountain-to-mountain raytracing program 
HAWAII (see Chapter IV). 

The program traces raypaths, beginning at the transmitter and con- 
tinuing to some point at the range (150 km) of the receiver, thereby 
simulating the configuration of the Hawaii experiment (Figs. 16 and 17). 
The starting angle of the ray at the transmitter is varied by 45 mrad, 
in 1 mrad steps, thus generating an image of the raypaths for signals 

leaving the transmitting antenna across its main lobe. 

>v 

The HAWAII program output is a computer-generated plot of the ray- 
paths between the two stations (at 3.05 and 3.35 km altitudes and sepa- 
rated by 150 km in range), and plots of received power and integrated 
phase defect (excess path length caused by bending and retardation), 
both as functions of altitude at the range of the receiver (see Figs. 

22 to 37). If the program detects a multipath, it will determine the 
magnitude and phase of the resultant phasor at the range of the receiver. 
This program was run on an XDS Sigma-5 computer, and the plots were drawn 
on-line on a 10 in. Calcomp plotter. 

A total of 77 refractivity profiles were available for raytracing. 
The profiles were transcribed from analog tape onto graph paper. Values 
of refractivity were read onto computer cards at 40 altitudes between 
6500 and 13,000 ft (corresponding to 2 and 4 km). These cards were used 
as input to a data-conversion program on an IBM 360/67 computer which 
gave 40 values of refractivity as output, interpolated at 50 m intervals 
between 2 and 3.95 km. This interpolated data formed the input for the 
ray tracing program on the XDS Sigma-5 computer. 
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As seen in Figs- 26 and 27, the refractivity profiles made by the 
aircraft produced only relative values of N 0 To generate realistic 
values of excess path length caused by retardation, absolute values 
were required. The excess path length resulting from bending, and the 
shape of the raypath, are both functions of the gradient in N and, 
therefore, can be computed given only relative values of refractivity. 

The value of N at 4,0 km was taken arbitrarily as 0.0, and a 
standard value of refractivity was added to each of the 40 points on 
the profile. This standard value was determined from the temperature 
and pressure values at 4.0 km in the U.S„ Standard Atmosphere Supple- 
ment, 1966, 15° N Annual profile [9], along with a value of water-vapor 
pressure obtained from a standard profile [10], using Eq. (2.2) to de- 
termine N given the values of T, p, and e. 

As a result, the refractivity data progressed through several stages 
of transcribing, reading, and interpolating before it was incorporated 
into the raytracing routine. Visual comparison was made between the 
data sent from Boulder and the final refractivity profile to ensure that 
the input presented to the computer was reasonably accurate in that no 
gross features had been overlooked or spurious features introduced along 
the line. 

As described in Chapter IV, the received microwave signals were 
divided into four fading types, according to the amount of fading pres- 
ent during a 30 min period : 

Type A: peak-to-peak fading range <5dB, one profile 

Type B: peak-to-peak fading range >5dB, 49 profiles 

Type C: 1 to 5 V-shaped fades > 15 dB, 12 profiles 

Type D: >5 V-shaped fades > 15 dB, 15 profiles 

The atmospheric effects generally held responsible for radio fading 

are (1) multipath propagation in which more than one path is available 

to the radio wave between the transmitter and receiver, resulting in 

destructive interference at the receiver, and (2) defocusing in which 

the radio signal is "spread out" over the transmission path (the loss 

2 

is greater than the 1/R loss associated with spherically propagating 
waves ) . 
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The results of this portion of the ray tracing study are listed in 
Table 17, in which the incidence of multipath and/or defocus ing is com- 
pared to the signal fading type. (For this study, defocusing is defined 


Table 17 


FRACTION OF PROFILES EXHIBITING MULTIPATH OR DEFOCUSING 


Flight No. 

Fading Types 

A 

B 

C 

D 

Multi- 

path 

Def ocus-^ 
ing Only 1 

MP 

DF 

MP 

DF 

MP 

DF 

1 



1/2 






2 



— 


1/2 




3 



— 

1/2 





4 



0/2 






5 



0/2 






7 



1/1 


0/1 




8 







1/3 

1/3 

10 



— 

1/2 


1/1 



11 



— 




4/4 


12 



— 




6/6 


13 



1/6 

2/6 





14 

0/1 


0/1 






15 



' 


2/2 




• • ' -■ v ! 16 •• 



3/4 

1/4 





17 





0/1 


1/1 


18 



2/2 






19 



3/4 

1/2 





20 



1/5 






21 




1/3 

2/2 




22 



3/4 


1/3 




23 



1/3 


1/3 




24 



3/6 






Total Count 

0/1 

0/1 

19/49 

7/49 

6/12 

1/12 

12/15 

1/15 


Fraction Profiles 

Fading Types 

A 

B 

C 

D 

With multipath 

,0.000 

0.388 

0.500 

0.800 

With defocusing only 

0.00 

0.143 

0.083 

0.067 

With multipath and/ 
or defocusing 

0.00 

0.531 

0.583 

0.867 


^Refers to profiles lacking multipath but exhibiting defocusing; defocus- 
ing is also present when multipath is present. 
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as a power variation greater than 15 dB when measured as a function 
of altitude at the receiver range.) There appears to be a definite 
correlation between the amount and intensity of fading and the exis- 
tence of multipath or defocusing. Although both multipath and defo- 
cusing are present in slightly more than half of the type-B profiles, 
they are present in 87 percent of the type-D profiles. It appears 
likely, therefore, that multipath and defocusing are responsible to 
some extent for the fading observed in the Hawaii experiment. 

The "shape" of the fades in the record illustrated in Fig. 18 
suggests that multipath is at least partially responsible for the 
fading; this shape very strongly resembles the result obtained by 
adding two rotating phasors, as illustrated in Fig. 38. The distin- 
guishing parameter in Fig. 38 is the relative amplitude of the two 
adding phasors. The two rotating phasors whose amplitudes differ 
by 0,1 dB or less will produce a "fade" that is very similar to the 
deeper fades in Fig. 18. This suggests that the observed fading in 
the Hawaii data may have been caused by rotating phasors of very 
close amplitude adding at the receiving antenna (a multipath situa- 
tion) . 

It is important to understand the limitations of this portion 
of the study. First, we are treating the atmosphere between the 
Hawaiian mountains as spherically symmetric which, of course, is 
not the case. The validity of this approximation is uncertain for 
the scale sizes considered. 

Second, the refractivity data taken by the aircraft represent 
a "smoothed" value over a considerable horizontal distance (the 
plane flies essentially horizontally, not vertically); therefore, 
the "vertical" profile of refractivity is not a vertical profile in 
the sense of a rawinsonde profile being vertical. In addition, the 
measurements go through various stages of reading and transcribing, 
as described above, before they reach the Sigma-5 computer. As a 
result, although the refractivity profile probably is a good repre- 
sentation of the Hawaii atmosphere on some gross scale, individual 
features of a small scale are necessarily eliminated. It is uncer- 
tain how important these features may be in explaining the behavior 
of the microwave signal but, if they exist, their presence is not 
inferable from the data, so their effect can only be hypothesized. 
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Third, it must be emphasized that the raypaths obtained by tracing 
through this "smoothed" profile are "average" raypaths (the actual ray- 
paths between the two mountain tops in Hawaii probably never looked ex- 
actly like the raypaths drawn by the computer). in addition to being 
spatially smoothed, the refractivity profile also represents a time 
sample because it takes the airplane a finite amount of time to generate 
the profile (on the order of half an hour). It is futile, therefore, 
to attempt to match any specific feature in the microwave signal, such 
as a particularly deep fade, with any one raypath configuration. In 
addition, because the raypaths are essentially an average picture of 
the situation that existed in the atmosphere, it is quite reasonable 
to treat the entire raypath pattern as a whole and not be overly con- 
cerned with multipaths or defocusing that occur very close to the re- 
ceiver altitude, A slight vertical shift in the refractivity profile 
was found to change the position of multipaths and focuses by a much 
greater amount than the actual shift in the profile. 

The results of the raytracing program indicate that phase defect 
is a very strong function, not of path configuration or of path length, 
but of altitude at the range of the receiver. In other words, even in 
a multipath situation, the signals that arrive at a particular altitude 
at the receiver range tend to have phases that are very nearly equal, on 
the scale of several centimeters (see, for example, Fig. 33 or 37), 

At first, one might be tempted to attribute this phenomenon to 
Fermat's Principle but, in fact, this is not the case. According to 
this principle, the optical length of an actual ray between any two 
points and P 9 



is shorter than the optical length of any other curve that joins these 
points and lies in a certain regular neighborhood of it. (Here, n is 
the index of refraction of the medium and ds is an elemental distance 
along the raypath.) "Regular neighborhood" refers to a region that may 
be covered by rays in such a way that one (and only one) ray passes 
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through each point [11]. Fermat's Principle, therefore, specifically 
excludes media that give rise to multipath situations. 

Figures 39 and 40 illustrate the dependence of phase and power on 
altitude at the receiver range and are expansions of Figs. 36 and 37 
in the region of 3.4 km. Note that the phase difference between the 
multipath signals is on the order of a wavelength. Figures 41 and 42 
plot the power difference and phase difference between the multipath 
signals as a function of altitude at the receiver range. The curve in 
Fig. 42 has a slope of 0.1 wavelength/m in the region between 3.402 and 
3.405 km and an average slope of -0.3 wavelength/m between 3.405 and 
3.4085 km. Figure 42 shows that a vertical shift of 6 m in the antenna 
position, or in the index-of -ref raction profile relative to the two an- 
tennas, causes the phase difference to vary only 1.2 wavelengths. 

Inspection of the amplitude record of the microwave signal during 
periods of deep fade shows that incidents of deep fade (<-30dB)are oc- 
casionally interspersed with an incident of signal enhancement. It Is 
much more common, however, for incidents of deep fade or enhancement to 
be clustered in groups, which indicates that, if the received signal is 
the sum of two rotating phasors, the phasors tend to oscillate relative 
to one another rather than swing through full circles of 360°. This is 
consistent with the situation illustrated in Fig. 42, in which the phase 
difference between two multipath signals is less than a wavelength. Note 
that the phase of the detected signal may vary considerably relative to 
some arbitrary value, while the phase difference of its two component 
signals remains small. 

Figure 41 shows the power difference over the same region. The 
average slope over the interval below 3.407 km is 0.86 dB/m, and the 
maximum variation over the 6 m interval between 3.402 and 3.408 km is 
~ 4.5 dB. 

The intensity distribution of the deep fades indicates that 27 per- 
cent of the fades are caused by signals with amplitude difference greater 
than 0.8 dB, 51 percent are caused by signals that differ by more than 
0.5 dB, 67 percent by signals that differ by more than 0.3 dB, and 75 
percent by signals that differ by more than 0.2 dB. Thus, 25 percent 
of the deep fades are attributable to the sum of two signals that differ 
by not more than 0.2 dB in amplitude. 
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The results of the ray tracing study indicate that there is a definite 
correlation between observed multipath and defocusing in the spherically 
symmetric profiles and the fading observed in the microwave signals. The 
very brief statistical study of the incidents of deep fade corroborates 
this hypothesis. Clearly, opportunities are plentiful for further statis- 
tical analysis of these data. Thompson et al [83 present some interesting 
studies of amplitude and phase distributions . 
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N-UNXTS 


Fig 0 22 . N-PROFILE FOR FLIGHT 03A1 S AS DRAWN BY THE SIGMA -5 COMPUTER 
The number 03A1 refers to the first ascent made on the third flight 
This profile was taken during a type-B period. 
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RANGE (KM) 

Fig. 23. COMPUTER-DRAWN RAYPATHS FOR PROFILE 03A1. Profile of Fig. 22 was expanded spherically symme- 
trically and used with the ray tracing program to develop this diagram. Box on left « Haleakala trans 
mitter; box on right = Mauna Loa receiver. Vertical scale is greatly expanded . 


ALTITUDE (KM) 



POWER CDS) 

Fig. 24. POWER VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 03A1 . This 
figure illustrates the power that would be detected by a receiver if the 
receiver moves vertically along the right-hand axis of Fig 0 23. 
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ALTITUDE (KM) 



PHASE (M) 


Fig. 25. PHASE VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 03A1 „ This 
figure illustrates the phase that would be detected by a receiver if the 
receiver moves vertically along the right-hand axis of Fig. 23. 
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ALTITUDE (KUl 



N-UNITS 

Fig. 26. N-PROFILE FOR FLIGHT 04D1, AS DRAWN BY THE SIGMA -5 COMPUTER. 
The 04D1 refers to the first descent made on the fourth flight. This 
profile was taken during a type-B period. 
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RANGE (K m 

Fig. 27. COMPUTER-DRAWN RAYPATHS FOR PROFILE 04D1 . 


ALTITUDE (KMJ 



POWER COB) 

Fig. 28. POWER VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 04D1. 
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ALTITUDE (KM) 
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Fig. 29. PHASE VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 04D1. 
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Fig. 30. N-PROFILE FOR FLIGHT 1XD2, AS DRAWN BY THE SIGMA -5 COMPUTER. 
11D2 refers to the second descent made on the eleventh flight. This 
profile was taken during a type-D period. 
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RANGE (KM) 


Fig. 31. COMPUTER-DRAWN RAYPATHS FOR PROFILE 11D2. Note the extensive multipath. 


(n*> 3aniinv 



POWER (OB) 


Fig. 32. POWER VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 11D2. The 
double-valued nature of this curve clearly indicates the presence of 
multipath. 
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Fig. 33. PHASE VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 11D2. 
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ALTITUDE (KMJ 



N-UNITS 

Fig. 34. N-PROFILE FOR FLIGHT 12D2, AS DRAWN BY THE SIGMA-5 COMPUTER. 
12D2 refers to the second descent made on the twelfth flight. This 
profile was taken during a type-D period. 
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Fig. 36. POWER VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 12D2. 
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Fig. 37. PHASE VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 12D2. 
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Fig. 38. RESULT OF ADDING TWO PHASORS. These curves represent five 
values of relative phasor amplitude. Reading upward at the 180° 
abscissa value, the phasors differ by 0.1, 0.5, 1.0, 2.0, and 5.0 
dB, respectively. 
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Fig . 3.9. POWER VS. 'ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 12D2, 
An expansion of a portion- of Fig. 36. 
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ALTITUDE (km) 



Fig. 40. PHASE VS ALTITUDE AT THE RECEIVER RANGE FOR PROFILE 12D2 . 
An expansion of a portion of Fig. 37. 
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Fig. 41. POWER DIFFERENCE VS ALTITUDE AT THE RECEIVER RANGE FOR 
PROFILE 12D2. Derived from Fig. 40.. 
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ALTITUDE (km) 



Fig. 42. PHASE DIFFERENCE VS ALTITUDE AT THE RECEIVER RANGE FOR 
PROFILE 12D2. Derived from Fig. 41. 
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Chapter VI 

FURTHER ANALYSES OF THE HAWAII DATA 

"One good turn deserves another." 

— Gaius Petronius 

. In the foregoing computer analysis of the Hawaii N-prof iles, each 
vertical profile obtained by the aircraft was expanded spherically sym- 
metrically and then used as input for the ray tracing routine; this was 
the most extensive analysis performed on the data. It is possible to 
use these data in other configurations, and these further analyses will 
be described in this chapter. 

A . Non symmetric Ray tracing 

Because the aircraft recorded several vertical profiles on each of 
its runs, it is possible to combine the profiles to construct a single 
nonsymmetric N-prof ile for each flight (an N-prof ile with nonzero hori- 
zontal gradients). Hopefully, such a profile would be a more realistic 
representation of the conditions that existed along the transmission 
path when the flight was made. 

To obtain N-prof iles, the airplane flew two patterns (Chapter IV). 
Flight -pattern I was H-shaped, with two vertical profiles recorded on 
each flight; flight-pattern II was sawtoothed, with several vertical 
profiles recorded (see Fig. 19). It should be noted again that, even 

in pattern I, the ’'vertical” profile is not truly vertical because the 

plane cannot climb at an angle greater than a few degrees from the hor- 
izontal . 

N-prof iles from five representative flights were chosen for nonsyra- 
metric analysis. 

Type-A period: one flight, flight pattern I 

Type-B: two flights, one each of patterns I and II 

Type-C: one flight, pattern I 

Type-D: one flight, pattern II 
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A type-A period is the "quietest'' from the point of view of the number 
of incidents of intense radio fading, and type-D is the "noisiest," 

The vertical profiles were used to construct a nonsymmetric 
atmosphere for each flight. This "atmosphere" consisted of a matrix 
of N-values, with 40 vertical points (spaced 0,05 km apart, as in the 
symmetric case described in Chapter V) and either 31 or 161 horizontal 
points (5 or 0.95 km, respectively), depending on the flight pattern. 

In the pattern-I data, each vertical profile was entered as a column of 
the matrix, the choice of column determined by the location of the air- 
craft when it made the vertical profile. In the pattern-II data, each 
vertical profile was read into the matrix in such a way as to simulate 
the sawtoothed configuration flown by the aircraft. This is the reason 
for the greater number of horizontal points used to represent flight- 
pattern II data (161 horizontal points are required to accommodate four 
40-point profiles in a sawtooth arrangement; as in flight-pattern I data, 
31 points would provide sufficient horizontal resolution.) 

The points in the matrix not occupied by original data were filled 
by constructing simple linear gradients between points of the original 
profiles. A separate horizontal gradient was established for each of 
the 40 vertical levels and, in the flight-pattern II atmospheres, be- 
tween each of the profiles. The result was an N-matrix with nonzero 
horizontal gradients. These Donsymmetric profiles (Fig. 43) were used 
as input to the HAWAII ray tracing routine. The HAWAII program, based 
on a Snell’s Law raytracing method, has the capability of ray tracing 
through nonsymmetric atmospheres . 

The results of this analysis are illustrated in Figs. 44 through 
48. As can be seen, the results of nonsymmetric raytracing do not dif- 
fer drastically from the results obtained in the spherically symmetric 
analysis . 

Although this case of a simple linear gradient is interesting, it 
has certain shortcomings. it is very unlikely that the atmosphere along 
the transmission path in Hawaii consisted of one or two simple linear 
gradients. A more realistic picture would account for small pockets of 
high or low values of N, corresponding to small convective cells. The 
model described thus far is more sophisticated than a simple symmetric 
atmosphere; nonetheless it is a bit too simple to be real. 
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Another problem is that the N-data used to construct the vertical 
profiles prior to computing the horizontal gradients probably will con- 
tain features that look like horizontal layering but are, in fact, spu- 
rious. These "layers” are generated when the airplane hits an updraft 
or a downdraft and suddenly rises or falls. This vertical motion ap- 
pears in the N-profile as a sudden decrease or increase in N; however, 
it merely represents the fact that the aircraft was not climbing or de- 
scending at a uniform rate. When linear gradients are constructed from 
this contaminated data, the spurious layer "stretches" ouf over scores 
of kilometers when, in fact, it might have a much smaller horizontal 
extent if it exists as a layer at all. 

A statistical approach was applied to overcome these difficulties. 
Each of the vertical profiles was smoothed with a 0.1 km lowpass filter 
to remove the effects of spurious layering, and these smoothed profiles 
were used to construct linear horizontal gradients,' as described above. 

A random-number generator was employed to construct a matrix of 
random variates, with the same number of rows and columns as the N-ma- 
trix. This random-variate matrix was smoothed, first horizontally and 
then vertically, to eliminate any large gradients, and then adjusted to 
a mean of 0 and a standard deviation of 1.00. Each row was multiplied 
by a factor obtained by taking the rms value of the high-frequency com- 
ponent of the N-profiles for each flight.* The resulting matrix of 
smoothed, adjusted, and scaled random variates was added to the N-matrix 
constructed from the smoothed N-profiles. 

Figures 4S and 46 are samples of the N-profiles from this matrix, 
before and after the random variates were added. Note the similarity 
between the N-profiles from the statistical atmosphere and those from 
the simple linear gradient atmosphere (Fig. 43). The difference, of 
course, is that the deviations from a smooth profile are regular, con- 
sistent, and relatively large-scale horizontal deviations in the case 
of the linear gradient; in the statistical case, the deviations are 
irregular and of a smaller horizontal scale because of the random na- 
ture of the added term. 


*These high-frequency profiles, two of which are illustrated in Fig. 

20, were provided along with the N-profiles by the Boulder group. 
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The results of ray tracing through these non symmetric statistical 
profiles (Fig. 47) are very similar to those obtained from the simple 
linear gradient and, by extension, to those obtained from the simpler 
spherically symmetric expansions of the vertical profiles. 

The similarity of these three cases illustrates the small effect 
that horizontal gradients have on ray propagation. Equally intriguing 
is the fact that, even in the presence of horizontal gradients, when 
a multipath situation exists, the excess path length (phase defect) 
appears to be dependent on the endpoints rather than on the raypath 
itself. Figure 48 presents some examples of the phase vs altitude 
record for nonsymmetric atmospheres that yield multipath. 


B. Ray tracing Satellite-to-Satellite 

An operational occultation system will transmit signals between 
two satellites rather than between two ground stations. To determine 
how an occultation system would respond if it were to transmit through 
the N-profiles measured in the Hawaii experiment, the satellite-to- 
satellite raytracing routine THRUWAY was adapted to accommodate the 
Hawaii data. The Hawaii N-profiles provide data only between 2 and 4 
km above sea level. Because THRUWAY requires data from 69 km down 
to the minimum altitude of the ray, it was necessary to "extend” the 
vertical range. 

The U.S. Standard Atmosphere Supplement, 1966, 15° N Annual tem- 
perature profile [9] combined with a standard high-humidity water-vapor 
profile (the high humidity curve in Fig* 53) generated a "standard" 
profile of refractivity vs altitude (Fig. 49), and the Hawaii data 
were then merged into this profile. In the nonsymmetric data, the 
horizontal gradient existing at the 2 km level was carried through to 
the lower altitudes. The atmosphere was considered to be spherically 
symmetric above 4 km. 

Four representative profiles (one from each of the "fading types") 
were chosen for symmetric raytracing in the satellite-to-satellite mode. 
Results of this experiment are shown in Fig. 50; only those portions of 
the raypath below the 16 km level are illustrated. The three parabolas 
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indicate the 0, 2, and 4 km levels, (Original data existed only between 
the 2 and 4 km levels . Between O and 2 km and above 4 km, the measurements 
were obtained from the Standard Atmosphere Supplement, as described above.) 
The satellite altitude is assumed to be 1110.0 km, and the transmitter is 
located 31.98° from the vertical, as measured from a position corresponding 
to the center of the diagram. 

These diagrams illustrate the relative preponderance of multipath and 
defoeusing, even in those cases that correspond to the very quiet well- 
behaved results obtained when raytracing in the stat ion-to-stat ion mode. 
Apparently, multipath is an unavoidable phenomenon in the satellite-to- 
satellite configuration, at least for those rays with closest approach 
altitude in the 2 to 4 km range. (N.B.: the pressure-reference-level 

system will have a closest approach altitude in the 6 to 9 km range.) 

In addition to the four symmetric cases, the five nonsymmetric cases 
considered in the stat ion-to-stat ion mode were also analyzed with the sat- 
ellite-to-satellite program. The results are shown in Figs. 51 and 52. 

The satellite altitude and position were the same as for the symmetric 
profiles. Once again, the measured data existed only between the 2 and 
4 km levels; however, horizontal gradients at the 2 km level were carried 
through to the lower altitudes. In addition, because data exist only for 
a horizontal range of approximately 150 km and the ray travels approxi- 
mately 500 km horizontally in the lower 3 km of its path, it was necessary 
to "extend" the Hawaii data horizontally. This was accomplished merely by 
"flipping" the atmosphere at each end of the 150 km range, changing nega- 
tive gradients to positive gradients (and vice versa), but maintaining the 
statistical properties of the profile. 

Once again, the preponderance of multipath and defocusing can be seen, 
as in the symmetric case. Similarly, as in the nonsymmetric station-to- 
station raytracing, there appears to be very little difference between the 
results obtained from the linear gradient profile and those obtained from 
the "statistical profile." 
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a- Flight 2: flight -pattern I, type-C fading 

Fig. 43. NONSYMMETRIC LINEAR GRADIENT N-PROFILES. These vertical samples were generated from the data 
accumulated by the instrumented aircraft. Linear gradients were used to fill in points not represented 
by original data. For clarity, each of the succeeding profiles is shifted slightly to the right. 
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Fig. 43. CONTINUED 
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e. Flight 16: flight-pattern II, type-B fading 


Fig. 43. CONTINUED. 
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44. COMPUTER -DRAWN RAYPATHS FOR THE NONSYMMETRIC PROFILES IN FIG. 43. 
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h. Flight 4 
Fig, 44. CONTINUED. 
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Fig. 44. CONTINUED 
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Fig. 44. CONTINUED. 
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a. Flight 2 : flight-pattern I, type-C fading 


Fig. 45. NONSYMMETRIC SMOOTHED N-PROFILES. Generation of these samples was the same as for the profiles 
in Fig. 43, with the exception that the aircraft-originated data were run through a smoothing filter 
prior to generating the linear gradients. 
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a. Flight 2: flight-pattern I, type-C fading 

Fig. 46. NONSYMMETRIC STATISTICAL N-PROFILES. These samples were generated by adding a random variate 
to the profiles in Fig. 45. 



eGo-ii-ms 





SEL-71-059 100 


4- 0 



Fig. 46. CONTINUED. 
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d. Flight 11 



Fig. 47. CONTINUED 
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Fig. 47. CONTINUED. 
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a. Flight 11, linear gradient (see Fig. 43d) 


Fig . 48. PHASE VS ALTITUDE AT THE RECEIVER RANGE FOR TWO NONSYMMETRIC 
PROFILES. 
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bo Flight 11, statistical gradient (see Fig. 47d) 


Fig » 48. CONTINUED . 
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a. Flight 2 

Fig. 52. COMPUTER-DRAWN RAY PATHS FOR SATELLITE-TO-SATELLITE PROPAGATION THROUGH NONSYMMETRIC STATISTICAL 
PROFILES. The lower 16 km of the raypaths generated by transmitting between two co-orbiting satel- 
lites through the nonsymmetric statistical profiles of Fig, 46. 







Chapter VII 


IMPLICATIONS OF THE HAWAII EXPERIMENT FOR 
AN OPERATIONAL OCCULTATION SYSTEM 

"By a small sample we may judge of the whole piece." 

— Miguel de Cervantes 

The Hawaii microwave signal was characterized by periods of intense 
fading which, in an operational system, could result in loss of lock and 
subsequent loss of data. The Hawaii test results, however, were obtained 
at a relatively low altitude (3 km, as opposed to the operational closest 
approach altitude of 6 to 8 km) and over a tropical ocean where water- 
vapor content could be greater than over virtually any other terrain. 


A . Water-Vapor Adjustments 


To develop a more realistic picture of how an operational occupa- 
tion satellite system would behave on a worldwide basis, a sample of the 
Hawaii test data was adjusted to simulate conditions typical of regions 
with lower water-vapor content than that of the Hawaii atmosphere. Al- 
though information was not available from which to state categorically 
that the anomalies in the measured N-profii.es responsible for multipath 
were generated by water-vapor inhomogen it ies, it is certainly not un- 
likely. An attempt was made, therefore, to correct for the effect of 
water vapor by smoothing the N-profile over five vertical points (equiv- 
alent to 250 m) . The smoothed profile was then subtracted from the un- 
smoothed profile, and the resulting high-frequency residual was assumed 
to be caused by water-vapor inhomogeneities. Of course, water vapor 
also contributes a low-frequency component to the N-profile but, for 
this study, only relative values were important. 

To simulate different water-vapor conditions, the high-frequency 
residual was multiplied by a factor and then added to the smoothed pro- 
file. This factor was determined by assuming that the Hawaii data rep- 
resented a typically "wet" atmosphere and then by taking the ratio of 
the water vapor on a typically "dry",, day and a "normal" day to that on 
a "wet" day. The values for "wet, " "dry, " and "normal" water-vapor con- 
centrations were obtained from the standard curves shown in Fig. 53. 
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Fig. 53. STANDARD WATER-VAPOR PROFILES. The two curves on the left 
are the same as those in Fig. 10. 

Note that the "dry" and "normal" curves correspond to the profiles in 
Fig. 10 and were derived from the same source. It should be noted also 
that the curves were derived from data applicable only to the middle- 
latitudes [10], where water-vapor concentrations are generally lower 
than over tropical seas, so the assumption that the Hawaii data repre- 
sents conditions for a typical "wet" day is probably conservative. 

The data sets chosen for adjustment were the same as those used in 
raytracing satellite-to-satellite in a symmetric mode (see Chapter VI). 
The four data sets (one for each fading type) were adjusted as described 
above and used for raytracing in the satellite-to-satellite symmetric 
mode. A brief examination of the results (Figs. 54 and 55) and those 
for the original data (Fig. 50) indicates that, as the high-frequency 
component of the N-profile is decreased, the incidence of multipath is 
slightly reduced. Figure 56 presents the smoothed and unsmoothed pro- 
files for the type-D example which is the "noisiest" of the four. 
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Fig, 54. COMPUTER-DRAWN RAYPATHS FOR SATELLITE -TO-SATELLITE PROPAG 
ADJUSTED HAWAII PROFILES. Lower 16 km of the raypaths generated 
orbiting satellites through symmetric N-profiles that have been a 
conditions for a "normal" dav . 
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FLIGHT 11A1 



d. Profile 11A1 


Fig. 54. CONTINUED 
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Fig. 55. COMPUTER -DRAWN RAYPATHS FOR SATE LLITE-TO-SATE LUTE PROPAGATION THROUGH SPHERICALLY SYMMETRIC 
ADJUSTED HAWAII PROFILES. Lower 16 km of the raypaths generated by transmitting between two co- 
orbiting satellites through symmetric N-profiles that have been adjusted to simulate the water-vapor 
conditions for a "dry" day. 







There is reason to be concerned with whether this smoothing process 
actually eliminates the effect of water-vapor inhomogeneities in the N- 
profile. Too much smoothing could remove just those effects that should 
be revealed by adjusting the data; too little smoothing could result in 
unduly pessimistic results.. To assess the efficacy of this process, the 
smoothed N-profile of Fig. 56 was "inverted 1 ' by assuming a dry -atmosphere 
and a standard -atmosphere pressure curve, resulting in a temperature pro- 
file that would generate the smoothed refractivity profile if water-vapor 
content were zero. " ■ 
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Figure 57 is the temperature profile corresponding to the smoothed 
curve in Fig. 56. It is clear that this equivalent temperature profile 
is not realistic and that a five-point smoothing, although obviously re- 
moving most of the effects of water-vapor inhoniogeneity, does not remove 
them all. The results of this analysis, therefore, are rather conserva- 
tive. If the curve were smoothed over a greater range so as to remove 
more of the water-vapor contribution, the multipath incidence would de- 
crease further. 

This point becomes significant when the next logical step is taken 
and the behavior of the system at 3 km is compared to that at 7 km. By 
following the "wet" curve in Fig. 53 from the 3 km to the 7 km level and 



TEMPERATURE (°K) 

Fig. 57. EQUIVALENT TEMPERATURE PROFILE FOR SMOOTHED (5 -point) N~ 
PROFILE 11A1 . 
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by taking the ratio of water-vapor concentration at 7 km to that at 3 
km, the N-profile data can be adjusted, as described above, to simulate 
the water-vapor concentration at the higher altitude. 

When raytracing satellite-to-satellite through this simulated high- 
altitude profile, the incidence of multipath is greatly reduced; in fact, 
in all four of the profiles tested, only eight satellite separations with 
multipath were recorded, as opposed to 60 multipath recordings when the 
profiles were run at the lower altitude without adjustment for water- 
vapor content. These high-level multipaths are caused by the unrealis- 
tic temperature profile in Fig. 57, which resulted from the insufficient 
five-point smoothing of the N-profiles. It appears, therefore, that the 
results of the Hawaii experiment produced a somewhat pessimistic picture 
of probable system performance. 

It should be noted that there is a possibility that a system operat- 
ing at the higher closest approach altitude (6 to 8 km, as described in 
Chapter III) might be affected by signals being deflected upward from a 
Hawaii-like atmosphere at a lower level (the antenna beamwidth will be ' 
such that the lower levels and the surface of the Earth will be illumi- 
nated). When the results of raytracing through the Hawaii profiles are 
examined, however, it can be seen that no signal is deflected to the ex- 
tent that it would be detected by a receiver placed at the operational 
satellite separation (assuming that the greatest separation is 7888.0km, 
as described in Chapter III). Although the possibility of interference 
from a low-level multipath does exist, none of the simulations indicates 
that the probability is significant that this situation will occur fre- 
quently or pose any real problems for the operational system. Should a 
low-level multipath occasionally interfere with system operation, the 
raytracing results indicate that the amplitude of the unwanted signal 
would be significantly less than that of the "direct, " or desired signal 
and, therefore, would have little or no effect on system performance. ' 

B . Other High-Altitude Effects , 

There are other reasons to believe that the results of the Hawaii 
experiment give ^ somewhat conservative picture of the performance 
to be expected from an operational satellite system. First, multipath 
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in the satellite-to-satellite mode will not necessarily have the same 
effect on the system (long and frequent deep fades) as it had in the 
station-to-station mode. Figure 58 plots the excess path length (phase 
defect) as a function of satellite separation; in other words, it plots 
the phase defect that would be measured by a receiving satellite located 
a given geometric distance from the transmitting satellite. 

Once again, as in the station-to-station mode, phase defect is a 
function of endpoints rather than raypath. Phase defects of two signals 
received at the same satellite position due to multipath are likely to 
match to within 20 or 30 cm. Because satellite phase-path results are 
required only to approximately 50 cm, two such signals will be indistin- 
guishable from one another (they are the same within the accuracy re- 
quired by the experiment, as described in Chapter II). 

If two or three of these signals are received at the daughter satel- 
lite, the deep fading associated with the station-to-station mode miRCht 
occur also in the satelli te-to-satellite mode. Inspection of the satel- 
lite-to-satellite results reveals, however, that 5, 7, or more signals ap- 
pearing at the receiver is not uncommon; in fact, this situation exists in 
49 percent of the resulting multipath configurations. Under this condition 
the signals add as random phasors and, rather than resulting in deep fades, 
the power received is increased by several decibels. Because the phase 
defects are within a few tens of centimeters of each other, regardless of 
which signal is locked, the desired information is obtained. This multi- 
signal multipath, in which more than three signals with comparable power 
appear at the receiver, is extremely rare in the station-to-station mode. 

Second, there is a time factor that must be taken into account. The 
Hawaii receiver and transmitter were, of course, stationary . They sat still 
while the weather blew gently past them. The duration of a deep fade in 
the Hawaii data was typically much less than 30 sec, and the maximum wind 
velocity was on the order of 50 km/hr. If fading is ascribed to inhomoge- 
neities traveling into the raypath with the wind, ah estimate of the dura- 
tion of a deep fade in the satellite-to-satellite configuration cars be made 

The satellites, in effect, move past the weather with a velocity on the 
3 

order of 25 X 10 km/hr; therefore, the ratio of 'Velocities" betwesw the 
station-to-station and satellite-to-satellite modes is roughly 0 . 002 . 
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a. Flight 14; linear gradient b* Flight 11; linear gradient 


Fig. 58. PHASE AS A FUNCTION OF SATELLITE SEPARATION FOR FOUR NONSYMMETRIC PROFILES. The phase of the 
received signal appears as ‘it would be recorded by the receiving satellite at the given separation from 
the transmitting satellite. The vertical scale shows separation in kilometers, with 7900 km at top 
and 8050 km at bottom. The scale is arranged this way because the greater separation corresponds to 
the lower closest approach altitude of the ray. The horizontal scale shows phase, in meters, with 
200 m on the left and 1600 m on the right. 




The fading duration could be .decreased by a similar factor, which would 
mean deep fades with a "duration 1 ' of perhaps 0.06 sec. 

Assuming that this rough calculation is significant, the crucial 
question is the number of 60 msec fades that can be expected per second. 
The prevalence of multiple-ray multipath in the satellite-to-satellite 
ray tracing suggests that these periods of deep fade are rather uncommon 
and occur far less often than in the case of the stationary Hawaii re- 
ceiver and transmitter. 

' A third factor is the magnitude of the phase defect experienced by 
the satellite and ground systems. A typical phase defect in the satel- 
lite-to-satellite mode is on the order of 140 to 240 m (corresponding 
to approximately 8 and 6 km closest approach altitudes, respectively), 
whereas the phase defect experienced in the stat ion-to-stat ion mode is 
closer to 40 m; thus, there is a factor of approximately 5 between the 
phase-defect magnitudes of the satellite and ground systems. 

This becomes significant when one realizes that the satellite sys- 
tem is not only sweeping through the weather at a much greater speed 
than is the ground system but also that, as it sweeps, the difference 
in phase defect is likely to be greater by a factor of approximately 5. 
This extra factor further reduces the expected fade duration from « 60 
to « 12 msec . 

To the extent that the Hawaii results can be extrapolated to the 
satellite-to-satellite configuration, there appears to be little reason 
to be pessimistic as to the potential dangers of fading caused by multi- 
path in an operational microwave occultation system. 
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Chapter VIII 


SUMMARY AND CONCLUSIONS 

"it is better never to begin a good work than, having begun it, 
to stop.” 


— Bede 


This research has investigated the feasibility of using a two-sat- 
ellite occultation system to supplement the infrared sensor system (SIRS) 
by proyiding an accurate altitude reference that will serve to fix, as a 
function of height, the tempei'ature-pressure profile of the SIRS instru- 
ment , 

The final stage of the error analysis used simulated SIRS data to 
determine the behavior of the system over a realistic range of inputs. 
Based on this analysis, the microwave occultation system was found to 
be capable of establishing a 300 mb pressure reference level to within 
approximately 24 m. The effects of water vapor can be corrected by 
suitable climatological profiles. This work has shown that even the 
most simplistic climatological approach results in a dramatic improve- 
ment in accuracy; more sophisticated methods promise to be even more 
effective. 

Improvements in the accuracy of the SIRS instrument will yield im- 
provements in the performance of the occultation system. Most of the 
error in the output of the system (neglecting the component attributable 
to water vapor) was caused by error in the input temperature profile. 

As these input errors are reduced, system performance will be signifi- 
cantly improved. Even with the present SIRS instrumentation, however, 
the microwave occultation system produces results that are acceptable 
to, and badly needed by, the meteorological community. 

Ground tests of a microwave link similar to that which would be 
employed in a pressure-reference-level system were described. These 
tests, conducted in June 1970, consisted of a continuous two-week mon- 
itoring of transmission between two mountain tops in Hawaii, with peri- 
odic flights of an instrumented aircraft along the radio path to measure 
the refractivity profile of the atmosphere. 


PRECEDING PA® BLANK NOT FILMED 


143 


SEL-71-059 



The microwave signal was characterized by periods of intense fading 
Using the refract ivity profiles measured by the aircraft, an extensive 
computer analysis was performed to determine the probable cause of the 
observed fades. Results indicated that multipath interference, originat 
ing from anomalies in the N-profiles, was primarily responsible for the 
fading . 

Multipath is a low-altitude phenomenon assoc iated . with water-vapor 
inhomogeneities. Extrapolating the results of the Hawaii experiment 
(performed at an altitude of ^3 km) to the higher altitude of an • opera- 
tional occultation system (w7km), indicated that multipath was greatly 
reduced. it was concluded, therefore, that the fading observed on the 

r 

Hawaii link will have little or no effect on an operational occultation 
system . 

A microwave occultation system appears to be a feasible and useful 
technique for acquiring an important piece of data. Further work in 
this area should be directed toward the operational aspects of this 
system. For example, although this study has demonstrated that water- 
vapor effects can be corrected by the simplest of climatic methods, it 
would be of value to have a more accurate estimate of water vapor along 
the raypath. Several techniques have been proposed, such as passive 
microwave sounders or attenuation measurements of the signal. Tank [12] 
has suggested a dual-frequency occultation system that would determine 
water-vapor density by matching attenuation of the two signals to the 
water-vapor absorption line at 22.3 GHz. 

Another area deserving further work is the inversion program that 
shifts the temperature-pressure profile to the proper altitude. The 
technique used in this study was a simple curve-fitting routine. The 
program can be speeded up somewhat by taking into account the climatol- 
ogy of the geographic area being sampled. Empirical orthogonal func- 
tions may be of some use. 

Although the Hawaii study was intended to provide useful informa- 
tion for the occultation system design, it has generated some scientific 
results of its own. For example, it would be interesting to investigate 
how the periods of deep fade correlate with large-scale weather systems 
that were in the area at the time. It has been demonstrated that fading 
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is probably the result of multipath; however, no attempt was made to 
define the actual mechanism that produces the particular N-profiles 
responsible for the multipath. Presumably, water-vapor inhomogenei- 
ties are the cause, but this study has not attempted to correlate 
such factors as atmospheric conditions, wind, and sea state, in order 
to describe the physical phenomena that may be associated with the 
inhomogeneities. 

It should be pointed: out that the recent studies by Pirraglia and 

Gross [13] of errors inherent in the occultation technique seem not to 

’ ' ' ) . i 

apply to our proposed system. As they state in their conclusion, two 
orbiting satellites in a common plane will not be affected by the er- 
rors described in their paper. This is noted here to avoid possible 
confusion on the part of readers familiar with their work. 
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Appendix A 


THE PRESSURE -REFERENCE PROGRAM- 


l. 

?. 

i. 

•5. 

7. 

8. 
9. 

!•>. 
u. 
12 . 
1 3. 
1 V. 
15. 
15. 
1 7. 

n. 

I p. 
?: . 
21 . 
2’. 
2 ?. 

24. 

25. 
2f . 
27. 

2a. 

2Pi 

33. 

21. 

32. 

33. 

34. 

35. 
3ft. 
37. 
3*. 
39. 
4.). 
41 . 
42. 
5 ?. 

44, 

45, 
4ft. 
47. 

46, 
49. 

*5°. 

51. 

5?. 

53. 

54. 

55. 

56. 
f-7, 
53. 
59. 


f 

C 

C*d9 f *.»*»*j»4i»«*««*(.»*»«»*t»4*e9*>*4*t*tM*t**»4***».*» w********************** 

C 

c 

/ / LUO Wl ft JOB (N529, 240.3.21 ,• STEVFN G. UNGAR* 

// Fit EC FlIRT HCO.PARM. FOP-! - 'CPT-Z.HAP* 

//FORT.SYSlN d') * • 

c ‘ 

C********** i. 0C :t.0.a»**« o tr OO5 .** ***0 •«■**** ** *#•♦**<*<>**■»**»«*****'********* ♦■**** 

C j : ‘ . 

c lud*!g--a program rn fix a pressure level in the amusphere 

c ! 

C MEANING CF ARRAYS: 

c pin is thf no pressure levels listed in the si*s data; Tin is the 

C ARRAY OF CORRESPONDING TEMPERATURE VALUES "MEAiUKc3“ St THE SIRS INSTRUMENT; 

C TPI II IS THE ARRAY OF TEMPERATURE ERROR S ; THE "TRUE" TcMPERATljRE, 

C TU)-TP(n. IS STORffi IN TRIM; 2(11 IS THF ARRAY OF ALTITUDE VALUES FOR 

C EACH OF THE 100 PRESSURE LEVELS; PW (I) HOLDS WATER VAPJR PRESSURE 
C VALUES FOR COMPUTING TH* PEFRACTIVITV PROFILE; cwli.Jl HOLDS THE TwG 
C STANDARD water VAOQR PROFILES (HIGH AND LOW); cP«(ll GIVES RAWtNSONDE VALUES 

C OF WATER VAPOR PRESSURE MEASURED AT PRESSURES iP«41l; APWfH IS EQUAL TO 

C IPWUI, ft J T IN REAL. RATHER. THAN INTEGER , FORMAT i TTIII AND PP ( ll ARE 
C STANDARD ATMOSPHERE VALUES OF TEMPERATURE AND PRESSURE; Ml I I IS THE REFRAC- 

C TIVITY AO ft 4 Y ACTUALLY USED IN RAYTRACING; NU [ I IS THE “MEASURED" ATMOSPHERE 

C ft c FRACTIVITY PROFILE; NN(II IS THE "REFERENCE" A TMD SPHcRfc M-PRQFItEt Pill) 

C IS THE ARRAY OF STANDARD PFgSSiJPES FOR WHICH PktSSJkE rttIGHTS ARE TO BE 

C COMPUTED; TLU» AND PLwIIl ARE CORRESPONDING VALUcS OF TtMPERATURE AND 

C wAtfC VAPOF PRFSSUPFj 5(11 AND 6(11 ARE CORRESPONDING HE I&HtS , FOR THE 

C "F F PER ENCE ATMOSPHERE AND "MEASURED" AtMOSPHEkt-. KESPtCT 1 VfLY ITHE LATTER 
C VALUE MEASURED P° I OF TO SHIFTING!; W<I I . X 1 1 I ) . Z i ( * | . V l (1 ) , 1 1 1 I » ,N2 ( I) AND 
C EKWII! APE WORK VECTORS; LATIIl, LONG( I I . STN&NEU). AND STNOIIJ ARE USED 
C FOR IDENTIFICATION OF DATA SETS. 

C 

C MfANING CF SUBROUTINES: 

C AT SO ANCi All AP r A SORTING AND AN I NTFRPUtAT I JN ROUTINE IN THE IBM 
C SCIENTIFIC SUBROUTINE PACKAGE (SSP). SEE IBM MANUAL H20-U2O5-3 FOR DETAILS. 
C S E 1 3 IS AN SSP SMOOTHING ROUTINE. ATBILO, N, KA V Tk A, ST GEN. AND TALLY ARE 

c described und c r thfir RESPECTIVE LISTINGS. 

C . 

C ' 

UAL PlirO! ,TU Ou >,2(1001 , TT(9I |, Pw( 1 0 J ) , W ( 1 jj I ,X 1 1 13 j > , Z 1 1 100 > 

DIMENSION T Dll 0D> , TP I 100 ) , L AT I 2 > , L ONG ( 2 I 

DIMENSION Y1 (100) , Erfl 25 ,2 I , Pll « » , TL I 91 , A(9t,U(9J 

DIMENSION PLW(O). IPWI151 ,5Pw( 1 5 1 , E HW( 2 5 1 , APWl 1 5 1 , 1 1 ( 25 1 

RFAl*8 SS ,STN4MFf fl!,RE, STNOIIOI 

REAL PP(Q11 .Mill •MNIlQOI .N2U0I. ITYPE .NLIIOUI 

common 2 . pp. TT, NL , MAX 

COMMON /A/ PE.IH 
COMMON / R 4 Y T/ N 

c 

C INITIALIZATION OF CONSTANTS 
C 

DATA PL /5U2. ,400. ,300. ,200. ,14U. , 1U0, , 70. ,60. ,3U./ 

L cV EL S= 9 

DATA PI /3.14jSR?6/,4H/637R.15/,eR/6356.77/,LtVcLS/^/,Pw/lU3*0,0/ 

F=1.0-?R/AR 

C 


A.’ 1 . 

ft). 

n. 

63. 

64. 
15. 
6ft. 
67. 
6 8 . 


C HEAD OATA COMMON TO AIL DATA SETS 
C 

READ (5,11)1 (P( 11,1=1, 100! 

WRITE 16,116) IP! I). 1 = 1 ,1001 
WFIGHT =1.0 

C IG IS THE CONSTANT IN THF EXPONENT OF THE HYPSOMt TRIC FORMULA 1 
G= 28. ft 644* 9. 3 366 5/8.314 32 
BEAD 15,102) PP 
WRITF (6,102) PP 


Ceding 


PAm bla kk not FILMED 
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69. 

79. 

71. 

12 m 

73 . 
74. 
76. 

76. 

77. 
76. 
79. 
BO. 
81. 
82. 
8 3. 
«4, 
95. 
86. 

87. 

88. 
89 . 
99. 
91. 
97. 

93. 

94. 

95. 

96. 

97. 

98. 

99. 
100 . 
1 01 . 
102. 
1C 3. 

104. 

105. 

106. 
1C7. 
108. 
U'9. 
1 10 . 
Ill . 
) 12 . 

113. 

114. 
1 1 8 . 
116. 

117. 

118. 
U9. 
120. 
121 . 
122 . 

123. 

124. 

125. 

126. 

127. 

128. 
129. 
137. 
131 . 

132. 

133. 
i 34. 

136. 
l>b. 

137. 
1 33. 
179. 
140. 
1 41 . 

Uo 


READ <5,1021 rr 
WRITf <6,1021 TT 
Dfl 920 J=l,2 

R FAD <5,114} <EW< !»J>, 1-1,81 
WRITE <6,1141 (FW(I,J),I=l,9) 

00 920 1=9.25 

C <0.5476 IS THF INVERSE OF THE WATER VAPOR SCALt HEIGHT) 

9 20 t W( I,J)=EW<9,J)*EXP(-.5476*I I -8 1 I 

WRI TE <6,n4M<FK(l,J», 1*1 ,25 J, J*l,2> 

C (MIX IS THF NUMBER OF PRESSURF LEVELS IN THE I.mPUT DATA U30U H IS THE 
C STARTING HF1GHT F HR THE ft AY TR AC ING ROUTINE*, cftS IS T He ERROR IN SATELLITE 

C POSITION; F P PH IS THf MAXIMUM PEPMI SSAB.LE FRROR IN PHASE DEFECT ; S,S IS THE 

C SATELLITE SEPARATION.) 

BEAD (5,100) Mf X,H,£R5,ERPH,SS 
W P I T E < 6 , 103) MIX,H,EPS,ERPH,SS 

1 TOP = .7 

ERS-fBS/7481.0 
SS 1 $S / ? 

I H= l 

c 

C PFAD INPUT OttTA F^p 5PFCIFIC DATA $ FT 
C 

27 READ <5,U7,ENn = 999> STNO 
WRITE <6,113) STNO 

READ (6,117) 14T,1£1MG,MJNALT,STNAME 
WRITE <6, 117) LAT,LONG,MINALT, STNAME 
PHI=<IAT<?>/60.D*HT< ! ))*p 1/180 
C (RE IS THF RADIUS nF THF EARTH AT THE GIVEN L4T) 

RE*AR*U.0*< F-3.C7?*F*F )*5IN<PHI)*SlN<PHm 
WRITF <6,119) PH!,RF 
MAX = MIX 

C < DH IS THE AMOUNT QF VERTICAL SHIFT IN KM) 

0H=1 

7300=9 

R E AO <5,1021 I T f i) ,1=1,100) 

WRITE <6,1021 <T<I),l*l,luO) 

READ <5,102) <70(0,1=1,1001 
WRITE <6,1021 (TOO >,1-1,100) 

C l 1 TOP IS THE NUMBER OF DATA SETS READ) 

ITCjP=ITOPM 
DD 47 A 1=1,100 
400 TR < I ► =T ( I I -TO ( I ) 

WRITE < 6 , 102 ) |TM I >,1 = 1,1001 

C < N P W IS TH C NUMBER OF WATER VARQR DATA POINTS. Jw*i UR 2. CORRESPONDING TO) 
C LOW Oft HIGH STANDARD WATFR VAPC1R PROFILE.) 

READ <5,1221 NPw.JW 
WRITE <6, 122) NPW.JW 

READ 15,123) U PW ( 1 I , EPW < 1) , l =1 .NPW ) 

WRIT? (6,121) < I PW< ! I ,EPW< 1> , l=l,NPW> 

C 

c 

C "REFFRFNC5 A TMnSPHE RE ” PORTION OF PROGRAM 
C 

C INTERPOLATE NPW WATER VAPOR VALUES AT 100 INPUT PftfcSSUiUS 
C 

DO 80 1 = 1, NPW 
8U APW<I)=!PW<|) 

PWMIN=APW(NPWI 
DO B1 1=1 .MAX 
IPWMAX-I 
IPWM-lPWMAX-1 

IF (PUl.LT.PWMIN) GO to 82 

CALL AT SO < t> <J>,APW,EPW,W,NPW,l f Zl,Xl,NPWl 

CALL ali <pm,zi,xi,pwm ,NPW,.01, IER) 

IF UfB.Gf.2) WPJTF 16,106) IER, I 
91 CONTINUE 

82 CONTINUE. 

C 

C CALCUIATF RFFPACTJVI7Y AND ALTITUDE PROFILES FUK " kcF t K tNL fc ATMOSPHERE” 

C 

Z<1)=0 

PLl-AlOGlPIin 

TS1 = TR { 3 ) } 
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143. 

144. 
146. 

146. 

147 , 

1 4R » : 
140 . 

1 b'J. 

1st. 
1-*?, 
15 3 . 

164 . 

15 5. 

156. 

157. 
15R. 

1 5° • 
160. 
161. 
16?. 

16 3. 

164. 

1*5. 

146. 

167. 

168. 

16 7. 
17j. 
171. 
17?. 

17 3* 
174, 

176. 
17«-.. 

177. 
17R. 
17*. 

1 Pii. 
18). 
1 “?. 
IP*. 
IP4. 
l p S. 
1H6. 
187. ' 
133. 
1«. 
190. 
1 C 1 . 
19?. 
153. 
19'.. 

1 96 . 
15o. 

1 c 7, 
1^. 
190. 
20 '. 
?n. 

20 
Z--'\ 3. 
?04. 
2 ' 0 . 

?'>7. 
? t a, 

?l 

21 «, 

?1 1 . 
212. 
’1 3 . 
214. 


NN( 1 1=77.6*1 Pill -PW( 1 » 1 /TS1+3. 75E 5/T$l#PW< 1 ) /I SI 
Vi R I T 6 (6,1001 IH.Zm.MNm 
0(1 1 Is?, MAX 
Pl2*ALDGIP<I 11 
T S2= T ft < I 1 

Z ( T I = Z < 1-1 1-MTS1*TS?1/2*<PL1-PL2I/G 

pww=o.r 

C (wATFP VAPna=4 A P n v f. 22 KM, f ALCULA T ED ACCORDING Ij L XP JNfc N T 1 A L 
C BF TWEEN MAXIMUM HEIGHT OF DATA AND 20 KM. I 

IF (I .FQ. 1PVM) Z PW= 7 ( I) 

IF ( Z m.LT.20.0. AND. I. GE. IPVMAX) ; PK( I l = PV( IP«M J*fcXP( - .>476*,( Z( Jl“ 

*/pw nil ’» 

IF 1ZU1.LT. 20. 01 PN=PWII) 

NN( 11 = 77. 6* ( PH I 1 - PV>'I / T S2 + 3. 75F 5/T S2* PW W/ T S2 j 

WP.ITF -16.1 VM; i .211 l.NM U 1 

p 1 1 = P L ? ; 

1 T S 1 = T 5 2 ; i 

C . i ■ 

C JNT r P PJLA TF VALUES OF TEMPERATURE I Til, ALTITUDE (A), *ND WATER VAPOR 
( PRFSSU» i; I PL W ) FUR 9 STANDARD PRESSURES I PL > 

C 

no 2 1 =i ,i fvei s 

CALL AT SO I PI il) , P.TR ,W ,MAX, 1 .Z1 , XI, MAX ) 

CALL Al I I PL U I . Zl.XI, TIM 1 .MAX..01 , IER > 

CALL AT SO (Pldl.P, 7 , W , MAX , 1 , U , X 1 , M A X 1 
CALL AL1 < PLt 11.71 »X1, Al I) .MAX, ,01 tlfR I 
CALI AT SG |Pl III. P.PW.W, MAX, 1,71, Xl.MAXI 
CALL AL 1 IPH M .Zl.Xl ,PLW(1 1 .MAX,. Jl, I£R1 

2 CONTIHUF 

D(l 467 1 =1 .LEVELS 

667 WRITE IMTAI !,»L( Tl ,A( II, TU !|,PLV( I I • 

c <z*u i- is rue hfight rr thf 300 ma lfveli - 

7 300* Al ’ I 

4 WRITE (6.101 I ?30‘ 

1mA Xb »<mAx» 

c 

r Sf T N< [ I *P £'F C ACT'IVI TV r 'F STANDARD OPY AT M jSPHfcKc 
C 

on 4 i s j , qi 

N( I 1=77.6 *Pt»( I |/TT( I | 

4 CONTIHUF 

C 

C INTERPOLATE DATA R F c « ACT I V I T Y PROF 1 L E AT 1 KM INTcKwALS, AMU REPLACE 
L STANUA-Y ATMOSPHERE REF RACTIVTTY VALUES WITH lrH.Sc VALJcS 

r ... 

DO 7 1=1 , 1 W A x 
X = I - 1 

TAIL A T $0 ( X,Z,NN,W,"AX,1 ,Z1 ,X1 .MAX) 

CALL A L I IX, 21, XI, Ml I »,MAX, I.E-4.IEK) 

TF (IFM.0E.2) GO T° 303 
7 CONTINUE 

C 

C SMnPTH PR 3 F I L f I* NF I GH^noHODD OF TRANSITION BETwccN C«T« KcFRACTlVlTY 
f. PRDFlLf AMO STANDARD ATMOSPHERE REFRACT 1VITY PKJHlc IIKANSITIUN 
C OCCUR S WAR 70 KM ALTITUDE) 

( • 

DH y 1*1 ,10 

I I s I v AX-5 ♦ I 

3 j N2in=NUll 

CALL SEls (N2,N?,10,IFR| 

IF ( JFF.Nt.OI GO TO 3< - 4 
f'h 3! 1=1,19 
! I* I^AX-SU 
3i MIlDWm 

c 

C P7YTKACL THROUGH R c c f 0 E Nf F ATMOSPHERE PROFILE 
C 

CALI ’ AY T R A (SS,!"S,H,SR 1 
C (PHASf 0FFFCT.5R) 

C 

C 
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21 5 . 
2 ) 6 . 
217 . 
213 . 
210 . 
2 2 !-. 

22 1 • 
222 . 
??3. 

224. 

225. 

226. 
227. 
223. 
??o. 
2 30. 
231. 
2 32. 
2 3 3. 
2 34, 
? 35 . 

236. 

237. 

238. 
230. 

240. 
? 4 1 . 
24?. 

243. 

244. 

245. 

246. 

247 . 

248. 
240 , 
2 5u. 
251 . 
? 52 * 

253. 

254. 

255. 

256. 

257. 

258. 
257. 
26 I. 
261. 
262 . 
263. 
244. 
265. 
266 . 
?rr. 
26 «. 
269. 
27C. 
271. 
?7 2 » 

273. 

274. 

275. 
275. 

277. 

278. 

279. 

280. 
281. 
2A2. 
233. 

284 . 

285. 

286. 
?B7. 
288 • 
289. 


C "MEASURED AT M 0$PHERF" PORTION nF PROGRAM 
C 

C SFT 2KW=APPR0PRl ATE STANDARD WATER VAPOR PROFILE 
C 

00 87 1—1*25 
EWW(I»*EW( I*JW1 
.3 7 Z Z ( n * i - 1 

WRITE I 6 ♦ 1 1 4 I EWW 

C CALCULATE HtFRACTLViTY AND ALTITUDE PROFILES FOR "MeASUREO ATMOSPHERE 1 
C 

Z ( l 1 = 0.0 

PH*AU)G(P!l n 
TS1=.TU 1 
P W ( L) * F WW ( 1 1 

NL(il=77.6*(PU)-PW(llWTSl*3.73E5/TSl*PW(ll/T6i 
WRIT.E 16. 1001 I U *7(11. ML III 

on a 5 i = 2 , max 

PL2*AL0G(Ptl )>. 

TS? = T( I 1 

Z< IlrZU-1 ) + (TSUTS2)/2*iPLl- 0 L2)/G 
PL 1 =PL2 
85 T S 1 = J S? 

on 86 1=2, MAX 

PW< 1) ='J.C 

IF (7(1 1 * GF .20.0 ) GO TO 86 

CALL ATSG ( MU ♦ZZ,Fvw,w,25,i,Zl t Xl,251 

CALL ALI (Z( 11,71, Xl,»W( II , 2 5 , .01 , I ER ) 

IF (IfiF.r,P.2| WRITE (6,104) IFR.I 
IF IIEP. SO. 31 WR1TF (6,116) ll 
IF ( IER.FQ.31 WRITE (6,116) 1 1 
IF ( IER.EQ.31 STOP 

8 a NL( I l = 77.6*(Pt I )-PW( 1 > l/T | 1 10#T3E5/T I 1 I*NI1)/II 1 1 

00 83 t«l, LEVELS 
C 

C IMFBPCJLATE VA1UES OF TEMPFRATUPE ITU, AL T ITuOG (8), AND WATER VAPOR 
C PRESSURE (PLWI FUR ° STANDARD PRESSURES I PL I 
C 

CALL ATSG IHIH.P, T,W,MAX, l,Zl,Xi .MAXI 
CALI ALKPK I),7l,Xl,TL(I),M4X,.01,IER I 
CALL ATSG (PUU.P, Z»W,MAX,l»Zl»X L .MAX I 
CALL AL.I (PL( 1), Zl.xl ,81 I ) ,*AX,.01,IER> 

CALL ATSG (PIU >,P,PW,W,MAX,1,Z1,X1,MAX) 

call ali (PLm.zi,xi,PLW(n*M4x,.'ji,cefl) 

88 CONTINUE 

C (ZZ300 IS ALTITUDE OF THE 30'J MB LEVEL I 
ZZ3D0=B( 3) 

DO 16 1 = 1, LEVELS 

16 WR I T c (6,100) l,PLin,3<t>,rim,PLW( II 

c 

C MERGE KEFRACTIVITY PROFILE WITH STANDARD ATMOSPHtRE REF RAC 1 1 VI TV 
C °ROF I LE , AND RAYT°4Cf 
C 

CALL AT8ILD 

CALL RAYTRA ( 5S, FR S.H , PDU 
C (PHASE DEFECT* PCM ) 

C 

C SHIFT PROFILE UP l KM (DH*1 > 

C 

C MERGE REFRACT IV1 TV PROFILF WITH STANOARD ATMOSPHERE KEFRACTIVITY 
C PROFILE, AND RAYTRACF 
C 

DO 83 1=1, MAX 
33 ZU)=Z<l>»nH 

ZZ30C =ZZ3C0+DH 
CALL atbild 

CALL RAYTRA ( SS , f R S ,H , PD2 ) 

DPH()Z=P02-PDI 

WP ITE (6,202) OPHDZ 

DPH!=SF-PD? 

C (IF ERROR IN PH AS F DEFECT IS SUFFICIENTLY SMALL , CEASE I TfcR AT I ON) 

IF ( ABSIDPHI ).LE.EPPH) GO TO 91 
C 

C CALCULATE NEW DH , SHIFT °ROF ILF, MERGf WITH STANDARD ATMOSPHERE 
C PROFILE , AND RAYTRACE 
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c ' 

12 DH= DPHI/OPHDZ*^riGHT 

WRI Tf. ’l 6, 201) nn,0PHI ,z IIH) • 

no n 1 = 1 , max .... 

i i z t n = z < i > ♦ oh 

ZZ3D"=Z7 3:>.3+PH 
c ( phase 111 F c CT = pn? J 
C . 

C C A L Cllt. ATE UPHDJ, D^HI 
C 

TAIL AT9IL0 

CALL 0AVTR4 I S S , F R S ,H . Pf)2 ) 

1 S PPh1 = SK-°D? 

f. ( ] F ERF 'IP TN PHA?F DEFECT IS SUFFICIENTLY SMALL* LtAiE ITERATION! 

JF 1 A° S ( f) D H | ),LF«FP p Hl GO TO 91 
f ^ 

C 

f Phi NT PFSULTS 
C 

C ( ZOIF =E RRtJB , I A uc T R R S * AT 3 30 LEVEL I 

r GO TO 12 

e i zotf =1000* t zz3 n C'Z3co » 

rfPJTr (6,1131 SIM] 

write i6«mi znr>c, n ^, zdif 

wP I T E (6,111) .... 

*AX=MIX • 

DH»l 

AVS*Z?3l.A 

77 ' CONTINUE • 

WRITE (6,113) STNU 

WRITE (6,119) STNAHF,LAT,inNG,RE 

C • 

C P* i Nit ALTITUDE ERRORS Fno r ACM OF 9 PRESSURES 

C . 

CALL $TG?M A, R.Pl , ANS, ! TOP,tEVFlS ) 

WM Tr (6,111) 

GO TO' 21 
939 CONTINUE 1 •" 

C 

C CAL Clll AT P ANO PRINT STATISTICS FOR ENTIRE RUN, "„GUNT" IS AN ENTRY 
C IN »*STC,€N.*' 

c. 

WKlfF (6,120) I TOP , SS 
CALL COUNT ( 1T0P,LFVFLSI 

STOP . * 

3u3 WRITE (6,107) IFF,! 

STOP 1 

3u* WRITE (6,1>9» IFp 

• stop • " 

I DU FORMAT (I4.5F9.3I 

101 FORMAT ( 3FS.3 ) 

I I ? FORMAT <lnc«;?j 

1 j 3 FORMAT |»y3Dt* **9. PRESSURE LEVEL: CALC. VALUt = » ,F8. J,» KM. TRUE 

* VALUE = *,FR.3,» KM.*//» DIFFERENCE =',F12.3,' METERS*) 
l'J4. FORMAT (• ERROR IN ALT » F W . I F« = ■ , K. * I = * , I 4 ) 

l.'E FORMAT (* FFRCIR IN ALT ,1400. IER=»,U,* I = ',14) 

1 06 FOR MAT (' ERROR IN ALl.PW. !£R**»I4,' I 

K-7 FORMAT (' ERROR Hi ALT . 1FR=«,I*,' l=',I4> 

UP FORMAT (' »,?FR.3) ' r 

H 9 FORMAT (• FRPOR IN SE13.N. If P« • • 

110' '- FORMAT ( 9 F 3 « 3 ) 

, 1 1 V FORMAT ( • 1 ' ) 1 ' 

112 FORMAT ( 4 1 SX , F 1 ? . 6 1 .) ' ' 

113 F Oh. MAT (I OAR I 

114 FORMAT ( 0 F 1 0 , 6 ) ' ! 

116 • FORMAT ( * * , Q F.R« 3) 

117 FORMAT I 2 ( 2-1 3 , 1 X ) , I 5, 7 AR , A5 I 

US. FORMAT (• LAT«* ,F 10-4.;* RahTaNS KE=»,Flu.4 t * KM. * I 

113 format |* C * ,7A°, AE3//2X, 213,* N ',213,* W RADIUS UF EARTH', 

*' AT -THIS LATITUDE' ,F]0. 4,' KM.') 

120 FORMAT ( » ISTATTSTICS FOR*, 14,* STATIONS'/' SAlcLLlTc SEPARATION' 

*.Flo.4,' Km.ij < 

121 FORMAT (I 5, FID. 6) 
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3 66 . 
767. 
3 f fc. 
3 6 9. 
■*7f . 
371 . 
377, 

373. 

374. 
773. 
37^.. 
3 77. 
3 7H. 
37«. 
3 1>0. 

361. 

362. 

383. 

384. 
395. 
3P6. 

39 7. 

388. 

389. 

390. 
3°1 • 

3 9?. 
708. 
394, 
3«5. 
306. 
397. 
399. 
399. 
4C^. 
41 ' 1 . 
47?. 
4G3. 
■40'.. 

40 8. 
40,6. 

417. 

4 r » . 

4t 9. 

41 «. 
411. 

m;. 

41 3. 
414. 

418. 
414. 
417. 
4)8. 
41°. 
4?i. 
<*2 1 . 
4??. 

423. 

424. 
424. 
4 24. 
4 ? 7. 
423. 
4? ° » 
4?:,. 
431 . 
43?. 
433. 
4 34. 
436, 
4? t. 


12? POP mat {2131 

123 FOPMAT (14.P11.61 

2CC friCMAT <• lDS0H=*,t2l .9 1 

?,1 FORMAT (4F70.3) 

2v-2 fDKMIT (* ,627.9 l 

(-MC 


SUOFOUT T N c iTPIlO 


C This SU* 0 OUT INF MfoG c S THE "MEASUPEO A Tmq SPHt kb" i^-PRJF ILt WITH A STANDARD 
C ATMOSPHlCE RFFPACTIVITV PROFILE 

c l 

REAL 7(10f I ,PP(91 1 ,TT (Qll ,N2( l' , l,N(91 I ,PV.(lu^l 
REAL T{ lO.i) ,PUOol .ZUIVI ,XlUO3),NL(l0(-l 
REAL WC1 A 0) 

P F A l * 3 C c 

C.CMHiTM Z , P° »T T » til » **AX 
COAM IM /A/RF,ih 
COMMON /RAYT/N 
I MA X= J ( M A* 1 
nc 12 I « 1 , 7J 

-12 N( I l=77«6*PP< I )/TT{ I 1 
J M A K= MA X 
PC 3 I « 1 H ( I MAX 
X = I-1 

CALL AT SG ( X, 7 , NL . W , M A X , 1 , 2 1 , XI , JM* X ) 

CALI A L I (X.Zl ,Xl ,N( I I, JMAX, 1.P-4.ICR 1 
I M It R*G C • 2 ) WRITF (6, 100) 1E*»1 
IF IIEP. OF.?) wo I TR (6,102) 21 

C woJTf (6,108) I , X ,N( I 1 

3 CONTINUE 

DC 4 1=1,10 

I 1= |«ax-*M 

N?UI--N(!U 

CALL 8=1 i (N2,*-?, P, IFR 1 
IF (1RR.NF.0) W P 1 T c (6,10l> JEP 

nri 8 1 = 1,10 

I I = I v ax-8 ♦ i 

8 nuu--n?(ii 

RETURN 

1C.) (• cuppq T M ALI. 

101 FORMAT <* ^PROf IN S613. 1FR=',IS1 

1 j ? FORMAT (« ' , B F 1 0 » 4 ) 
r. o format ( 1 3, c 9.4,8x,Fi?. ei 
FNfl 


0 F A L FUNCTION N<“3(H 


C THIS SU8S03TINF IS fALLFO BY RAYTRA. IT I HfE RP0L4 1 ti IN THL REFRACTJVITY 

r afkay ai a give*' value of c (distance from th c Cc'xTt* jf the earth). 

c 

M AL*8 R 

REAL N2< 31 I 

COMMON /PAYT/N? 

c CM MTN /a/Rf,lH 

«a=R-pc 

1=1 NT l G. A > + 1 

IF U.LT.l) WRITE (6,1) R A 
! = MAX'M 1,1) 

P AP =A 1 NT ( o A 1 
J*U1 

I = M IN.M 1,91 1 

J = M IN 3( J , 9 1 ) 

N*N?( I |+(P«-«AP|» (M2( J)-N2< I) » 
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407. 

43“. 

* 35. 
4 4“. 
4*2. 
44?. 
443. 

4 4/, . 
44«j. 

446. 

447. 

4 4 (■ . 

44 9, 

4 c ''■ . 

4*1. 

45*. 

4T3. 

45*. 

45s. 

454. 

447. 

453. 

459. 
4(_v. 
4t 1 . 

4 1 ? • 

46 1. 
4>.*4 . 

4 7.4, 

47 6, 
46? . 
tfeS, 
4M>. 
47 j. 

471. 

472. 
471. 

4 74. 

475. 

476. 

477. 

478. 

4 79, 
4?<1. 
4*11. 

40?. 

483. 
4P4. 
4E5. 
400. 
407. 
4P«. 
■469. 
4 90. 
4*1. 

44>. 

4*3. 

454. 
4Q5. 
496. 
457. 
499 . 

499. ' 
506. 
5C1. 
5u2. 
5P3. 
504. 
57 5. 

506. 

507. 

500. 
57 9. 


N = "1*1.00-6 ♦ I.l^ 
fifr'ic/j 

1 FORMAT (» RA LE 5 S THAN PF. KA» # ,F15..7I 

FMr 



SU wc-, '.in IMF PtVTRf I ACL ,ACC*H,«»MQFT1 

C 


C THJS IS l f c OMAT* S OhlNtlPtF QAYTMClNG PLMJTINfc Hc.» l5 Y uR . JOSE 
f PC^AL tkf 4-0I A7. ' ACl IS HALF THfc SATELLITE STPAnATluN, 1 \* KM; ACC IS 

C TH“ e ?*.>4 IS SAT C U ITF SEPARATION; H IS A START iNu VmLJc f OK THf 

C CLOSFST .APPROACH A l T | T| jr>* OF T HC RAYS PHOFT IS I Ht RfcSJLTJNG PHASE 

C DEFECT. : T M p PkOOFA" 0 C G INS AT He I GMT H, : R A YTKACt 5 TO SATELLITE 

C ALTITUUC. A*fO TH^-i ADJUSTS H SO TH4T fHf RESULT*^! 5 ATfcLL I T c HALF- 

C SPACING EQUALS 1 *CL T* RffHlN ACC. ; 

C 

Q f 6L* J rC^S.^STM, rv.TAIN.OSORT.OAffSIN 

RF AL* 0 P 1 t 07 t M, C, *N,RN t nb LR,o ,RR. DELS, DlLSN.OlL *. Ri »P l » SET, GAM 
BF-U* 1 ' NP ,TTHp, xi ,Y1, X2,v?,THP, TH,THI , T*<2 , Y, Ru . KUP t Ot LT , ACL , OELTH 
“ t 'UH PI LTD I ,t U 1, PWU001 ,PU9» 

RFH** PMOF T , E ♦ L 
REAL*P93, C ‘F,7 
CDPMt N / A /R F » I H 
P 1= 3. 141 59269 359 
c R ? = ’ ■• 

, E«l=G 

F W * 0 

0? LtM *r.’.UU05 
RHP - U >u. ■"> 
r0P*- , C'P + 63 70 . 3 5 
o I s oz *h ' 

£ = -.&47F-3 
11=0 

if i n.Fj.oj «i*Ri*i.rnc 

KKO = (- 
Klu-J 

ni'ii 

HR = P I - R f 
II 3= II 
r?*4 

IF ( U.6T.*) r?= 3 
6 j RMA<=ct40 7.9 
SFP= 100.0 
SFN--IC0.0 
FACT 1=1.0 
*80 

STUl=0.f> 

Y =D6l TH 
K»G 

2 3 STl.iaO*UD*' 

Cl-STU 

4 C =0.f>0t! 

AtTN*C-.'J 
AT TN? *0. *> 

AT06 = 0. ' 

AT 

AT WV6*Q. 0 
AT 0=0.0 
ATWV'C# 0 
ATw=O.D 
CAT1 *0.0 
SN=o.ono 
TH*0.9F>0 
CONST = P c 
I WDIC F=0 
R ‘PI 
J »o 

Yl =F I — C. ?. 5 
Y 2 = R 1 40 . 25 
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51,'. 


DNDP = V f ¥2 I-NIY1 ) 

511. 


0NDP= •:*. “■ /'». 4 

5*2. 


RHP*M 4 1 1 / "‘NOR 

5)3. 


RN=ISM o n»e 1 

514. 


RHU=a35(«H^> 

515. 


RR = R HO—R I 

51* . 


GA^=FR*OSIM °I-DFLTH l/RHO 

517. 


GAIDARS I N f G AM | 

515. 


BFT*DFLT H-GA M 

519. 


DFlR=RHn*OSl N( 3FTJ/0SIMPI-DEITM 1 -R I 

52c . 


9F L 4 = 3HQ*fl F T 

621. 


R = 9*r)FLR/2« n 

522. 


R K I 

5 ? 3 . 


NP=NI o 1 

6?*. 


GO TH 51 

526 . 


CHNT INiJF 

52b. 


NP*N( P ) 

52?. 


DFLP=R*DSQRT( (R*NP/KN) **2-1.000 l*DELTH 

5 29. 

20 

RK=R 

529. 


J = J*1 

53C. 


R = P *0 tl R / 2 

531, 


NP*N ( R ) 

522. 


0ELK = R*DSQPT1 ( 1*»<P/RN 1**2-1.001* DfclTH 

533. 


ofls=r*r*np*pflth/rm 

534. 

51 

CONTI NUF 

535. 


th=th*delth 

536. 


IP 1 II1.GT.P! GO TO 69 

537. 


OEl SN= { hP-1 .l" DO 1 *DFL S 

538. 


0FIC*9? LS-l PSQRTI Ott S *DF IS-OE LR*DFLR 1 *UCOS ITH 1 *U.l * *JS 1 N< Trt M 

539. 


H= R - R F 

5 40. 

1 8 

CONTINUE 

541. 


c*c*nuc 

542. 


SN*5N*OFLSN 

543. 

69 

CONTINUE 

544. 


OELtH»Y 

645. 


R-»RP 

54b. 


R=8*05LR 

547. 


If (fc.GF.FMAX) GO TO 6 

548. 


GO TO 5 

549. 

6 

N2*IN01CE-1 

55U. 


IF < INDJCE.NF.OI GO TO 7 ' > 

551. 


R1 = R 

55?. 


TH1 =TH 

553. 


1 NO 1 C F* INOICE + I 

554. 


GO TO 5 

555. 

7 

IF (N2» 8,9,9 

556. 

8 

t F-'O IC-F * I NO I C F ♦ l 

557. 


GO TO 5 

558. 

9 

R2 = R 

550. 


TH2=TH 

569. 


XI = R 1 *0S I N< TH! | 

561. 


Yl«Rl*OCOSiTHn 

562. 


X2=R2*DSIN( TH2 » 

563. 


Y2=R2*DCOS(TH2J 

564. 


TTHP=-( Y2-Y1 )/< X2-XI » 

565. 


THP»DATANITTHP| 

566. 


Xl*( Yim*TTHP|*OCOS< THPI/ROP 

567. 


XU(050RT1-X1*X1«1.00CI»/X1 

563. 


TH c OA TAM Xl)*THP 

569. 


X]»R0°*0SIN(TH1 

57C. 


C=C*(X1-X2)*(1 .C'OU/OCOS ( THP )-l,ODOI 

571. 


Y1=0SURT( P0P**2 — Xl* XI I 

572. 


TTHP*ACt./<OSORT ( ROP*R0P-ACC**2 »1 

573. 


THP*OATAM TTHP) 

574. 


IF J 1K.GT.R»,0R.(F4CT1 .LT.P.1M FACTl-1.0 

575. 


IF (K.LT.49I GO TO 55 

576. 


K*0 

577, 


Il-o 

578. 


S T f J= C 1 

579. 


Kfi=0 

580. 


Sf P*100.0 

581. 


SEM-100.0 

582. 

55 

CONTINUE 
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583. 


IF (STU1.GT.0I K=0 

58*. 


IF (ST'Jl.GT.O) 11=0 

585. ■ 


IF (K.GT.O) GO Til 24 

586. 


IF ISTJ.NE.C1 1 GO TO 19 

587. 


IF (I I .GT.OI GO TO 2 5 

588. 


7 = TH 

589. 


Rp = R J 

5 <50. 


*!=PI-1.UD0 

591. 


$TU=1.001 

59?. 


GO TO 4 

'593. 

19 

STU=1Z-TH1/(C0-9I 1 

594 . 


GO TO 26 

,595. 

2 5 

ST1J=F 

' 5®6. 

26 

EP=STU 

597. 

24 

CONT.T MUE 

598, 


E = FP 

59 ). 


EP S = T H“ THP 

601. 


Ip <K8.FO.OI ! 

bV ! . 


X=AHS(8PS> 

tc?. 


$t»=fp 

603. 


STlJ] =STU 

6( 4. 


Ko=K4H 

6C'5. 


18 ( III. GT.OI r,r. TP 68 

6t b. 


Ohr.FTs?00’J.)*(r*SN| 

6? 7. 


• ’ 04?*2.G*0AT) 

665. 


K«ITF16.75> PHO p T i OA 7 

6'. °. 


HR 1T8 (*. ,1Q) 5TL),HP,MK.S C P,SEN. EPS 

61H. 


K)G=510»1 

61 J . 

.8 8 

CPNTINUF 

6 1 2 . 


l c ( X.LF.ACC 1 GO TO 22 

M3. 


f ACT? = ARS( f P5 | 

614. 


8tCT3=A8SlFP| l 

618, j. 


IF 1 f AfT?.GF.FACT3l F AC TJ =FAC T 1*0. T5 

616. 


H*H J -Pf 

617. 


I c (c'PS.GT.M GO TO 67 

615. 


I c ISFN.lT.fPSl HN«H 

il l c . 


IF ( SFS.LT.6PSI SFN-EPS 

6?'-’. 


GO TO 62 

6 ? 1 * 

67 

IF <SFP.GT.FPSI MP=h 

6 22. 


IF < SFP.GT.F d$ ) SFP= E ®S 

6?3. 

67 

CONTINUE 

6 ?4 . 


IF ( JSFP.LT.10F-.01.AN0.I SFN.GT.-109.0M 

6 ? 5 . 


R 3 =R I 

6 76. 


FI *» f-F4CTl*FPS/$t«l 

6? 7 , 


K* K* 1 

6?8. 


GO TP 61 

6? 9. 

oO 

FAC* (HP-HN1 /(SEP-SEN 1 

6 30. 


K*KM 

631. 


F3=« I 

637. 


X P = A3S ( S F° ) 

633. 


XN = A 6SI SEN 1 

634. 


IF< XP.CT.XN1 R I = F.F^HP-S9P*FAC*FACT1 

635. 


IF (XN.LE.XP1 »l =RF*HN-SEN*FAC*FACT1 

r>?6 . 

61 

CONTINUE 

637. 


H = R I E 

633. 


KK«J=KK0*1 

'639. 


IF (XKC-.GT.501 R I =-l 0 . ODO 

640. 


PHOI = PHf)F T 

641 . 


EP1*FPS 

64?. 


DA 1 = 2. 0*PAT1 

643. 


H= R 1 -R F 

644. 


' IF H KKC.LT.49J.PR. IH.LE.O) > GO TO 5U0 

645. 


WR I T E t 6 ♦ 4 1 1 H»l, SFN,HP . SEP 

64 6. 


WR I TT < 6 1 1 4 I 

.647. 


HR I T E < 6 » 151 FPSiSTUl ,H*F4CT1 

64 8. 

5 'jiff. 

CONTINUE 

64 9. ■ 


IF (X.GT.ACCI GO TO 73 

651. 

?? 

CONTINUE 

651. 


IF <Klf. GT.OI GO TO 70 

b52. 


111 = 0 

653. 


GO TO 23 

654. 

70 

CONTINUE 

655. 


MO=ROP*DCOS < T H 1 -R F 

6 86. 


Rr = flO»>-Re 
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6*7. 
65.0, 
6*7. 
66 '. 
6.61 . 
66 1 . 

66 3. 
664. 
666. 
666 . 
667 ■ 

6£.a, 
666 , 
67j. 
67 ) , 
6 7?. 

67 3. 
6 6 4. 
D 7 S . 
6 76. 
677. 
6 76. 
6 7' 5 . 
6 0 .. . 
6*1 . 
6P 2 * 

683. 

684. 
686 . 
686 . 
687. 
680. 

689. 

690. 
69|. 

692. 

693. 

694. 

695. 

696. 
6*7. 

698. 

699. 

700. 

701. 

702. 

703. 

704. 

705. 

706. 

707. 

708. 
705. 

710. 

711. 

712. 

713. 

714. 

715. 

716. 

717. 

718. 

719. 
7? ). 
721. 
7??. 

723. 

724. 

725. 

726. 

727. 

728. 


ThsTH*57. 2957800 
H*R I-Pc 
Ocl H=H-hO 
A T0= 2 • 0*ATO 
ATWV=2.0*ATkV 
ATW=2.fl*ATW 
DAT 1 = 2. 0*OAT1 
AT TN = AT WV 
AT TN2 = 0A T1 

WPITt (6,371 AT TN * ATT M2 
11=113 

oisr» r *.o 
0 E L L * ? f - u 0 • 0 ♦ C 
P0= 2 "i ».0*SN 
PHOF T =OFL l* R D 

I c (tPS.NF.FPH Fkl = ( PHDFT-PHDl 1 /(ePS-EPl ) 

IF ( ( I I.EO.Ct.AND. IEPS.NE.EP1) I Ek=E R 1 

IF l (Il. c Q.9).AN n .(FPS.NF.fPl> I F«2= (CAT1-CA. I/I tPS-cPl I 
0A1 1 =[iAT 1 -E PS* = P ? 

ATTN?=DAT1 

WR I T F. (6*37) ATTN,ATTN2 
pH[)PT = pHnpr-EPS*pR 
Cf. 1 = Ab? ( C 0 1 1 
wRlT=<6,34> c ei,FP? 

WRIT? <6,111 of t TH,ACL* r »c LL ,R0,PHnFT,H,H0,nELrl,ii f ArTM.ul6T 
RETURN 

II FORMAT! * ORB ! T AIT ITUDE «• . FI 1 .4, « KM.'. LUX* 'CENTRAL ANGLE = * , F 1 2 
1.5,' DEG. • .10X, 'ST. LINE PATH-* ,f 12. 5 ,* KM.'./' oENDl NG-* »F9. 3 , ' 

. 2M«' «l 9X , ' RE T AftDAT ION- *.F9.3»* N . • , 17X , • PH ASc OtFfcC T «• .F9. 3. • M. % / 
3' ACTUAL RAY HF IGHT* * .F 8 . 3 » * KM«* ,10X . 'APPARENT kAV HE IGHT- * , FB. 3 . 
4' KM.', 9X,'0FLTA H«*,F9.3,» KM.*,/' 0ENOING ANGL E * • , F 1 i . 8, « OEG. • 
5, 9X, 'SIGNAL ATTENUATION-* ,F8.4.' 08. * * 19X , 'SURFACE D1ST ANCE= »,F9. 3 
6, • KM.* | 

14 FORMAT!* •,*fPS*.'STUl , ,*H*,»FACTl«) 

15 FORMAT ( ' »,F12.7,E16.6,2F12.7> 

16 FORMAT < * *,016.61 

17 FORMAT ( 8016.71 
30 FORMAT I 1 12 1 

35 FORMAT! 2X, • PHASE DEFECT-* ,f 12 .6 ,1 X , » AT TN= • ,F 12 . 6) 

36 FORMAT!' ' , ' ERR- • ,F 1 2 . 5, IX , *€«2- * ,E 1 2. 51 

37 FORMAT!' * , • ATTN6-' ,F 1 2. 7 ,2 X, * ATTNl* ' , Fl 2 . 7 1 

39 FORMAT!* *,6*16.7) 

41 FORMAT!* *,4711.71 

END 
C 
C 

(>*»***•••* *•*••*****«•*«*•*•******»•*•*•»*•*«- ***««•»,«*»* ***«•*»***«•* 

c 

c 

SUBROUTINE STGEN ! A, 8,0, ANS, NUM, LEVELS I 
C 

C THIS SUBROUTINE CALCULATES AND PRINTS THE AVERAGE, S1ANDARD DEVIATION, 
C RMS, MINIMUM AND MAXIMUM VALUES OF ALTITUDE ERkuR AT THt 9 PRESSURE 
C LEVELS. ACTUAL CALCULATIONS ARE PERFORMED BY SUBRjUUNe TALLY. 

C 

REAL A!9I,C! 1 20, 9) , D( 9) , 0 | 9) ,RMS! 9) 

DIMENSION AVER! 91 , S 0! 9 I ,VMIN!9),VMAX(9),R(9,1) 

DATA C/1 C80*0»0/ 

DIF-ANS- 8(3) 

00 1 J-l, LEVELS 
C!NUM,J)-(0!J)*niF-A! JM*1.E3 
R! J,1 )=C(NUM,J> 

1 CONTINUE 

WRITE <6,1071 
00 3 J-l, LEVFLS 
WRITE ! 6, 106 ) 0! J ) ♦ R!J,1I 
3 CONTINUE 

CALL TALLY (R, AVER, SO, RMS, VMIN,VNAX, LEVELS, 1 ) 

WRITE ! 6 , 1091 LFVEIS 
WRITE (6.101 I AVER 1 1 I 
WRITF (6.102) SO(l) 

WRITE (6,1051 RMSI1) 
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‘ ; V?29. 

731. 
■;* 7 32. 
^>.3 3. 
*'•. 734. 
•''-:---73S. 
•'* 73 6 ’. 
V"-737,' 



. 7? 6. 01 



•7.39V 

1JC 

NO & 

740. 

lul 

’•741. 

' 10 2 


742. • 

10 3 


743. 

104 


744, ’ 

•74 5. • 

U-5 

'•••■ loft 


•7;4ft,. 

.. 1 J7 

. -.1 

747 , 

108 


:7A8. . 



•> * 'O',' ’• 


74 - 5 . 
750. 
/ 75 1 ; 
■it 752 . 

753 , 
7 754 . 
^'755. 
4-7.56,.; 
'•4757. 

■' 7 59 •' ' 
■■'. •• 7 . 5 < 5 . 
V/ 760 . 

7. 7.6 1 * 

\V': 7 1 2 '.,. 
?'^763.' :• 
'..•■••764. ■ 
“>7 6 5.7 

r 76 7 ... 

' 76 3. 
- /,!V 769 ." 

: ' 7 7 ). 

■ ’771 

' :' 7 7 ? , • 
4:777. • 
^ 774 . ■ 

77775 '. 

7-776. 

7 77 . . 
" 1 7 78 ..; 
779 . 

-y 780 . 

’■' 761 .. 

782 . 

■,, 4783 .! 

"Vi 784 .- 

■ :' h 765. 
f.tf?E6i ; 
..■* 767 , 

. :.. 78 «« ' 
■-' 7 & 9 * 
79 J. 
7 ' 791 . 
.■-7792 , 

' .' 793 . \ 
' 4794 ... 
s;"-' 7 ' 95 . 

7 79 ft. 

•-'J, 797 *> 

/' 793 . 

799 . 

‘,-' 7300 . 

: :';V 801 . 


- WRITE' (6,103) VMINI 1 7 

• WRITE (6»104|-;;.VWikXm 7'. ' 7 ' ‘ • T 5 *. : 

return ■ .■ .. -V-. - • ; \.\ 

ENTRY COUNT ( NUW, LEVELS) , • 7 ' . . ( 

CALL TALLY <C, AVER, SO, RMS, VM IN , VM AX ,NUM» l E VELS » 

'V- WRITE, (6,108) •, ... 

00 5 J c l.» LEVELS' " 7 : ; 4 -;' 7' - 1 ' 

WRITE (6,106 1: 0<JI , AVER < J I » $0 ( J I , RMS LJ ),VM{N(J 1 ,VHAX( J 1 
-CONTINUE - ‘ 

return ’■ v 

FORMAT <//lJX, ‘STATISTICS 1 FOR IOWER',14,* PRESSURE LEVELS*/ I 
FORMAT (• MEAN* , 6X,ICUF9. 2,2X| | , 

FORMAT (• ST. OEV.* ,2X,lOtf9.2,2X)| • v~.. '■ 

FORMAT- (« MIN*, 7X,10<F9. 2,2X1) -7 

FORMAT (« MAX' ,7X,10( F9«Z,2XI 1 , • >, •' 

FORMAT (• RMS*,7X,1A(F9.2,2X») ' 

FORMAL (* • ,F10.3,ICC F9.2,2xn ••.' •■■ 

:> FORMAT; (//•■ LEVEL (MB1 */l ' i : 

..FORMAT; (//• LEVEL (Mp| * , 5X , * MEAN * ,6X , * S T * DEV.*, 

*5X, *RMS », RX, *MIN* ,8X, *MAX* /» 

END ' 


C , . 

C SEE 
C . 


-SUBROUTINE TALLY ( A , A VER , SO , RMS ,VMI N, VMAX , NO»NV I. 

ST.GtN • FOR - OFSCR I PT I 3W . 7' •••_ . • 7 , <- 7'. 

REAL At 12 0,9 I, AVER « 9 1 , SOI 9 ) , VMI N! 9 > , VMAX t 9 I , IOT ALti 0 1 ,RMS(9J 
00 1 K = 1 , NV , 

TOT AL IM *0*0 . - ■- ‘ 7 ' 

AVERIK )=0.0 

SO(Ki = 3.0 . . 7 . 7 - 

> MS I K» « 0.0'; -• , 

VMINIK I =1 , 0 E 75 / • 

VKA X ( K) *-l . OE. 7 . 5 ; 7 '■ 

CCNTINUE - -..V-.. - •' 7.'. J... : ■- 

Of! 7 J = 1 ,nh 7 

00 7 . 1=1 , nv 77 *■' '■ ' ••: ;• 

AA=A( J, II 

TOTAL! n*TOTALUI*AA 
If (AA-VMIN(II) 3,4,4 

VM1N< I 1 *A A ! 7 7 ‘ ; 4 7; •;'■ •. ;4'-' ‘'' ! .^-7 

IF UA-VMAXt ! Ml 6,6,5 V 1 : 7 :ir.-7' *' '■' > ••■?•• .-V"''-’ 

VMAXt I ) = A A ' 

SOI M = S0II)6AA*AA ■■■}'.■ 

CONTINUE " '•';' J . . ■••••:,.■'.•■ , v * " 

00' 8 1 = 1, NV ' 

AVER m»TPT ALII »/N0 ' 

SOI 1 1 = SORT I A6 St. I SOI I ) -TOTAL! II*TOTAU l l/NOI/CNU-i III 
RMS 1 1 I'sSQRT t AVER I I l*AVER.I I I^SOt II *S0( III 

CONTI NUE •; . , 7 '' 

RETURN ' •; • ‘ 7 

END i ' • . •••' 7'- - ‘ - 'j '7' 
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Appendix B 


SATELLITE -TO -SATELLITE RAYTRACING PROGRAM THRUWAY 


THRUWAY raytraces between two satellites. The "atmosphere" is 
divided into four regions. 


(1) The region between 4 km and RO (69 km) consists of 
spherically symmetric standard-atmosphere data. 

(2) Above RO km, the medium is assumed to be vacuum. 

(3) Between 2 and 4 kra, the atmosphere consists of non- 
symmetric Hawaii data. 

(4) Below 2 km, standard-atmosphere data are again used, 
but the horizontal gradients at 2 km are carried 
through to the lower altitudes. 

Because the XDS Sigma-5 computer does not have sufficient core storage 
for all the Hawaii data, it is stored on a disk (called the RAD, for 
Rapid Access Device) and pulled when necessary. 


l. 

i. 
a. 
5. 
4 • 
7. 
3. 
9. 
1 >• 
11 . 
17. 
n. 

14. 

■is* 
16. 
17. 
1 D . 
19. 
2 J. 
21 , 
??. 
2 3. 
2*. 
25. 
25 . 
27. 
2?. 
? 9 • 
3-:«. 
31 . 
3? • 


C 

c 

t «*******%#»« «****#**«:*t***»*l'.**««******* • **(.*****«*« *#***««**• ******4 ********** 

C 

c 

J.1A 36 *>7. UNGAP 

ALL »6T l MIL .GO) , I FSI2 , | FORMAT »B 1 # 1 5 SI lb , 3C 1 , ( Vt_ I 

co&TRAn Jin 
C 
C 

C 

r. 

C *****thmj*AY*«"»«"!’ 

C NON-SYMMFTRTC a T mo$ ph£ b t 

IMPLICIT Rf Al*«f«»CM? > 

REAL* a 5TNPHJ ,TH,COSPHj t 6‘ , C t SN» SN 1.C1 . PHOF T T K JN ,SL JPt . fcP. THP.TH} 

Rfc4L*3 M*PH»THl,THIfTHPKECf TH2«U»lHn v PSI » rHOI&f^N 
RtAl*A r'UM,r>FLH.%' ; p( ? ,3 r .) 

P£AL*4 B< 531 »?'>?!£«?!,'>*< 93) 

REAL N2(4%20l 
INTEGER ONT.TIMes 
LOG 1C AL HOP AO - 

: LfMO.AL PLF)TR,<3FGJN r PL OT P , D E R I V . PL OT PH, P L ul Dd 

COMMifj.- • TH0IS»THM1.N/FLY/If LY, PLOTR , PLOTP , PL JTPh .PLOT DB 
CCMMC-N- fl? 

COMMON FFP,'1,F02,U8 

OS C'CII))*l/!>Cr>S(U) 

OAPSIN1U1*OATA*)(U/DSOPT J1-U*U1 1 

c 

C INITIALIZE CONSTANTS 

c 
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33. 

34. 

35. 
■*6. 
37 o 

3b. 

. 

(.r., 
4 L. 
43. 

43. 

44. 
4?. 
4 Ci » 
47. 

4 ^ • 
Si' . 
* 1 . 
t 3 . 

63. 

64. 
66. 
*6, 
67. 
firi. 
69. 
6J. 
7 1. 
7 ?. 
63. 
6 <-. 
t 6. 
6c. 

67. 

t-R, 

69 . 

T-J. 

71. 

7?. 

73. 

74. 
76. 

76. 

77. 
7». 
7‘ . 
4 ' . 
PI . 
P?. 
3i. 
44. 
BE. 
46. 
°7 . 
RB. 
99. 
90. 
°J . 
93. 

93. 

94. 

96. 
u fc. 

97. 
9 P . 
99. 

1 0f>* 
l-’l. 
10?. 
103. 
1 04. 
HO, 


C rosi 16 THE I Mr R FM C *'T IN STARTING ANGLE PS1. Kc la THE KA01US OF THE EARTH. 
C R3P IS T He SATELLITE A LT l TUOt = DHIG IS A BIG NUMBER, UstFUL TO HAVE AROUND. 
OPSIO.O'WJl/3 
C ^ L L PI or. SI 1 4 > 

0610= l • PSD 
Rf =6^71 , 

PI?=3. 14159267? 

3np = U10.0 
KDP = hC1 D + c E 
f 

C (-EAO INPUT IT A T A 

r 

C M'l IS THE H=IGHT OF THE "TOP" OF THE ATMOSPHERE, SS IS SATELLITE SEPARATION, 
C PS I IS THt INITIAL STARTING ANGLE, AND MUM IS T«c NUHbi.* UF RAYS TO BE 
C L AIJ'JC Hr'' 1 o NREAP IS A SURP.CUJT INF THAT R 9 AOS THr N “PR JF I L t DATA. K IS THE 
C NUMBER 3 F RAYS FAR LAUNCHED, 

Rf AO ( 6 * GOT 1 RD.THPFEC 
WRITE (6,6091 «0,THPPFf. 

READ |& f *1GI SS.PSI 
WRITE! 6, SS.PS1 
READ (5,lfej) .J'f 
WRITE (6,169) NO* 

CfLL NRfcAD 
CO* «.'■)*« E 
R PS I = P S I 
M=RE*4 0 9 
R?=RF*I 6.0 
K = T 

1 COST I NUF 

ACLl=9.n 

ABC 1=0.0 

L NCONT IS AN ENTRY IN NR PAD. 

IF (K.NF.91 TALL NCONT 

°SI =R PS I 

K«0 

C (EL TH IS THC STEPPING VALUE IN THE A AY TRAC 1NG KsMjTiNt. 

OL L T w = . • .* '.,.-16 

W R l ! e (0,6,1) SS, 90, THPfltC.PS (, OclTH 

C THO IS THE L p F T-HAI.F CENTRAL ANGLE MtASUStJ Al The. SUrfAte OF THE EARTH. 

TMO= «»ARSIN(SE/ ('"P*?! | 

C 

C L-CTERMINF c A Y p A T m PARAMETERS *£IWcE\I TP4NSM I T T I i»o StTcLlITE ANO TOP EDGE OF 
f AT M ')S°Hf PE 
C 

4 CONTINUE 

P FG IN= • T R UE . 

OL = Of MSI TM<T“P«1I 
wPITl (R.4-*| bob, R( ,TH0,»S1 
b A NG 0 = •> np/pn«-nr 
WRITE (6,44) PA. "GO 

43 Fl'KMAT (• • , •P'»P«' t F12o7, • <0*',F12.7,' THO* 0 ,r U o 7 , 

* • PS!=»,Fl2.7t 

44 FTRM.M (» Tt>»7=»,Fl? 0 7> 

thi »-pi ?- pst ♦pars rnRO»/»o*r.ci 

WF I T " (6,40) THI 
4 5 FORMAT (“ TMl e* ,F! 3.7 1 

a R=kO +1 5 IN ( THO* TH 1 ) /9C 
f 

C RAYTRACJ'IG ROUT I HE 

C 

c 

C INITIAL! Zf R AY Tf AC I NG PARAMETERS 
C 

C P A NO TH of-F POLAR Cf ORO I MA Tf S WITH ORIGIN AT THc CcnltR OF THE fcARTFo H IS 
C THt CLPSPST APPFPACH ALTITUDE Of THE PAY. SN IS THc cXCtSS PATH LENGTH OUE 
C TO RETARDATION. C IS ThF EXCESS PATH LENGTH OUc TO bc.NJlNG. PH IS THE 
C ANGLE BETWEEN THE POSITION VECTOR » AND THE RADIUS DF CUKVATURE VECTOR PCM. 

C NSTART IS AN ENTRY IN SU6FPUTTNE NNRFAD, NS cT la AN efriT k V IN SUBROU? 1 nM: 

C NR E AD. 

C no 
PH=*THl^PSI 
WRITE <6,6161 PH 
6X6 FORMA! (« PH=*,F12.7) 
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Us. 
1C 7, 
lUB. 
109. 
12 '?. 
111 . 
112 . 
in, 

1 14. 
11 ?. 
ns. 

117. 
UP. 
1 19. 
12 :.-. 
121 . 
122 . 

123. 

124. 

125. 

126 . 

127. 

128. 
129, 
1 3^» 
131 . 
1??. 
133. 
1 34. 
136 . 
136. 
1 37. 
13*. 
139. 
UJ. 
141. 
I 4 ? • 
143. 
1 4 A . 
l 4 *. 

1 46 . 

147. 
l 4 R. 
1 49 . 
.1 5 ( • « 
1 M. 
152. 

163. 

164. 
J 65. 
156. 
1*7, 
156. 
l e 9 . 
If 6. 
161. 
162. 
1 63. 

164. 

165. 
1 6s. 
167. 
ISA. 
16 7 . 
1 7 i . 

171. 

172. 

1 73. 

174. 

175. 

176. 

177. 

178 . 


R = R0 

1H=TH1 . . . 

h=r o ■ - ' - • 

S N = <;.0DP 

DEI 4 = Q.ODP 

CALI NST ART (TH.THIMDX ) 

cgsphi = ocns (ph> 

SIWPHI = PS l N (PH) ( . 
ftS«'3E*l.60 
CALL USf T 
C 

C MAIN R4YTPACING LOOP 

C VSHIFT IN An ENTRY IN S'JBfl OUT I Nfc NN° E AO. THINUX Is A PARAMETER ASSOCIATED 
C WITH CALLING DATA FRO* ThF DISK. N(P,THI IS A CALL TJ FUNCTION SUBROUTINE 
C M, WHICH COMPUTE S T H e . I Nl»F x nf RCPoaCTION AT PUiJTIlNV IR,TM IF R IS GREATER 
C THAN 4 KM (IP L^SS Than 2 K M. NN(«,THI PEK e nkMS T H t SAME FUNCTION FOR 
C 2 < c <4 KM. 

2 J CONTINUE 

IF (TH.GF.THINOX.ANp. P . (.F . t 3 7* . > t CALL f.'SHIFI (THlNJX) 

OF R I V=. 4 ALSE. 

ROISO.V* 

C 

C COMPUTE no* AND On(jTH« DNDTH*'J ABOVE 4 KM. 

C 

NGRfln = d.LP. 637E.ANO.R.GE.6373 
IF (HGRA')l GO T 0 40 
NI=N(3,TH» 

ONOTH=Q .0000 
IF U.GF.CSI GO TO 

OHDK= (N(K,TH)-N(0 +->OfS, THI J/ROIS 
GJ TCI <i 

4. N|*Yi(H,H) 

rNl>TH = INN(B,TH*THPIS)-'INIf , TH-THDI j) ) / (2 ‘ThOIS) 

5; CONTINUE 

IF (k.r,k.M7E.5.0w.R.LF.,37?.5) GO TO Si 
I F I ,J .L F. 4374 . 0 . AMD.H, GF, 6373.1 I GO TO 52 
1 F (3.GE.4374 ) GO TO 54 

-4 01 S.7-M -NNt B+PDl s. THJ l/(2*pnis ) 

GO T 1 4 

f.t. OM*h= ( iMS.nf IS t TH)-Nlfi*sCIS t TH| | / ( ? » P O ! S ) 

GO TO ?. 

52 FOIS=‘'.r e 

ONiWsf NN(R-HDI S, TH1-NM R*RO| $ , TH) t/(2*RJl S) 

Gj TO s 

51 3WH«I NU -RMS, THJ- ri(P*KpI S,TH) )/l2*R0IS) 

8 C ( N T I NllF 

OND TH=u. )PP 

C 

f CCMPIITF ®f*M 

C SPECIAL CASES EXIST etU? ONG«*o. IN WHICH CASE ROH=udlG. A.< APPPuKlMATION TO 
( INFINITY, A'lO Ffo PH *'• , h HORIZONTAL RAV, IN »HlCH CASE OcL*. THE INCREMENT 

C IN 8 , mijst 3 F- CfMPllTfO ACCORD! NO TO A S tCl \G-ilRjc K APPROXIMATION. VARIABLES 
C T AM) ° arc nUP*Y PAOAMF T fcp S FOR THE PLOTTING RuJTlNt PlUI. 

IF ( OARS I PNOP ) • L F . 1 . Op- ?<i I 0£RlV=.TRUE. 

BN = o«-NI 

IF (C.NF.i) c-o T o 2 
IF ( r> = 9 I v ) GO TO ¥'■ 

R OH = NT /(L.Nn«*CL'5PHl-sjr|PHt*DN0TH/R) 

C-0 TO 61 
f.) 

iMT.: ( f , 73 j ) R.TH 
M CONTINUE 

c = ah * i l.o— ornsi psn ) 

WPITr (6,180) C 

2 CONTI NUF, - 

IF (PH.NF.C.OPO) GO Tn 2? 

IF (PEF I V) GO TO 62 

ROH = N I / ( ONCR ♦C0SPH1 - S 1NPHI *DNOTH/R 1 
GO TO 63 

62 *0H«D9!G 

WRITE (6,70^J R,TH 
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179. 

63 


continu c 

1 8tJ . 



HELP = IR*P8LTH**?f/2*n .ODO-R/BOH) 

181. 



s - p. + PFt» 

1 u 2 » 



TH = TH ♦ OFLTH 

183. 



IF < .NUT. PL'lTR o (jR.p o 0T.R2) GO TD 2'J 

ieA. 



I =R*PSIN( TH) / c 1 

185. 



P=R •OCOSt TH1 -RF 

1 86. 



IF < p .GT ,, 1 2» Co'lPo Po L T, - 9. U» GO TO ?<j 

157. 



r. ALL PLPT (7,T,t*) 

188. 



GU TO 2C> 

1 89. 

c 



180. 

c 

CALCULATE DELPHI, QQH, OfcLSN, OfrlR* AND OtLC 

1«1. 

c 



19?. 

22 


DELPHI = OFLTH*U 0 on(i-R/( C0FPH1 4R0H( ) 

193. 



[F (FT, RADI GO TO ?i 

! 98 . 



OEITH** .roc 5 

185. 



GO Tn 221 

196. 

21 


Pf LTV) =m fl COOl 

157. 

221 


CONTINUE 

193. 



PH * PH *■ DELPHI 

189. 



cosph’ = dchsiph > 

2 Cm* . 



S1NPHI = OSINIPH ) 

20lo 



IF ( TER IV1 GO TO 64 

20? . 



MOM = M/IONHOPCOSPHI o SIMPH! «*ONOTH/R 1 

2(3. 



GO TO 65 

2'*4. 

64 


ROH=P M IG 

21 8. 



WRITE (6,730) R.TH 

7 U6 . 

65 


CUNT! NuE 

20 7. 



DELS = R*nFLTH/COSPH! 

2('8. 



OtLC = PFLS8(l.jro-DrOS«?H-PHI } 

2t.Q. 



DEL R « OFLS*SINPHl 

?!'■. 



OtLSN = ( M I - I.^OO) o DELS 

211 . 



if <n?i s.eo.c.oj G n to 23 

2 12. 

M3 


FORMAT (• **3018.7) 

213. 



TH = TH ♦ OFLTH 

214. 

615 


FORMAT (' TH= *,F1 2.7, • P=*,F12.7I 

215, 



C = C * OFLC 

216. 



r = p ♦ oriR 

?17. 

C 



213. 

c. 

1 F 

THE PAY STRIKES T-t~ SURFACE OF THE EARTH, 

21Q. 

c 



220. 



IF (R.Gt.PF) GO T n 300 

221. 



PS I = PS 1 +PPS 1 

222. 



NUH=NUM-1 

223. 



WRITE (6,103) O.NIJM 

??4. 



GO TO 4 

2 25. 

c 



226. 

c 

CALL PLOTTING ROUTINES 

22 7. 

c 



2?3, 

?v»0 


CONTINUE 

229. 



IF ( .NOT.PLOTP.dRoBoGT.P26 GO TO 5 

233. 



IF (.NIT. REGIN) GO TO 10 

231. 



T=P*0SIM(TH)/R1 

232. 



P«R *i»COS( TH ) -8 F 

233. 



IF I °.GT. 12c.OoOR.PoLT 9.0) GO TO 5 

2 34. 



CALL PtOT<l*T,P) 

235. 



BEGINS. FALSE. 

236. 



GO TO 5 

237. 

10 


CON T T NUF 

2?8. 



T«P *n S| N( TH ) / p 1 

230. 



P«f*.-vc(JS(TH)-RF 

2 4<>, 



IF (P.GT. 12.O.0B. P.LT.“ 9.0) GO TO 5 

241. 



TAIL D LflT l?,T, p ) 

24?. 

5 


CONTINUE 

24.3. 



SW * SN ♦ OELSN 

244. 

c 



245. 

c 

IF 

THE PAY IS AT THE TpP OF THE AT MO SPHERE , 

246. 

c 

•JITH LINF4P INTERPOLATION 

247. 

c 



24 8, 



IF (R.GF.Rfl) GU TO 6 

2 48. 



IF (R.GF.H8 GO TO 20 

?5f . 



THC*TM 

251 . 



H*P 


STOP kmYIRwcING 


STOP RAYIkAclNi ANO P*G;'. ‘riO 
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GO TH 20 


252. 

253. 

254. 

255. 
2 56, 
2 57. 
258. 
2 5Q. 
261'. 
261 . 
262 . 

263. 

264. 
?t 5. 
2th. 

267, 

268 . 

269. 

270 . 
27 ) , 
2 7 ?. 

273 . 

274. 

275. 

276. 
2 77. 
77*. 

277. 
2fi'i. 
201. 
?H?. 
233. 
?H4. 
266. 
2 fib'. 
2fi7. 

2 P o , 
?fiy. 
2<M. 

2- n. 
?.-2 . 
20?. 
2 ^ 4 . 
? c 5. 
2 *:>. 
2 ? 7 . 
2?fi. 
? ot ;. 
3 r 0 • 
3" l . 

362. 

363. 
3‘%. 
305. 

3- .6. 
31 7 . 

30 a. 

i'-i. 

310. 

311 . 

312. 

3 13. 

314. 

315. 

316 . 

317. 
3 1 3. 

313. 

J ?1. 

371 . 

?22i 

323. 
? 24. 
375. 


C fNO Hf CAYTPACING ROijtiNF 
C 

c 

C CALCULATE RAYPATH PACAFFTFRS PY UNFA* 1 NTtF POLATl UN FROM POINT AT WHICH 
C PAY L E A V c S ATMHSPMERF TO POINT AT WHICH IT INTERSECT i SiTELLlTF ORBIT 
C 

23 TH = TH ♦ OFITH 

WRITS (6,723) OF L $ 

R * c .♦ OfcrLR 
C. = C ♦ Of LC 
SN = SO *■ 0F15N. 

6 R 1 = F-QF L r 

R=«-0F1.3 
TH=TH-OFLTH 
SN1 = 5N-0F L SN 

ci=c-oetc 

T HI = T H-pF l TH 

SLOP: =( PFl R*f COSl TH|1 -o«=OELTH*OSIN( THI ) ) 

SLOPE =Sl n P c / 1 R *0f IT H*OCOS ( TH I )*DFLPPOSIN( THI I I 
WRITE (6,175) Pl.THJ 

17* FORMAT (* Cl = ' , F12.6, 53, *7H!= 1 ,F1 2»6» 

A6f =()F Lfi *OCOS ( Th 1 ) 

ARF*R*Of I.TH*nS!NI THI ) 

ARf,»Q*i)t L TH* OF OS ( TH I I 
A PH = 0 ELP*OSI M(TH1 ) 

ABI *A‘»e-AFF 
AOj = Artr,*AhM 

WRlTf (6.1731 A5f ,iY>F, APGc AHH.AM ,AttJ 
OFl«P= (O-c i 

Ofl TH°= OF l TH* ( D C L R°/O t L R I 
FP=kl *0f LRP/2. 

OFl $o = RP<-°°*r <co, THI ) * OFL T HP/R-N 
C»M.'P,TM1 ) 

Otl S'! s *. I M ( RP , Th l )-1.70^)*OElSP 
THP * T ♦ * I ♦ Of lthp 

HEl CP*Of LSP-(0S3RT(f'H SP**2-0t LQp**2)*DC05»(TrlP) *0 t L KP*OS I N I T HP M 
£P* f)A T AN I “5L0°F ) 

KCM=< PO*(' , COSI THP) -SI f«P t *OSrN(THP l ) ) F RHP 

THp-OAKSlM (nSO«T ( l.no*SL0P e **2-K 0N**2 1 -K ON* SLOPc I / ( 1 . L 0PE**2 I ) 
OLP=(*i)P*OS I N(THf J-fin^f^tf^THP) )*m.D0/0CUS(EP))-i.wi>0) 

WFlir (6.150) C 1 » If 1C P , 01 p 

nFLL = n-3'i.*(r i *n[ LCP^DIP) 

WRITE (6,1731 TMO*rp»«tO'»F , ,TH P , Dc L THP 
A 0 A * 0 S I N( THO) 

APP=nS| M THP) 

ABC = 40P* A^A-R'i* APB 
ARt>«l.DI)0/0C( SII P) -1 .CDO 
*RITF (6,173) AR6 ,Ap0, AeC , APD,R ,DfcL5 
ApF=OrnS( THI ) 

A0F=ObiN( THI ) 

WRIT? (6,173) ARf , APF.OF lrp 
173 FORMAT (6(012.5,5X1) 

153 FORMAT f * PC l * • , fM 2 . 5 , • DFl C® = • ,012 . 5. • OL P = * . J 1 1 . 3 ) 

PO=l:- )••.*( SN1*0E l $NP) 
p = FP 

IF (F.L c .').0) OtL L*?.n 
IF |«0. L«=. 0,0) or>»0.0 
IF (R.Lf.O.H F«?»Q 
TP1=THO**7.29*70 
T0? = THr»*57. 20578 
FD1=PSI*57.2F‘ ; ‘ , 8 
F.n?*F*57.2 Q 57B 

ACL=7*R0 O: *0SPJ(0. 5*(TH0*TH01) 

Ht'«ROP*DCOS<»'.5*( THO+THDl ) - rc 
PHOF T=OF L L*RP 
H=H-3 c. 

OELH=H-HO 

ET« IFfPSI)*) If 0. 

I =K ♦ 1 
C 

C SAVE ANO WRITE RESULTS 
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32 £ '« 
3?7. 
? 28 * 
329. 
3V-. 
321 . 
33?. 
33,j. 
? 34 . 
335. 
3 36. 
337. 
3’8. 
3 39. 
34 

341 . 
34? . 
343. 
34*,. 
345 , 
?46. 
34?. 
34*. 
340. 
350 . 
351 o 

352. 

353. 
3 64. 
364 c 

?56. 

357. 

358. 

354. 
36C . 
361 . 
36’. 

363. 

364. 
3 66. 

366. 

367. 
36P. 

364. 
37.1. 

371. 

372. 

373. 

374. 

375. 

376. 

377. 

378. 
3 73. 

380. 

381. 
3 8?. 
363. 
324. 

365. 

366. 

367. 
366. 
3ti9. 
39U . 

391. 

392. 

393. 

394. 

395. 

396. 

397. 
393. 
399. 


B< 1 l^HO^T 

RD? (I »*ACL 

D£LR=aCt-ACLl 

ABC =ACL* DFl L / 1 9CG . * 

IF MCL 1 • FQ. 01 GO 10 24 
OFl S= ( A6C ♦ ABC 1 ) /? 

DHIU* “1 0*01 0G1 Of TABS (OttRl/MELb^DPS 1 ) ) 

WRIT? <6,1021 1,06(11 
24 CONTINUE 

ACL 1 = ACL 
ABC1 = ARC 

WRIT? (6,100) 1 c L Y , PS I 

WPITt (6,113) TDl,102,£CL,f.)ELL,RO,PHi3FT,H,Hn,utLH ( LOi,fcU?i>cT, THC 
IF 1 .MOT. PL0TP1 GO TO 3 
Sf PU ♦K*1»*H 
SFP(?,K4l 1 = ACL 
3 CONTINUE 

P$ I =P S 1 40°S I 
TM2 =>ThO 

K = K 4 1 

C 

C REPEAT RAVTnAClMG A T NEW STARTING ANGLE PSI 
C 

IF (KoLTofJHM} GO TO 4 
IF ( .NOT . PLOTt ) GO Tn 7 
TH = 9. 97044/6506 
P=6410,0 
T=R*0S1N(TH)/P! 

P*R *OC.IJ S ( THJ-PF 
CALL PLOT (!,T,P1 
7 CONTINUE 

r 

f SEARCH FOR mul t 1PATH 
C 

till PHASO" (P0?tD6,6,NtJM| 

IF ( .NOT .PLOT PI Gl> TO 11 
CALL RAMP 
<S=*9.R/|?U 
VS = B.C;/5 

CALL SCALE ( XS, VS ♦7990. 0,-1. 2 5) 

CALL CHAR (750v. 0,-1 .9.3.10,0.10,01 
WRITS (14,101) IFLY 
CALL PLOT U,*»«OO.t)O.N) 

CALL PLOT (2,7900,0,5.'') 

CALL PLOT ( l, S:P( ’.NOW) ,SFP( l ,NUM)) 

OO Q | = i , N 1JM 
K * NU 1 ♦ 1 

9 CALL PLOT ( 2 , S 5 P ( 2 , K ) , SEP ( I , K ) t 

CALL PLOT (1,6)40.0,9.05 
11 CCNTINUF 

IF ( P LOT PH ) CALL PHPLOT (RDZ,R,NUM) 

IF (PLDTnB) CA1L 06PL0T ( RO? , D 6 »MU M ) 

GO TO 1 

10) FORMAT ( 0 OF L I GHT Nn a 0 ,A4 , 0 STARTING ANGL E *'»F 1 u o 6 , 8 RAD. * ) 

U-l FORMAT (A45 

1J? FOPMAT (!5,F15.6) 

1C3 FORMAT I® ft LESS THAN P E*//® R=°,F12„6 / 0 NEW MUM**, 16) 

110 FOPMATCOCFNTBAL ANGLE 1«»F19.5,« 0fcG.*,9X, * CENT R AL ANGLE ?=*,F10 
«.5,» OEG. • , 1 (‘X, * ST • LINF PATH** ,F 12, 5 , * KM.*,/* BEN&1 NG= 1 , F 9 . 3 » • 
*M.» ,19X,*RFTAP0ATI0N«»,F9.3,» M. • , 1 7X , ‘PHASE utEtCT«* ,FV. 3, * M. » , / 
#* ACTUAL RAY Mf IGH7 = • , F 8. 3 , • KM. • , 10X , • A® PARENT KAY Me i GMT* • , F8. 3 , 
#♦ KM.*, °X » • DEL T A H* * , F9.it 1 KM.*,/» BENDING AN GL t i* • ♦ H l . 8, • DEG® 
#',9X, »3EN0lNG ANGl£2=*,FH.ft,« D6G. 10* *• TOTAL bcNGiNC ANGlt=*,Fl 
#1.5,* MR. %/• ANGULAR EPROR= * t F8»4 5 6 MR.*/) 


160 

FORMAT 

(12,2131 

530 

FORMAT 

(5FI i.fS) 

501 

FORMAT 
*• PS! = 

( *1$S=» RO c0 , FIB. 8o ‘ THPREC=* ,E15o3, 

* , F 15. 8, * DFLTH=%F15oR//5 

502 

FOPMAT 

(• ITERATION NO. 0 , 14,* U=°,F15 0 8,* THl=»,Fi5.d) 

503 

FORMAT 

15F1 5.61 

720 

FORMAT 

( ■C0ELS»** ft 20.13J 

73c 

C 

FORMAT 

END 

t* 0N0*» = r, » UDO, RoO,F15 e Q s 9 TH* 0 , F 1 5 „ 8 ) 
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409 . 

‘t’l. 

4t ?. 
4 ''''i. 
4f 4. 
4f'5. 
.406. 
4t 7. 
4 08. 
4f<*. 
41 

Ml, 

412. 

413. 

414. 

41*,. 

416. 

417. 

<.ia. 

413. 
420. 
421 . 

422. 

423. 
4?4. 
425. 
4 ?6. 
4?7. 
4? 3 . 
420. 
4 3.'. 
4 ? 1 . 
432. 
4 3 3. 
434, 
475. 
4?fc. 
4 7 7. 
43^. 
434, 
44 '. 
44J , 
<42. 
4< 

444. 
444 . 
4 4 (■ . 
44 7. 

44 4. 
440. 
4 50. 
451. 
45?. 
453* 
4«4. 
45*.. 
4 56. 

45 7, 
458. 
45«. 
46'-. 
4M. 
4/ ?. 
46? • 

464. 

465. 

46 fc. 
467. 
46 8. 
4< P. 
4 70. 

471. 

472. 
47 , » 


C* 4* ***** 

c 


StlP.RnuTINF K’R F AO 

THIS SUIVIOUTIN* INPUTS DATA TO THF PPOGR4M AND INITIAL l Z L i> SOME 
VARIABLES 

REAL *3 RROI S 

RFAl Ml U2S>,H(3),HP(?1 

REAL R I. ( 4 ) . 

I NT r GFR I H ( 3 1 

LOGICAL f, L0TR,PL n TP, PLOTPH.PL CT OB 

COMM , 1 N /FLV/lFLy.PLOT°.Pl0TP»PLllPH»PL0TDB 
COMMO N /MTRAN/'I] ,H, HP, SMJN.RDIS.P, !H, INDEX 
DATA R1/637Q.O/, pf/6371.0/ 

DAT A PL/4.3 71 .0,63 7 7.0,6375. ‘>,638 7.7/ 

READ STANDARD A TMOSR HE a F N-P30FILF, N 1 1 I > 

Hill r,IV c 5 THF tnwfS ALTITUDE AT WHICH THF VCRTlCftL RzGOtUlION OF THE 
DATA CHANGES. HMII GIVES THfc RESOLUTION ABUVc mLIITJDl Hill. 

I Hi 1 1 G1VFS TH£ INDEX NUM&eR OF ARPAY N1 ( I) UT *hICH THc hcSOLUTION 
CHANGFf . 

D FAD l s , 1 J 1 J ( H < I 1 , HR ( 11,1 Hi I 1,1=1,31 
WRITt {6,1011 IH( II ,HP { ! 1 , | Hi 1 >, I* 1 , 3) 

READ (5, 1 11, N 1(1), 1=1, 1235 

WRIT- {6,1C? 1 {'11(15,1 = 1,1251 
A =N 1 (1 71 
B=N | i 3) 


ENTRY NC V;T 

PLOTR, PLOT’, JIOTph AND "LOTOS CGNTFIL PLOTTING LUUTINtS, 
READ { 8, K 5 , fmd«2 ) °L OT R , PL QT P 
WRIT- (6,1751 0 L O 7 R , PL DTP 
READ (5,1 .61 "LOT^M, "LOTOS 
WPJTF (6.JC51 PlDTPH, PIOTD* 

SU3P UO T t h : ••■(NPFAn INPUTS A CD MANIPULATES T.H‘ HAWAII OA Ik. 
GALL '•■NRc APIA) 

(F ( • N >T, RlDTR 1 r. 1 Tn > 

CALL RAMI 

CALL SC A L F (P?. 4j,-?>. 77 >,12.05 
CAtL CHAR { -C. 364, «. l,». ,?r ,9. 27,0) 

WRIT* (14,1741 I r L Y 

no 1 ? j= 1 , ? 

THT=-7.(7j 

THTsIHT + U.O'H 

RT = Rl ( J1 

T = P T* SIM ThT t /R 1 

P*B T*r.')S< THT )-PF 

IF ( P.GT. 1?.P.DR.P.LT 9.D) GO TO 17 
TALL PL0T(1,T,P| 

DO 1? 1=1 ,147 
T HT = T HT 40 « rt :> 1 
T«PT*.S!N< 7HM/RI 
PsCTYFOSi THT1-PE 

IF < D .GT.1?.C«7R.°.LT • - 3. Cl GO TO 12 
CALL PLOT ( ? , T , P ) 

continue 

CONTI NUE 


c NT R Y f)S c T 
RMIN=H(3) 

EDI S =HR ( 7 ) 

INDEX* I H I 3) 

RETURN 
CONTINUE 
PFW1N7 « n 
STOP 

FORMAT 1 8 F10 * 6 1 
FORMAT (? FI'/. 6, 151 
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47 <*. 

<•75. 

476. 

477. 

478. 
47Q. 
48'). 
4 81 . 
462. 
4B3. 

4g4, 

468. 
466. 
487. 
486. 
4R0 e 
49 m. 
401 e 

49?„ 

483. 

484. 

406. 
4S6 o 

407. 

4 9 R , 
490. 
5C'i‘. 
5(1. 
5(2. 
5C3. 
5o4 ■ 
5C5. 

506. 

507. 

508. 
5v8. 

5 1 . 

511. 

512. 

51 i. 

514. 

515. 

516. 

517. 

518. 
517. 

52 ). 
521. 
c 22. 

523. 

524. 
5 ? 5 , 
5?4. 

527. 

528. 
5?o. 
530, 
531 . 
632. 

533 . 

534. 
5 3 5. 
5 36. 

537. 

538. 

539. 
5^0 . 
541. 
54?. 

543. 

544. 

545. 
548, 
547, 


1.-2 FORMAT I s »,<JFl%fc) 

154 FORMAT (‘FLIGHT ‘,A4| 

lu 5 FORMAT |?ini 

F NO 
C 
C 

C 

C 

RFAL FUNCTION N*fl(R,TH) 

r 

C THIS SUP^nUTtM? CO^PUTFS THF INDEX OF REFRACTION F ux T tic. RlGION 
C ABOVE 4 KM AND PFLOw 2 KM 
C 

REAL* ft P ,TH,0A,«AP,RROIS.FP/6373.0D0/,NN 
RFAL N?( 1251 ,HI 31 .HRI’l ,RE/6371oO/ 

INTEGER I HI 3 1 

COMMON /MTRAN/N?,H,HR ,R Ml N c RD l S, 5 , IH» INDEX 
PA-R-P? 

J=1 

DO 2 l=lg? 

11*4-1 

IF JR Ad ToHlI I) t GD TO 2 
J=1 I 
GO TO 7 

2 CONTINUE 

3 CONTINUE 
RMlNsHl J ) 

«OIS*HF ( J I 
INOEX = IHJ J) 

1 CONTINUE 

I * 1 01 NT I <RA-PMtW)/«ni $)♦ INDEX 
RAPaROIS*! I - I NOf X I f e M | N 

J*IM 

N = N?(I | *1 RA-PAP1 *<N?( JI-N2I t ll/ROIS 
C 

C IF ALTITUDE L K $S Than 2 km, OFFSFT N-VALUE WITH V~LJfc AT 2 KM 
C 

IF 1FA.GT.2.01 GO TO 4 
N * N - 8 I NN ( ft B f T H I - 1. ^0 01*1.06 

4 CONTINUF 

N*1 #0 ">0 *N *1 .0- 4 
UPTURN 
PND 
C 

c 

c 

c 

REAL FUNCTION NN*s(R,TH) 

C 

C THIS S UBR JUT I N F COMPUTES THE INDEX OF REFftAt T 1 C).« EUa THl KtGinN 

C BETWEEN ? AND 4 KM 

C 

REAL *8 R 9 TH.N21 , M?? ,THO I S , R D i S , RH I N, N 

RFAL-»3 GA,RAP p C5LF 

REAL N2 (40,701 *RR IS/Co n 5/ 

INTEGER PM T, TIMES 
COMMON T HD I $ i TH M I N » N2 

THAoTH-THMT N 
RA«R-.8373.» 

I»1DINT( RA/pR I S ) ♦ 1 
PAP=RR1S*I 1-1 I 
J«I*1 

I «M J N7( 1,401 
J*MINO( J.4C 1 
KalDlNTITHA/THDISlo-l 
K=MAXD(k,] j 

K«MINO(K f 201 

N2i=N2U,K} + (PA-FAP|*(N2J J(.R}-M2( I,K1 1/RPIS 

K1 b K ♦ 1 

K1 = Mf NO I K 1 , 23 ) 

N22=N2I I,Kll+(RA-R4Pl*(N2(J,Kn=N2II*Klll/RRIS 
N=N21 + I N??-N2I|M( THA-OC-llfTHniSI/THOIS 
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5*8. 

5 49. 
55 /. 
551. 
55?. 

55- >„ 
f- 5«* . 

555. 

556. 

557. 
556. 
559. 

56- .. 
561. 
5 62. 
5 53. 

564. 

565. 

566. 
5(7. 
568 . 
569. 

5 7M. 

571. 
67?. 
5 7 3. 
574. 
S7F. 
57*. 

577. 

578. 
57?. 
5 ft'.). 
5Pl • 
65?. 
6ft?. 
5P4. 
565 . 
586. 
5ft?. 
5ftS. 
659. 

5 90. 
691. 
592. 

59?. 

594. 

595. 

596. 

597. 

6 Oft. 

5 99. 
6f . 
6\.-l . 
66 2. 
6j5. 
6'J4. 
6-. 5. 
*>■.£’• 
6' 7. 
6uft. 
t'-9. 
61 

£11. 

612. 

6 1 ). 
6 » 4 . 
615. 
616 . 

617, 

618. 
619. 
(■ 2' . 
621. 


IOC FORMAT (3020.6) 

NN=1.000'N*l.D-6 
UPTURN 
E MO 
C. 

C 

C********************************************** ********** j»**«s** *«#***#* ********* 

r, 

C 

SUBROUTINE NNRFADINBASf » 

C 

C THIS SUBROUTINE INPUTS AMO MANIPULATES HAWAII DATA. IHe HAWAII 
C DATA IS READ FROM TAPE, WRITTEN ON DISK, AND ACCfcSStU FROM DISK 
C WHEN K c Qlj I R( D. UNIT NO. HD REFERS TO 9-TRACK MAG T«Pfc„ UNIT 
C NOS. ?0 AND 21 PEFfP TO DISK (CALLED "RAO** IN SIGMA 5 TchMl NOLO&Vj 
C 

PEAl*3 ROIS.THOIS,RMIN,TH 
REAL Nl( 4i!»?0) ,NKASF 
PEAL N(4.3,?0) 
l or, I CAL PRINT 

EQU I VALENCE (NI ,*') » I LRANGE »NUM> 

COMMON T HD I S , T HM IN »N1 

COMMON /FLY/ J FLY 
WPITC (I , l^Q | 

WRIT;' (6.10B1 N3ASE 

C IFIY IDENTIFIES FLIGHT NUMBER. MRANCE IS THE HORIZONTAL NUMBER OF 
C. DATA POINTS I* THF HAWAII OAT A ARRAY. LRAN&t IS THt TOIttL NUMBER 
C OF HORIZONTAL DATA POINTS. 

8 HE *0 (ftO,LC7,FNO=ft| IFLY 

WP 1 T F 16,107) IFLY 
R F A0< Sj » 1 00 ) MOAMGF 
WP l T c (6,100) MPANGE 
R F A D I S v> » 1 UD ) l R ANGF 
WRITF 16,100) L RANGE 
R0IS*.)5 

TH01S=l50/( 53 7V»(*M RANGE ) 

T HM I N = -TMO | S*LR ANG5/2 

WRITE (6,103) , 

RE AO M '11 01 P»!NT 

ft F 60(30, 101 ) I IMt I, J), 1=1, 40) fj-l,20) 

I F ( P-iINTJ WPITf ( 6 * 1 0 1 ) ( CNU, ,49I,J«*, 20) 

DO 3 1=1,40 

no 3 j*i , pm 

3 NU,J)*N( I , JI+NAASE 

C RADREW, RADWRI AND RAOPO REWIND, WRJTc ON AND fttAD FROM, RESPECTIVELY, 

C THE DISK 

CAIL ft A D ft T W ( 2 0 ) 

CALL 9A0WRT ( ?'' f K # 9DO,0.65C.r;,fi50I) 

I NOr X = ( NU'6-?- ) / ) 7 

If (1 U'Hx.LT.l I GO TO 1 2 

CALL PADPfWI?!) 

OH 11 &=l ,!NnFX 
00 21 1=1,49 
N(f ,1 I* N (If 18) 

N < I , ? I S N I I , 19 J 
21 N( 1 ,a)*N(I,?0) 

ft) AD( 30,1 Cl ) H MI ,J) ,1 *1 ,40) » J*4 ,ZJ) 

IF (PR I N T I WPITf 1 6 » 1 ( I I ((N( l*J),I = l,4j),J*t,2w) 

00 5 1*1 ,40 
O') 8 J=4,29 

5 MI ,J)=NII,J)*N3A?E 

11 CALL RAOWRT 121 ,N,6«? ),0,t 5.02, C503) 

12 KRf P = NIJM-(lNDFK*17*2DI 
IF (KPEM.EO.O) RETURN 
DO h 1=1, 40 
Nn.n»Mii,iP) 

N( I,?1=N( 1,19) 

4 M I ,3 )=NU ,20) 

KRf MaKPf H* 3 

RFAO( SO ,1 M ) ( C M ( I,J), 1 = ) ,40)tJ=y,KREM) 

IF (PRINT) WPITF (6,191) ( INI I , J» ,1*1 ,4?) • J = 4,K*i.M) 

KKK=KPfM+l 
DO 7 1=1,4? 

DO 7 J=KKK ,20 
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622. 
623. 
b?4. 
625. 
626 , 
62 7. 
6 29. 
6?^. 
f V'. 
6 31 . 
<32. 
633. 
6 36. 
6 35. 

636. 

637. 

638. 
6*9. 
64". 
641 . 
642. 
<4*. 
6 A*.. 

665. 

666. 
667. 
663 . 
66 6 . 
< 60. 
661 . 
65?. 
663. 
6 5'*. 
655. 
6 66. 
667. 
668 • 
65<=. 
6 6 r. 
661. 
66 ?. 
663. 
< 66 . 
666 , 
666 • 
667. 
669, 
669. 
6 7r . 
671. 
67? . 
673. 
67** . 
673. 
67fc a 
677. 
6 7«. 
67!?. 
66 0. 
681 . 
682. 
633. 
686 , 
6«5. 
o86. 
667 . 

6R6. 

689. 
69 >. 
691, 
6 C 2 • 
6 e 3 . 
696. 


7 N< I * JMNRASE 

no 6 [=1,90 

00 6 Js6,K,RFy 

6 M( I • J I *Nf I « J I 

CALI RAlwET .<21 . N, 6 »:>♦!,* fi5 rt 4* 6505) 
Rt T1J9 M 


5 >0 

MN= l 
GO T n 

4 O'*' 


501 

NN = 1 
GO TO 

60 1 



NN=2 
GO TO 

OtO 


50 3 

NN= 2 
GO TO 

601 


5 j6 

*T> Z 
iF 

— 1 ^3 

L3 

60 o 


505 

*N = 3 
GO TO 

iul 


6 )U 

WRITE 

STOP 

16 »?0P 1 

NN 

6M 

wf i r~ 

<6*231 > 

NN 


STno 

C 

C 

ENTRY NSTART ( TH c TH1 *'•40 X 5 

C N S T ART PEApS THf INITIAL . S F T OF N- VALUES FROM ThL UlSrf, iNTj CORE. 
THMl*! = -THP|5«-LRANf.E / ? 

THA = 1 7*T HOIS 
THl NOX-= THMJ N*T*-*4 
IF !TH 0 r,F o THI\nx> GO TO 1 
CALL W>fW <?*H 

CALL -1 AQRQ f ?0*N, 8uO, 0, C 41 ./C * 1 40 II 
CALL °4r>R£W <?11 
GO TO 13 

1 CALL *AD«*trf <21 I 

INC*< TH-THMINl/THA 
00 2 1>1« I NO 

? CALL RAPED ( ? 1 , N » 6PC *0 , £400* £4311 

TMINPX=THINDX»TMA*INP 
THWlN=THlN0X-T«6 
GO TO 13 

40 J WRITE <6*3301 

RETURN 

4C1 WRITF <6*331 1 
RETURN 
C 

c 


C NTRY NSHIETITHINOXl 

C NSHIFT RF AOS Thf N^xf SET OF N-VALUFS FROM THt DlSrt INTO CURE 
THMIN=THIN 0 X 
TH 1 NDX = 7 HIN 0 X «-THDIS*i .7 
CALL RAORD I ? 1 * U , t 80 ,0 * £ 506 , £ 50 ? > 


13 

CONTINUE 


5 36 

RETURN 
WE ITE 

< 6*2021 


5<-7 

RETURN 
WR | TE 

< 6*2031 


I UO 

RETURN 

FORMAT 

< 141 


till 

FORMAT 

«?<16F5.1/I*8F5.1I 


lu 2 

FORMAT 

< I5,5X,fi >,61 


103 

FORMAT 

C ■ 1* , 10V, MNPUT )ATA'//» 

106 

FORMAT 

I F 1 0 . 6 1 


106 

FORMAT 

<*1'*4P<* *, 5 1 FI 0. 6 

*5X>/)1 

1*7 

FfJRMAT 

( A41 


109 

FORMAT 

1 « NF A SE * * *812.61 


U 9 

FORMAT 

M PRINT? FORMAT L2 

e 1 

no 

FORMAT 

<L?1 


2 10 

FORMAT 

< ♦ ERROR IN f 4nWR T 

STATEMENT 

201 

FORMAT 

<« EOF c ivjc HUNT cR FD 

IN RAOWPT 

2 u2 

FORMAT 

( ! ERROR IN RADRO 0 1 


2 j3 

FORMAT 

<« EOF ENCOUNTERED 

IN RADR 0° 1 

30 (> 

FORMAT 

< • ERROR IN INITI AL 

P AQR Q° 1 


NU. * , 1^1 
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*.<55. 

69ft. 

6 SIT. 
tse. 

699. 

700. 
701# 

7 r-2. 
? r -3. 

704. 

705. 
7<6. 

707. 

708. 
7i 7. 
720, 

711. 

712. 
7)3. 

714. 

715. 
7)5. 
7)7. 

71 A. 
719. 
72?. 
72) . 
77? . 

72 3. 
7?4. 
72 5, 
72fc « 
727. 
72H. 
729. 

736. 
7 3T . 
73?. 
733. 
7?4. 
7’6. 
7 3ft. 

737. 
7 ?* « 

7?4. 

74 

741. 

742. 

743. 

744. 
74 5. 
74ft. 
747. 
74a. 
7*P. 
7 5.‘. 
751. 
75?. 
753. 

.754. 
755. 
75o. 
7 6 7. 
755. 
759. 
750. 

761. 

762. 

763. 

764. 

765. 
76fa. 
7f 7. 
76*. 


301 FORMAT <• EOF 9NC0UNT FP FD. IN INITIAL RADPD*) 

END ' " ; " . ' 1 

f 

f. 

c 

c 

SUBROUTINE PH A SOh <H, 09, PH, IT c-R) 

C 

C THIS SlJHR JUT IN C INTERPOLATES RECtIVED PH A St UcFcCT ANU POWcR AT 
C l KM If'jTEPV&LS IN $ AT C LL 1 T E SEPARATION. IF A MULTIPATH Ii DETECTED. 
C THE SEPARATE AND COMBINED VALUES AT EACH 1 KM INTERVAL AKc PRINTED. 

C 

DIMENSION HU TER ) iPBM T 5R) » PH < I Tfc R ) * PHD 1 1 A ) ,0t>0 II 4 ) 

CCMPl CX E * E T 

«fU 7W0P I/ft. 2932 35/ 

CCMMC N /FLY7 I Ft Y 
PHASE(P)=4MOnio,o.r?2)/0.032*TWOPI 
DEI 1 )=Ju( 21 
DP l't 1=1 ,ITFO 

l) writ? if, i?4) i,Hin,o>»in,PHin 

HM I N= 1 . E 5 ) 

UMAX*- 1.F.3 

on i 1 = 1 , it c r 

HH=H< I ) 

IF UM.LE.HMIN) HM1N=HH 

1 IF (HH.GF.H*AX) H ** A X bMH 
GH= A I M f I HM AX ) 

WRIT 6 (6,1 ?5) HMJN.HMAX ,0H 
WRITE <6,InH I F t Y 
WRITE 16,10?) 

WRITE IftU'U) 

7 HH>H1 1 ) 

VN*7 

ET= >,'? 

00 4 )=?, J T C R 
MMaMC I P 

IF <HH. IT. MW) c, 0 T° 5 

IF ( GH.GT.MH. nRiGH.LT. HM) GO TO 6 

8 CCNT 1 MU C 
NN=N'Jf 1 

PP1 = »JH( I -1 ) 

PP?=PH( | ) 

DM J = 0 3 ( I - 1 I 
DB2 =f'P ( I ) 

f AC = K H-) H)/<HM-HH) 

0*i'.)lNNl*0<U*lri?-PP1 ) *FAC 
vot T = 1 3 *»I')R 0 {NM) /20) 

PI=PHAS7 <o°l P 
P2»PH.».Sf I pp? P 

e?-°l |*CAC 

E »C MP LX ( V r LT * f.nS ( P) ,Vt il T*S IM I P ) ) 

tTjf) #c 

PHD INN) = PP1MPP2-PP1)*FAC 
ft )<H«H ■* 

GO Tn 6 

5 IT IGW.t T.HH. MR. e.H.Cf.HH) GO TO 6 

r.n t.j r 

4 CONTINUE 

IF <r AOSIETI.LS.O) GO TO. 3 , 

VOLT = 2 •* 3LPG1 01 P ArtSIfT ) ) 

GO Ti 9 

3 v^LT=l.r-« 

* CONTINUE 

P = 4 T 37)? ( 4 I WAG If T ), P c At ( cT) ) 

P*P* 1 A'.'/ 3.1 615 S?7 
IF (NN.LE.lt GO TD 2 

WRITE <6,1011 GH,VfUT,P,NN,<PHDUP,I=l.NN) 

wmt c (6, u-ft) npnm,i*i,NN) 

2 GH=GH-1 

IF IGH.GF.HMIN) GO TP 7 
RETURN 

IOC FORMAT I*l , .22X,'PFCEIVfO SlGNAL'»5X»AA/7) 

101 FORMAT {F6.1,5X,FU»5,4X,F6.2 t l9»5<,7ElU.3Z4&A,7FlW.32 
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769. 
77u. 
771 . 
772. 
??i. 
774. 
7 7?. 

776. 

777. 
776. 

7 7 C • 
78«». 
7B1. 
7??. 
7fx. 
764. 
785. 
7P6. 
7t 7 • 
7BR. 
769. 

79 >. 
701. 
79?. 
7»3. 

794. 

795. 

796. 

707 . 

795. 

?oQ. 

80 3. 
80 1 . 
80?. 
80 3. 
904. 
805. 
Hf 6, 
80 7. 
90 P. 
809. 

6 1 :: . 
811. 
812. 

813. 

814. 
816. 
916. 
617. 
616. 
619. 
820. 
821. 
822. 
e23. 

624. 

825. 

826. 
827. 
328. 
829. 

8 3,'. 
631. 

832. 

833. 

834. 
935 . 
836. 
637. 
838. 
8 39. 
640. 
841. 
84?. 


102 

FORMAT C SAT SEP PEL S1G AM P ANCLE nO. 

* 1 0 X » 1 PHASE DETECT (ABOVE) AND POWER (BELOW)') 

JF RAYS' , 

133 

FORMAT 

(• (KM) • ,7X, MOB) «.5X,Mi>FG) S35X, 

• IMETcRS ANO 

1-J4 

FORMAT 

( I 5,3F1 5,4,1 


H-5 

FORMAT 

(i HMJN*‘,F1^.4, • HMAX=« ,F10.4, ' GH=’(,Flo 

o 4 ) 

lu6 

FORMAT 

( 2 (48 X. 7 FI 0.3) ) 


H‘7 

FORMAT 
f NO 

( * ' » 



C 

r 

C * *#** * ************ <■«■*$ ********* as ***«* *<>*#* **«><t**4t*i>Oi*fr* *»$** *$***« <-.11 ****** 

0 

c 

SUB POUT l M c PHD^pT (Rn,B,J) 

C 

C THIS SUB»ni)Tli<f PLOTS PHASE OFFfcCT AS A FUNCTION OF SATtlLlTE 
C SEPARATION 
C 

PlMfiMSlbM *0<J) t H(Jt 
COMMON /PLY/ I FLY 
x:>= 8 . •.•/ 1 400 

Y C s 9 , / 1 50 

YY= 790joO« 805' p -. ) 

CALL P4«L 

CAU SCALE ( »0,YC,10J. 0,8050.0) 

CALI CHAP (30 A . 1, 8040.0*. 2». 2,0) 

WRIT* (14,104,) 1FLY 

CALL CHAP (1O0.O, 8049. C l v .t,OI 

WRITE (14,1(31) 

call grid i 3, 200 . j,eo5o<>u<> in o 0»i5> 

CALL CHAR (IPO. 2, 7898.0,,!, ,1 ,0) 

WRITE (14,1001 

CALL CHAR U5"/' ,7896. C,.l, .1,0) 

WPITF (14,10?) 

CALL GRIP (0, 2JJ.P, T90P.H, 109.0,1 4) 

CALL CHAR (1550. 0, 7895. 0,.l ,.1,0) 

WRITE (14,1031 

Y = YY-RO< 1 I 

IF (Y.lT. 7900.0) Y=7«0n.<) 

IF tY.GT.R050. 3) Y»R05C.O 
CAL L PLOT ( 1 ,8(1 1 ,Y) 

CALL ®0 1 NT (9) 

no i i=?,j 

r=YV-co( r i 

IF (Y.LT. 7900. 01 Y*79no 0 9 
IF (Y.GT. 8050.0) Y-P050.3 
CALL PLOT 10, BI I),Y) 

CALL P01NTI2) 

1 CONTI MUF 

X = 1 600,0*2. P/XO 
CALL PLOT U » X, 8050. J) 

RETURN 


1UP 

FORMAT 

( '6050 M 

101 

FORMAT 

( • 79C0 ' ) 

1 C 2 

FORMAT 

I *200* ) 

133 

FORMAT 

( • 1530' ) 

104 

FORMAT 

END 

( A10 ,5X, •PHASE' ) 


C 

c 

tt«»#(»*,H^#<.»»t»<.*»6»m*t**»*t*t*******t»»M**»»***MI*<iiM*»l**MMM<t*** 

C 

c 

SUBROUTINE OBPLOT (RQ,OB,J) 

C 

C THIS SUBROUTINE PLOTS PEf.FIVEO POWER AS A FUNCTION JF SATELLITE 
C SEPARATION 
C 

DIMENSION roi ji ,ob< js 
enHMON /FLY/ I FLY 
X0=B.O/30 
Y0«9.9/15( 

Y Y= 79 TO. 0*805p, T 
CALL RAML 
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64 t. 

344. 

349 . 

343. 
047. 
349. 

349. 

350. 
..35) . 

P5?. 

. a«3. 
B 54 . 
355. 
856. 
657. 
853. 
859, 
860* 
861 . 
862 • 
86 . 1 . 

864. 

865. 
8 66 . 
367. 
868 . 

869. 

870 . 
67U 
97?. 
P71. 

6 74 . 

876. 
6 7 • 
677. 

878. 

879. 
6 80. 
88 ). 
8 92. 

8 83. 
I 84, 
886 , 
RP 6 , 
K87. 
339. 
8P9. 
390, 
991. 
9 °? » 
057. 

344, 
896. 

9 96 . 
697. 
89?. 
399. 
903. 
9(1 . 

902. 

903. 

904. 

905. 
506 . 
407 . 
908. 
90°. 
41 J. 
Ml . 
M2. 

913. 

914. 
615. 

916. 

917. 


CALL SCALF (KC.Yd.-30. 0,8053. 01 

CALL CHAR |-25. 9, 8040. C‘.. 2,.?. 01 

WRITF ( 1 4 .104) IFLY 

CALL CHAR I -32. 0, 8048. O, . ! , . 1 , 0 ) 

WRITF (14,171). - 

CALL GRIG (3,-10.0,8050.0,10.0,15) 

CALL CHAR (-32. G, 7898. 0,.L, .1,01 

WRITS 114,100) 

CALL CHAP (-30. 5,7895,0. .1,.!, 01 
WRITF (14,1021 

C ALL GRID ( 0, -30. 0,7900 .0 , 5.0,61 
CA(L CHAR l-n.5,7895.0,.1,,1,0) 
W»ITF (14,193) 

Y=YY- 5 0( , 1 

If (Y.LT.790C.0) Y = 7900 . 0 
IF (Y.GT.805G.01 Y= 87 50.0 
CALL PLOT (1,03(2), V) 

CALL P'lINT <<M 
OP 1 1=3, J 
Y=Y Y-R 0 ( I I 

IF (Y ,LT. 7000.0 I V = 79 OC .0 
IF (Y.GT. 8050.(1 1 Y = 69 50 . (J 
CALL °L HT (:),0«m,Yl 
CALL POINT (?) 

1 CUNTINDF 


VjO 

iul 

10 ? 

103 

U.4 

C 

C 


X=2. )/xz 

CALL PLOT ( I , X,R050.<7) 
UPTURN 

FORMAT PB05V) 

FORMAT ('7900M 
FORMAT (•-3'' B ) 

FORMAT (• DM 

FORMAT ( A17.5X, ’POteER* » 

f NO 


c 

c 


JtOAO (GO) , I MAP, PROG), (UOCH, 21 
: INC. LOO; |P:XYPLT) 

: ASSIGN (FS?-.8T,<1 ) 

JASSIGN (F:21,RT,K?) 

ALLneT (FILE.K1 I ,(FSI7E,1 I , (RSI7 C ,B00I 
AUvJBf (FILL.X2), (FSI2F.46J, (RSIZ8,660» 
ROV 


69.0 


1.0-7 


7925.010 .(07 


3" 




J.JoJ "* 

.250 

1 


ll,i-07 0 

.501 9 

45 


32 .OOC 1 

.007 

87 


o.OwCn 

341.3928 

7.2570 

327.3831 

1 ,000) 

2 c ?.04C6 

1 . 2 500 

282.0127 

2. COOO 

253.15L « 

2.2509 

244.3950 

3.0099 

2l7.p?44 

3. 2 5 JO 

212.6913 

4. ('790 

1 9? o 6675 

4.269,9 

1 96.7897 

5,nf>i>7 

l7-.).?47o 

5.2500 

165,0793 

6.0 GOO 

151 . 3403 

6.2890 

146.1870 

7.00( 9 

1 3 1,4315 

7.2500 

129.4621 

P.L'Dl-J 

1 18.1977 

3.2500 

113.5425 

9.0( or* 

104,4268 

9.2509 

101.5135 

10.(000 

43.1453 

l J.250J 

49.4742 

1 l.i.-JC "1 

8?, 3 >77 

11.5000 • 

77.9737 

13,(09.' 

6^,7463 

I3,600o 

60.7433 

15.9070 

40,fl3or 

15.5007 

46.5888 

17.(000 

36.8649 

1 7 0 5oQ0 

33.4527 

1 9.C OjO 

25.1183 

14.5700 

22.8748 

21.(070 

17.3603 

21,5000 

15.8726 

?3.( -JG9 

12,2558 

23 • 5000 

11.2736 

25. (.090 

0.7970 

25.5090 

8.1059 

27.(079 

6.3657 

27.5000 

5. 8 65 3 

29.(990 

4.6205 

29.5909 

4.2704 

31. ( COM 

3.4531 

31.5000 

3.1962 

34,0030 

2.1 P?7 

35.0000 

1.8787 


J.5090 

3 i4 • 5 35o 

0.7503 

302.4658 

1.5000 

271. 52^3 

1.7503 

262.328b 

2.6-39P 

2^4.6477 

2.7593 

226.9707 

3 o 50u0 

2u5. 505b 

3. 7503 

199,4466 

4.5000 

1 B1.2236 

4. 7500 

175.7153 

5.5 900 

16J.2175 

5. 7503 

155.4059 

6.5000 

1 41. 7899 

6. 7500 

137.6058 

7.5099 

l 25. o l5o 

7. 7539 

121.9384 

8.5JOO 

1 iw • 4 io? 

8.75Q0 

107.4028 

9.5000 

90 . o 65 i 

9.7509 

95,875* 

10.5,'Ut 

07. 86i 4 

10. 7503 

85.3063 

12.990J 

79.3565 

12. 5003 

63, 9497 

14.0000 

56.9*48 

14. 5 (.00 

53.2913 

16.7030 

43.4 4* 9 

) 6, 5000 

40.4898 

13.900O 

3^.9727 

13.5009 

27.6086 

20.0y00 

2w » 8 5o 5. 

20. 5009 

19.0221 

22,9030 

14,517a 

22.5(00 

13. 3285 

24.0000 

1^.3747 

24. 5000 

9.5514 

26.0009 

7.4 71o 

26.5C09 

6. 8898 

28.7900 

5. 4140 

28. 5( 09 

5.0011 

30.9000 

2.8i5o 

30. 5000 

3.7321 

32.0 jJJ 

2.9595 

33.0009 

2.5401 

36.0000 

1.6190 

37. 9000 

1.3971 
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916. 

38. 0 jf J 

1.2^72 

39,0000 

1.044? 

4 JoOoOj 

jo9045 

41. 0000 

919. 

4 2, 10 IV 

"v.691 1 

43.0000 

0.5921 

44 , 0 O(j0 

uo 6 15 i 

45. JO 00 

QJ'l, 

46.1 U'O'J 

0.3917 

47. O' 00 

? o 3 4? 1 

48.000 i 

o.J UJi 

49 0 JO JO 

921. 

so. t r- io 

t.233! 

51. 1 li:JO 

0.2054 

52oCOOO 

Jo i 8 a 6 

53.OCO0 

92? • 

54. Ci flu 

To 14?7 

c -5 .0000 

0.1264 

56.0000 

J.i 116 

57.0000 

9 ?3. 

5 3 , f ' 9 f 1 J 

O.C«72 

59.0000 

C.076 9 

60.1/00 J 

0 w679 

6i. or jo 

9?4 . 

62. OOC •> 


63 .Tf'O.j 

J.C46 8 

fc4aJ'/-.>v- 

Jo u4l 2 

65.0007 

92?. 

66. r3')u 

3.031 7 

67. 00 r ^ 

0.0278 

fad.Ou'JJ. 

J.J243 

69. DC 00 


9 2fc. ’>.•>184 


9?7. 

PLOT 

NU 

9 28. 

PLOT 

PLOT 

929. 

PLOT 

M(1 

93... 

PLOT 

PLOT 

931 . 

pi nr 

MO 

93 2. 

PLOT 

PLOT 

933. 

PLC'T 

MQ 

0)4, 

PLOT 

&I.OT 

935. 

PLOT 

MO 

936. 

PLOT 

PLOT 

937. 

PIN 



0,7845 
0.4490 
0,2646 
0,1611 
0.0988 
0.0601 
J, 0362 
0.0212 
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Appendix C 


DESCRIPTION OF RAY TRACING PROGRAMS 
1 . Description of RAYTRACE 

The technique used in RAYTRACE is based on Snell’s Law, which can 
be written as [53 


A = 

P dp 


(C.l) 


where 7\ is the wavelength in the medium and p is the radius of cur- 
vature of the path followed by the ray through the medium. Now, 



(C.2) 


where n is the index of refraction of the medium, c is the speed of 
light in a vacuum, and v is the frequency of the wave. Therefore, 


dA = - -f- dn 
n v 


and 


(C .3 ) 


A _ c 

p nvp 


c dn _ dA 

2 dp ” dp 
n v 


1 1 dn. 

p n dp 


n 

dn/dp 


(C.4) 


The total derivative dn/dp can be written as (see Fig. 59) 


dn -» 
-r- = Vn 
dp 
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■ dn | 
'dr ' 
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rt 1 dn 
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r 



(positive d0 
is clockwise) 


Fig. 59. COORDINATE SYSTEM AND RELATIONSHIP OF p 
TO dn/dr AND dn/d0. 


where r and 0 are polar coordinates with origin at the center of the 
Earth, and 0 is the angle between ~p (the radius of curvature) and r* 
(the radius vector). Note that, for r* pointed outward, dn/dr is neg- 
ative for a normal atmosphere. Now, p can be written as 


n 

P ” |dn/dr[ cos 0 - (l/r)(dn/d0 ) sin 0 


(C .5 ) 


The first term in the denominator is the effect of vertical gradients on 
the raypath; the second term is the effect of horizontal gradients and 
is usually much smaller than the first term. 

The program operates as follows. An "atmosphere 1 ' is given to the 
computer in the form of a temperature profile (temperature vs altitude) 
and surface values of total pressure and water-vapor pressure. From 
this information, the program calculates values of N (the ref ractivity ) 
for a grid representing a vertical slice of atmosphere. Grid points are 
spaced vertically 1 km apart to a height of 90 km and horizontally every 
15 mrad for a total horizontal spread of 300 mrad, centered at 0=0. 

Starting at 0=0, the program calculates the values of the de- 
rivatives dn/dr and dn/d0, determines the angle 0 between the 
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radius of curvature p and the radius vector r, and uses this infor- 
mation to calculate the magnitude of the radius of curvature p. The 
program steps off an increment d(9 of 0.5 inrad, with p held constant, 
and computes the resulting increment in r. The new values of r and 
6 are then used to calculate a new p, and the process is repeated un- 
til the ray reaches a height of approximately 65 km. At this point, a 
linear approximation is used to complete the raypath to the satellite 
altitude. The program returns to the 9=0 position and repeats its 
calculations for the other half of the raypath. A running tally of path 
length and deviation of path length from that of a straight line is main- 
tained. These tallies are used at the end of the run to calculate the 
phase defect that would result from such a raypath. The above procedure 
is repeated in each run for starting altitudes (at 9=0) of 0 to 30 km, 
at 1 km intervals . 

The equations used to determine 0 , r, S (path length) and C (path 
deviation) are given below. 

For the first step (at 0=0=0), a second -order approximation is 
necessary: 



V -^> 2 

2 [ 1 - (VPo 5 . 


(C.6) 


The program now steps to the point (r + Ar, 9+Ar) and computes dn/dr 
and dn/d Q at this point (using linear interpolation between the points 
on the N-grid). The value for 0 at this new point is computed by as- 
suming that the raypath is a smooth curve, which implies that the new 
value of p to be computed at this point will pass through the center 
of curvature for the previous value of p (see Fig. 60). This leads 
to 


K 



1 - 


n-1 


sec 0 


n-1 


A0 



(C .7) 


175 


SEL-71-05 9 



As n ~ j 



Fig. 60. RELATIONSHIP OF 0 , p, A<3, r, AND Ar. Note p n and 
p pass through the same center of curvature. 

Ut J- 


The value of p is calculated from (C.5) for the next step. The 
n+1 

incremental path length AS n is given by 


AS = r A 9 sec 0 (C.8) 

n n n 


and the incremental path deviation Ap is given by 

n 


AC = AS [1 - cos (0 - 0 )] (C.9) 

n n n n 


Finally, the new value of Ar fi+1 is calculated from 
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(C.10) 


Ar = r A0 tan 0 
n+1 n n 

and, using r n+1 = r n + ^ r n+ i and ^n + i = ^ n + A©, the program begins 
the new iteration. 

These equations are derived in the next section. 

2 . Derivation of RAVTRACE Equations 

Derivation of the first-step equation [Eq . (C.6)J. 

In the beginning, 0 ^ = 0 q = 0, and p* lies along r*: 


Ak 



Because 0^ and ^ are very small, 

A* = r 1 6 1 = (C.ll) 

and 

Ay = p i - P x COS ^ = P 1 (l - cos i^) 

Because ^ is also very small, cos ^ « 1 - (iJr^/2), therefore, 
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Now 


Ax = (p 1 - Ay ) tan \|f 



Ax = pj cos \|r tan t 1 = P x sin 
so 


Ax = 

and 



Ac 2 


Now 


, if Ay = Ax (4^/2), then from Eq . (C.ll), 


'A 

* = ^7 
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As a result 
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Equation (C.10) is derived directly from Sketch 3: 

Ar = r tan 0 „ A 6 

n n-1 ^n-1 

Equation (C.9) is derived using 



D 


B v C 



The quantity £C represents the difference between A £ and the 

n n 

length AB which can be found as follows. The angle OBE is given 

by 90°.-e n , Because OBD = 0^ and CBD = 90°, by simple arithmetic 

ABC = e -0 . Thus , 
n n 

AB = AS cos (0 - 0 ) 

n n n 


and 


£C = AS 1 - cos (0 - 0 )] 

n n[_ n > n J 


WhlChlSEq - <C - 9) - ■ - BLANK. NOT 
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Description of THRUWAY 


3 . 

THRUWAY is a program that traces a ray path between two satellites 
in Earth occultation configuration. It differs from its predecessor 
RAYTRACE in that it begins the ray path at one satellite, follows it 
through space to the atmosphere, traces the path through the atmosphere, 
and then out to the second satellite. RAYTRACE began at what was arbi- 
trarily designated the center of the raypath, at a specific altitude. 

It traced the raypath in one direction to satellite altitude, returned 
to start, and then traced the raypath in the opposite direction out to 
satellite altitude. The intersections of the satellite orbit with the 
raypath gave the satellite location that would result in a raypath with 
that specific starting altitude. By starting the ray at different alti- 
tudes, it was possible to interpolate the results to find the character- 
istics (phase defect, bending angle, minimum altitude) of a raypath for 
a given satellite separation. 

THRUWAY does the same .job, but a bit differently. In RAYTRACE, 
the independent parameter is the minimum ray altitude; in THRUWAY, the 
independent parameter is the starting angle of the ray. A ray is 
"launched" at a particular starting angle, and its path is traced to 
and through the atmosphere and out to satellite altitude again. The 
straight-line distance between the starting point and the intersection 
of the ray and the satellite orbit is designated satellite separation 
(SS). The starting angle is then stepped 1 mrad, a new path is calcu- 
lated, and the derivative d(SS)/da is calculated (where a is the 
starting angle). This derivative is then used in an iteration scheme 
to calculate the correct starting angle ; 


a 


h+l 


da 

d (SS ) 


(SS 


true 


ss ) + a 

n n 


Iteration continues until the error term SS. - Ss is acceptably 

true n J 

small . 

The advantage of THRUWAY is that it is no longer necessary to 
stipulate that the center of the raypath must coincide with the center 
of the available data set. The satellites can be positioned anywhere, 
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and the resultant ray will be traced regardless of where the minimum 
altitude of the path may fall. 

A program such as THRUWAY became desirable with the acquisition of 
the Hawaii data. With such a data set, it is unrealistic to arbitrarily 
define the center of the raypath as the center of the profile; in fact, 
the information for which we were searching might be masked entirely by 
such an assumption. Thus, a program which could start at a well-defined 
transmitter and trace its way through to a well-defined receiver was re- 
quired. THRUWAY was developed to meet this need. 

Two forms of THRUWAY exist. One version (called HAWAII) is used on 
the XDS Sigma-5 computer to raytrace mountain-to-mountain through the 
Hawaii data. In addition to raytracing, this program also computes power 
densities at the receiver distance, draws plots of index of refraction, 
power, and phase at the receiver distance, and traces the raypaths them- 
selves. The second version raytraces from satellite-to-satellite through 
the same low-resolution atmosphere as does RAYTRACE, with high-resolution 
Hawaii data superimposed at the lower altitudes. 

4. Mathematics of THRUWAY 

All of THRUWAY is divided into three parts (Fig. 61): 

(a) the section of raypath between the transmitting satellite 
and the upper edge of the atmosphere 

(b) the portion of the raypath actually in the atmosphere 

(c ) the path traversed by the ray on its way to the receiver 
after it has left the atmosphere. 

In the first and third sections of the raypath, the ray travels in a 

straight line at the speed of light in a vacuum. . The ray slows down or 

bends only in the second portion of the path. 

Because the raypath in (a) is a straight line, no real "ray tracing" 

is involved. One begins at a given position in space, say (r ,6 ), and 

s s 

launches a ray at a starting angle a which is the angle between the 

tangent to the satellite orbit at (r ,0 ) and the raypath. The prob- 

s s 

lem is then to determine at what point (r does the ray encounter 
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Fig. 61. THE THREE PORTIONS OF THRUWAY, Portions 1 and 3 are 
straight-line segments that lie outside the atmosphere. Por- 
tion 2 is within the atmosphere, where the refractivity is 
greater than zero. 

the atmosphere and what is the initial radius of curvature vector (p, 0) 
at that point. These are the initial values for the actual ray tracing 
procedure of part two of the raypath. In addition, although the ray 
does not actually bend while it is in free space, it does acquire "bend- 
ing" which is defined as the excess path length caused by the raypath 
being longer than a straight-line path between the satellites. These 
parameters are illustrated in Fig. 62, and the equations relating them 
are derived in Section 5. 

In (b), the ray is in the atmosphere, and THRUWAY, like RAYTRACE, 
uses Snell's Law to raytrace through the atmosphere. Snell's Law can 
be written as [5] 


A = dA 
P dp 

where A is the wavelength in the atmosphere and p is the radius of 
curvature of the raypath. For this portion of the program, therefore, 
THRUWAY uses essentially the same equations as does RAYTRACE. The dif- 
ference lies in how the parameters are initialized. The important pa- 
rameters for raytracing using Snell's Law are 
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r, 6 = position coordinates of the ray 
p = radius of curvature 

0 - angle between the radius of curvature vector and the 
position vector 

The origin of the coordinate system is the center of the Earth. 



Fig. 62. SIGNIFICANT PARAMETERS OF THRUWAY. 

In RAYTRACE, initializing these parameters is not difficult. Because 
the starting position of the ray is an independent parameter, r and 9 
are determined immediately. Similarly, because the ray is constrained to 
start out horizontally, the angle between the radius of curvature vector 
and the position vector is zero (0 = 0). The radius of curvature is com- 
puted immediately by use of Eq. (C.5). Although the actual initialization 
of parameters poses no problems, the first step iri the raytracing does. 

It is necessary to go to a second-order approximation for Ar because 
the first-order expression used to compute Ar along the raypath yields 
a value of zero if 0 = 0 . 
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In both RAYTRACE and THRUWAY, the independent parameter in the 
ray tracing is the coordinate 0, and A© is set to some constant 
value (usually 0.5 mrad). Because of this finite step size, THRUWAY 
does not have to use second-order approximations for Ar, unless, of 
course, the ray enters the atmosphere at a value of 0 which is an 
integer multiple of the step size; this is an extremely remote possi- 
bility. 

In (c), the raypath is the segment between the edge of the atmo- 
sphere and the receiving satellite. As in (a), the ray travels in a 
straight line, and one must determine how much bending is accumulated 
as the ray approaches satellite altitude. This segment of the program 
is identical to the corresponding segment of RAYTRACE. 


5 . Derivation of THRUWAY Equations 

The equation relating r , 6 , r , 0 , p, and 0, and the expres- 

_ . s s o o 

sion for bending in part one of the raypath, can be derived from Fig. 63, 
which is a segment of Fig. 62. 



Fig. 63. A SEGMENT 
OF FIG. 62, SHOW- 
ING THE PARAMETERS 
IN GREATER DETAIL. 
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It is useful to establish a relationship between 6 , 0 , and ijr 
(which is the angle between the raypath and the line connecting the 
satellites). From Fig. 63, 


£ CBD = e - 0 
o 


•i ABD =90° 


Now, 


so 


Because 


we have 


Also, 


or 


£ ABC - £ ABD - £ CBD 


£ ABC = 90° - 6 + 0 

o f 


5* 90° - * ABC 


>jr = - 0 


a = 9o° - (90° - e ) - y 

s Y 


(C . 12 ) 


a - 9 — >jf 

s 


(C . 13 ) 


Now, from the law of sines, 

r^ sin (90° - + ty) 

r “ sin (90° + 0) 
s 

or 


sin (90° + 0) = -5. sin (" 90 ° - (6 
r 

sin (90° + 0) = ~ cos (0 - f) 

o 
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where 


Because 


and 


we have 


or 


90° + 


-1 

sin 


r 

s 

— cos 
r 

o 



0 = 6 - i 

o 


90° + 6 

o 



9 = -90° + + sin 

o 



(C.14) 


(C .15 ) 


where r is the satellite altitude, r is the altitude of the upper 
s o 

layer of the atmosphere, a is the "bending angle" or "starting angle" 
of the ray, and ^ is given by Eq. (C.2). 
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